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ABSTRACT 

Algebraic Methods are employed for frequency response studies of 
physical systems. The transfer function of each system's mathematical 
model is composed of at least two independent variables embedded in the 
coefficients. These, together with magnitude, angular frequency, and 
phase angle are the variable parameters to be studied. 

Computer programs are developed for graphically representing these 
variable parameters. Interpretation of these results is made to evaluate 
the utility of such methods in analysis studies. 
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1. Introduction. 



The dynamic response of linear systems can be found basically by two 
methods; experimental and analytical. Continuous advancements in tech- 
nology have precluded the feasibility of the former, relegating its utility 
to hybrid studies at best. Increased attention has therefore been focused 
on the analytical techniques available, on development of new techniques 
where necessary, or possible improvements in known techniques. A funda- 
mental dilemma to the analytic approach is the synthesis of one or more 
linearly independent differential equations which adequately describe the 
system to be studied. Whether the problem consists of multiple equations 
involving eigen values, or a single mathematical statement, the system 
limitations and characteristics are revealed by a thorough study of the 
roots of the characteristic equation. This equation is generated im- 
plicitly by matrix representation of a set of simultaneous first order 
equations, or explicitly by a rearrangement of a single higher order 
polynomial. 

Many sound advances have been made since the turn of the century in 
the realm of root location techniques for time invariant systems. As 
control engineering emerged from its embrionic stage and became embedded 
more firmly in all phases of engineering, added emphasis was placed upon 
the knowledge of systems behavior in a dynamic atmosphere consisting of 
at least one varying parameter. Among the more important achievements 
in system stability analysis are the Routh, Mikhailov, and Nyquist cri- 
teria. Frequency response characteristically has been greatly utilized as 
a fundamental system specification , and the graphical techniques of 
Bode and later Nichols have become widely used throughout the world for 
both analysis and design. Many other graphical tools have been formulated. 



9 



including sophisticated root locus schemes involving cascade compensa- 
tion, and various forms of feedback. Inherent in all of the above, how- 
ever, has been the severe restriction of evaluating systems behavior as 
a function of only one variable parameter. The number of systems which 
truly fall within this category are in the minority, and it therefore 
becomes necessary to employ single variable techniques on multivariable 
problems. The result is a non-optimum compromise between system sensi- 
tivity and computational requirements, i.e., the calculations necessary 
to accurately predict system behavior in multivariable environments be- 
come unwieldly. 

Several noteworthy achievements have been made to partially overcome 
this difficulty. Extensive application of the analog-computer has been 
made to problems suitable for its use. The advent of the digital com- 
puter has opened a magnitude of possibilities, one of which has been the 
emergence of state space theory. A large potential exists in a relative- 
ly new field known as Optimal Control Theory, which is founded upon the 
maximum principal approach of Pontryagen, and Bellman's dynamic programming 
method. Stochastic processing and predicting filtering theory have become 
very influential in the field of controls due in large patt to the accom- 
plishments of Wiener and Kalman. A still different approach has dealt 
with the various possible algebraic manipulations of the characteristic 
equation, hence the name ALGEBRAIC METHODS. It is this topic to which the 
following study is devoted. Basically, these methods employ techniques 
which permit engineering systems to be designed and evaluated by simul- 
taneous variation of two parameters. Extension to one or more additional 
parameter variations is possible, but again the computational requirements 
and intelligent interpretation of results become dominant factors, just as 
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in the expansion of single parameter techniques. Little imagination is 
needed to impress in one’s mind the impact that this tool will have in 
optimizing two parameter design problems as its capabilities become better 
defined, and engineers achieve proficiency in its use. 

The following section presents the rudimentary mathematical theory 
as applied to electrical control systems employing two variable para- 
meters. Primary emphasis is placed upon formulating the necessary rela- 
tions for frequency response graphical representations. Subsequent sec- 
tions are devoted to the generation and evaluation of these displays from 
the standpoint of analysis and design. Frequency response and damping 
criteria are investigated by means of various forms of derivative feed- 
back. Noise limitations in physically realizable systems preclude the 
usage of derivative feedback in excess of second order in most cases, and 
for this reason the study of such stabilizing schemes is not considered 
practical. Typical examples are presented and a comprehensive evalua- 
tion of each is endeavored. Graphical representations are utilized and 
analyzed according to their respective value for information content. 

The concluding section recapitulates the most important limitations and 
advantages revealed by each graphical analysis example, and postulates 
further extensions of specific developments where applicable. A general- 
ized version of computer programs employed, with an explanatory flow- 
graph of each is included in the Appendix. 
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2. Mathematical Theory of Algebraic Methods. 

Conventional network analysis techniques rely heavily upon pole- 
zero representations in the complex frequency plane, as well as real time 
transient and frequency response curves. These methods are not conducive 
to the synthesis, or design problems, however. System parameters are 
seldom specified precisely, and a slight variation in one might correspond 
to an excessive alteration in another while still satisfying yet a third 
design criterion. The single parameter variation is seldom a rewarding 
approach to this type of problem, and it therefore becomes necessary to 
transform response curves into another reference frame where system para- 
meters are the variables. Fundamental complex variable transformation 
theory provides the necessary conversion media, based upon Cauchy's 
Principle of Argument. 

2.1 Transformation Techniques Employing Chebyshev Polynomials. 

The ensuing development is based upon three fundamental steps, and 
permits transforming the complex frequency plant into a coefficient 
plane and a parameter plane, where the corresponding variables are arbi- 
trarily designated as a , a and alpha (06 ), beta ( (^ ) respectively. 

X y 

The first step involves choosing a contour to be transformed. The complex 
frequency damping constant zeta ( J ) has a radial line contour in the 
s-plane, and is selected as the transformation contour. A mapping func- 
tion must be considered next, and for this the characteristic equation is 
appropriate. The third stipulation requires a transformation variable 
which allows the mapping function to be represented by two or more inde- 
pendent equations. In this case, it was decided to separate the complex 
variable into its real and imaginary parts. Listed below are the neces- 
sary steps required to achieve this transformation. Considerable attention 
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is devoted to a complete and detailed development in order to provide 
continuity in the fundamental relationships which might be readily refer- 
enced as major divergence points for certain ensuing pursuits. The basic 
relations presented in the following diagram and adjacent equations serve 
merely to establish a consistent notation for all further developments. 

The approach, although exceedingly elementary, is considered necessary to 
avoid perpetuating existent ambiguities where similar ground rules have 
not been employed. 

S = O' -h J u; 



© = 



cos 






G — - J — 



cr = uUvi 6 




C05 ^ Q = 



sin ^ © 



in © — 



e 






u) 

lOy) 



UU 






uJr}^ I “ <^05 ^© 



uJyf sin O 



Substituting the relations for cr and ix) into the top equation; 

S = - I ^ i ^ ^ ^ ^ uJr)C 0 s 6 -h J s in & 

5 = ( CCS © j 5 r © ) 

Employing the results of the mathematician Euler, this last equation may 

jS 

be written: ; — cu-p 6 

And using DeMoivre's Theorem the relation for any power of S is found 
to be; 

S ^ z=z (cos kG-hjsinkQ) (2.1-1) 
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It is now necessary to define two new variables: 



( S) — cos k © — cos (" K cos i" 5 )} 



Uk i-l) 



s,i m k & 
s I n O 



si r) f K cos (- I )) 
( cos ■' (-1)) 



Substitution of these equivalences into (2.1-1) yields: 

Since engineering problems involving stable systems are concerned with 
positive values of zeta, it is advantageous to employ the following 
relations which are obtained by trigonometric identities: 



Tk (-V = -I" I, ; U,(-S) = -I ujv 

All characteristic equations can be reduced to a general form which is 
herein defined as: 



! — / ^ "77“/ y? /-N 

h ( s) = ao 4- a, s -K s 4 - — • 5 4- s == O 

Rewriting this in summation notation, and employing the identities 
developed above, we can write: 

F(s) = ^ IT-*' +• j Z, ou^ = 0 

^ H=0 ^ k=o ^ 

^ ^ ^ K , 

P r?) = £. ujy^ cos k& -h j £ sin ke - O 

K=C K^O 

A polynomial may be set to zero when its orthogonal components are 
independently zero; hence: 

£ (-/)'■ T, (i) =0 (2. 

H=0 
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=0 



(2.1-3) 



k ,N K+-I 

(-/) 



It would now seem plausible to inspect several values of these new vari- 
ables ( l) and ),; ( § ) > in order that any discernable character- 

istic traits might be utilized. 



T„ (O = Tc(-i) = To (-1) = cc(oe) = cc,Co)= I 
TTi) = (-/)' T (-S)= -T (-f) = fcos~'r-S^) = i 



U. (V ^ (-1) Uo (-1) = 



sin (oe) _ 
sir) ( e ) 



- O 
Siv (&) 



J,(A)= (-l)^U,(-§) 



sinCc o s (~.0) I 

Sin(cos“'('-5)} 



o 



Comparison of these polynomials with a table of Chebyshev functions will 
reveal their identity. Accordingly, the following recurrence relations 
may be employed for complete tabulation. 

(5) -Z ^ Tk (E) + Tk-.r?) - o 

U K*. (5) -z ^ U K ( ^) + Uk-,(D — o 

It can be shown that T|((Q and U^C^) related by the following: 

Tk (n = 5 Uk ( p - Uk-, (V> 

Substituting this into (2.1-2) yields: 
which can be rewritten as follows: 
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Factoring a (-1) out of (2.1-3) and multiplying that result through by 
( ^ ) in a separate step yields the following sequential equations, so 
arranged for easy inspection. 



This last form of (2.1-3) is identically equal to the first term of 
(2.1-4). Therefore the two fundamental relations become: 



The utility of (2.1-5) and (2.1-6) is not obvious, but their judi- 
cious use results in a powerful tool. Employment of the bilateral trans- 
formation allows the natural angular frequency ( ) to range from 

minus infinity to plus infinity, thereby eliminating it as an unknown 
variable per se. Choosing a specific value for zeta ( ^ ) allows the 

values of U C to be obtained by the preceding method or by refer- 
ence to precomputed tables [1]. Thus any two values of can be al- 

lowed to vary, thereby formulating two independent equations in terms of 
two variable coefficients of the characteristic equation. Mitrovic [2] , 
who pioneered the above development, chose the two lowest order coef- 
ficients, Qq and CLj , for his study. Elliott, Heseltine, and Thaler, [3] 
and summarily Siljak [4] extended Mitrovic *s approach to additional co- 
efficients. A brief summary of the coefficient plane derivation appears 
on the next page. 





0 



(2.1-5) 




0 



( 2 . 1 - 6 ) 
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2.2 Coefficient Plane 

For the generalized form of the characteristic equation, let two 
variable coefficients be defined as and (X y . The order of the 

equation ( T) ) is no smaller than X , and X is greater than y . 
The value of y is non-negative. )> Y^O) Transposing the 

and (Xy terms in (2.1-5) and (2.1-6) yields the following results: 



( I) Cl^ LVy^ “f* r ( ) Ciy UJ i^y^f ( \) 

= 1 (-0%, w/ u s) 

K = C 

( ''I ) CL UJ ^ y ( i ) ~ f ) y 5 ) 

= i rp 

K=0 

K^y 

Solving the above equations simultaneously for and CL y yields: 



£ (-1) «, 
k = 0 

K 5^ y 



Uy-,(0 U^a) -Uy(i) Uk-,(V 

U/0 -iV/ 
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( 2 . 2 - 2 ) 



33 f x-y 
a (-1) W 

y )f=o 

K^X 






U, (S)My-/5>-«,., <■!)«,«) 



/•*! 



It will now be shown that: 

= M y (n Ux-I (i) — Uy., (0 Wx(^5) (2.2-3) 

Let X = Z 

u,.xs) = u,(« un) - u^(v u,., (0 

From the preceding Chebyshev recursion relations: 

U,(U = I 

U2.(5) = z i u,(d-Uo(0 ^ Z^(I) -0^ 

Therefore : 

^^.-y n) = Uy ^ ^ f^y., a) (2.2-4) 

But it has been shown that: 

U_^ = - U ^ ( U (2-2-5) 

Substituting (2.2-5) into (2.2-4), we see the general recursion formula 
as previously stated for Chebyshev polynomials: 
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Uy(U - 2 SU,_, (S) + (S) “ 0 

It can be shown by inductive reasoning that the above relation exists 
for all values of K [5]. Therefore we can also write: 

( 1 ) = U.y ~ Uy., ( 1 ) U, (V 

Thus the general solution for two variable coefficients becomes: 



a 



/ 



n 




K = 

/ 



K 

W 






Uk-Y (1) 

( i) 



a. 



y — 



n 







uu 



K 9 ^/ 






^ K-X ( S ) 

L(x>y LV 



A graphical display of the characteristic equation is often formed with 
Cl ^ as one axis and as the other axis. Such a plot is termed a 

coefficient plane. With sufficient data and interpolation, it is pos- 
sible to derive a very satisfactory design in this manner for systems whose 
characteristic equations possess the desired variables in coefficient form. 
An excellent treatise on the applied coefficient plane can be found in 
the NASA Contractor Report by Thaler, Siljak, and Dorf [6]. 



19 



2.3 Parameter Plane, Linear Case. 



The coefficient plane provides a means of working with two variable 
coefficients of a polynomial, as long as each of these variables does 
not appear in more than one coefficient. It can be recalled that each 
coefficient within a polynomial is merely the sum of the product of the 
roots taken 7T) at a time, where yr\ equals the order of the equation 
minus the order of the particular coefficient's independent variable. 
Therefore severe restrictions are encountered when attempting to use 
this method for simultaneous, but independent variation of any two roots 
of a polynomial. Thus the coefficient plane method does not include 
the more general class of problems which have the desired variables em- 
bedded within several coefficients. Siljak [7] postulated that this 
class of problems could be handled in a manner similar to that of the 
coefficient plane, provided these variables appeared linearly. Desig- 
nating the variable parameters as alpha ( ) and beta ( (^ )> one 

can rewrite any coefficient of a general polynomial as follows: 

Using this expanded form, the two independent equations (2.1-5) and 
(2.1-6) become: 




+ 




K 



u 





d ou 

h V 



K 
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k~c 



1/-0 iC JJI) + ? 7=^-')'' 



n ^ 

+ £ fl) d 

K^o ^ ^ K 



UU 



u^ (s) = o 



This rather unwieldy notation will be slightly modified by the follow- 
ing substitutions: 

Let : 



£,a, 


) = 


1 

K-o 


('-/) , 




U.'/ 


HJ 


c,rs, 


tO‘y^ ) 


1 

K = > 








(5) 


o,(^, 


<^«) ~ 


1 






W,., 


rs) 


Ms, 


k) — 


7) 

£ 








(S) 






'n 

£ 

K^o 








a) 


r\ (S 


, = 


-1 


(->r 






aj 



Employment of these substitutions reduces (2.1-5) and (2.1-6) to the 
following form: 

c(/3, ) +(3C|f?^wy^) + D,(S^u>^)==0 

■< ^ I, ( S j ~ ^ 
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These two equations are quite easily solved for alpha and beta by Cramer's 
rule, or by employing the algebraic method of elimination. The resultant 
solutions are: 



^ C I C ^ ^ ^ ) ^>1 ) ^ S > 

— Bx (% ) ‘^•n) C, (l)Uiy,) 



a = Pi ( ~ Ox( B , ( j j CJj ^ 

2 1 f S I *^> 1 ) P 2, f Bx ( ^ ) P| C 5 J ‘Ml) 

The stability variables most commonly employed by control engineers 
are the damping ratio, zeta, and the undamped natural frequency, (uOy))’ 
Most of the current literature places great importance upon these terms, 
and interprets system performance by correlation with the standardized 
second order system. In some cases, however, the rectangular complex 
frequency plane configuration consisting of the imaginary frequency 
sigma ( cr ) and the real frequency omega ( L/-) ) is more useful than 
the polar coordinate system with damping ratio, zeta ( J ) and the 
undamped natural frequency ( UJy^). In second order system parlance, 
sigma represents the rate of exponential decay of the response signal, 
while omega is the angular frequency at which this signal oscillates. 

The necessary steps for representing (2.3-1) and (2.3-2) in terms of 
sigma and omega involve the use of the binomial expansion. A brief 
outline of this method is now presented. Recalling the basic relations 
between CT , (/J , ^ and U ) in 2.1, it is possible to rewrite (2.1-1) 

as follows: 
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where the individual components of and are seen to have a 

particular recurrence pattern as shown by the following tabulation: 

S° =1 

I 

5 —(T j tu 

S^=cr^ - 3<rco^ +-j( 3 <r^uJ - 

H <r uJ -V ) 

5^— \ 0 (T^ uf^ -h S (T -h J ( S' <r^ CO - fO<r^co^-i-uJ^) 

From which: 

X 1 

Xj - cr 

X^ — <r^ ~ (ju^ 

X 3 - cr^~3<ruJ^ 

X^ = (T^ ecr^uJ^ -¥ Uu^ 

Yc — C 

Y/ ^ 

72 .=^ 2 <r 

= 3 cr ^ UJ — ^ 

>V = ^ (T ^ UJ - ^ <r uJ ^ 

T 

which can be consolidated into binomial expansion form for the general- 
ized case. 

s«=[ (y .(i) +(y J 

+ - + --- J — 0 
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ZC 

uo 



J 



S ^ [ (-/) ( ^i) cr 

_L I ^ ^ \ ^ kf -Zl't t Zt -t/ n 1 ^ 

+ J L(-‘> Uih) uj J j — 0 



where 



and 



i<) = 

(n 



k [ 



_ 



(K-X) { (XJ! 



for K ^ X 



o 



whenever ^ ^ 



Recalling the generalized form of the characteristic equation 

H-o 

we can again separate $ into its perpendicular components and 

V' 

such that ^ represents the summation of the real terms, 






and 






denotes the summation of all imaginary terms. Thus: 



f % *, = i ^ o-,U-i)(£i) 

K =0 K =0 (=o 

^ fio io (zi + l) 



In order to effectively utilize these equations for computer applica- 
tions, it is necessary that the recursion relation be explicitly formu- 
lated. This can be achieved by observing the following: 

S = <r -h j UJ 5_cr'-ju;=rrrQ 

i 

S = O' — J UJ => S - <r^+-jwJ = 0 

multiplying these two equations; 

(S - <r--j Lo ) ( S - crw-j UJ) = O 

Expanding this, we have: 

5 ^ — Z <^5 (T^ -^(X'^ = 0 

Pre-multiplying this by ^ we may write 

s'^[s^-2<r-S -i-cr^^u>^] = 0 



24 



o 



which, upon expansion, becomes 

S ^ == 

This may be rewritten in terms of orthogonal components where 



II 

4- 


^k^2. 


4- 


j 


-hZ 


K 4- ) __ 




•i- 


j 




k _ 




4“ 


j 


Vk 



-h { cv^)(X)f -h j y^) = O 

Separating this into its real and imaginary parts results in 



^ K +-2- 


- X <rX^^^ 


-f- C tr" ^-huJ ^ ) 


'X 

II 

0 


-hi 




^ (<r^ + uj^) 


Yh^O 



Noting that CT is X, , and U>* is y, , we can make these substitutions, 
resulting in the following explicit recurrence relations: 

^K-t-2.~2 .X, Xk-^/ ~*~ ( X Yi ) Y ^ ~ 0 

-2X, -^rx, 

Thus one can solve a given coefficient equation with coefficients of the 
form = for the parameters alpha and beta in terms 

of the variables sigma and omega. The method is the same as for the case 
involving the variables zeta and uJ * The two independent equations 
beccwie : 













n 








OC 


K-0 




Yk ( 


k~o 








= c:? 


(X 


Yi 


h 


Yk ^ 


/* 1^. 


n 

-h A 

kso 


A 




= o 
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which may be written in abbreviated form as 



oi O') uJ) 4- ^ C, CoTj i/j) +- 0, (<r^ vj) = O 
oc Bx(°~)<^) P C^(cr^c^) -H uJ) = C> 



where ; 



V 




e, — £ fcif 





As in the preceding development, the solutions for alpha and beta take 
the form: 



These equations are seen to be of the same form as (2.3-1) and (2.3-2). 
Thus the method of solution is identical to the one using Chebyshev 
polynomials, with only the recurrence relation and the solution variables 
being different. 

An analogous approach to this problem does not involve the Chebyshev 
polynomials method at all, as can be seen by the following procedure 
which starts by employing (2.1-1) directly: 



C I Co") b, Ccr, laj) — _C Pi (cr; to) 

C<r^ uj) Ci. (cTjUj) — Bo.(<r^ 1^) C, (ayut) 



(2.2-3) 



l3= Pi ^2- ^ — p2.C<rjOj) B,c<n<^) 

B, Ccr^Lp) Cy_Cc 7^^uj) — Bo.C^,u>) C,C<r;uj) 




K 




K=0 



k=o 



+ J £ ujy, ^ sin k& = O 

K«0 
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Separating this polynomial into its orthogonal components we may write: 
^ CL^ 

^ Kz :0 ^ ' 

Irr){FfS)j=^ ^ IVy^ ^ Siri k & :=■ O 

Now substituting the integral components of d,^ into each of the above 

K 

yields the following general terms: 

yj ^ JJ H 

oC /r ou CCS kS /S £. ccske 

K = c ^ 





+ £ 


K 

l^'r, cc ^ 


kG -= 


0 












7? 

£^ k 

/i = 0 


^ in k 6^ *+■ ,6 




^ c!n k& 






71 

+ £ 


h . 


n k ^ = 


0 






Solving these equations simultaneously for alpha and beta, we have: 



t) ^ 3 k ^ ^ > 3 K 

^ uOy^ coskOZ. stnk&—£^c^(^y^'siYikO£.c^^%co^l<9 

K~o kf-o ff -0 ^2 3-5) 

^ n t\ , 5 k 3 k . 

u'j, coikG SIT) kG — ^ Siv k 9 £.^f^f^Yx S(r)kO 

. ^ • K--.> 






£^cl^uo\oskO £ \:>yi^Yi^sirikB in kO ^ k^^uJ^caskB 






Y) 

£ 



i coskB sinks - £.bf^uJ^^sin kO £, c oUy^^cos k& 

H=^C H=0 /I'-o 



( 2 . 3 - 6 ) 
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These two alternative parametric solutions are seen to be no less compli- 
cated in form than their Chebyshev polynomial or binomial expansion counter 
parts. In fact their application is more cumbersome since there is no 
analogous recursion equation for the trigonometric functions. Therefore 
the first approach resulting in equations (2.3-1) and (2.3-2) is utilized 
in all but the most specialized cases since it was the first method to be 
successfully developed, and its application is more familiar than the sigma 
omega form. 

EXAMPLE 1 (FEEDBACK COMPENSATION FOR LINEAR THIRD ORDER SYSTEM) 

At this point it will be advantageous to select a basic system for 
initial studies. The linear case is sufficiently unencumbered by mathe- 
matical manipulations to properly illustrate certain fundamental character- 
istics, and therefore serves as an appropriate introduction to the forth- 
coming discussions of section 4. Consider an uncompensated system with a 
forward transfer function of third order, whose gain is sufficiently large 
to make it unstable. Such a system is depicted by the unity feedback 
block diagram and characteristic equation shown in figure 2-1. The corres- 
ponding root locus diagram is included in figure 2-2 and it is seen by 
Routh^s Criterion that the limit of stability is at K = 30,000. For 
larger gains than this critical value, it will be necessary to introduce 
either cascade or feedback compensation to regain system stability. Selec- 
tion of tachometer and acceleration feedback results in the compensated 
block diagram and characteristic equation shown in figure 2-3. Choosing 
a forward gain value of 100,000 results in a highly unstable uncompensated 
system. The stabilizing effect of either one of these forms of deriva- 
tive feedback without the other is shown by the cross hatched root reloca- 
tion zones of figure 2.4. The effects of simultaneous variation of both 
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forms cannot be as easily predicted on the root locus plane. At best, one 
can generate several root locus families, each as a function of only one 
variable parameter, and gain from these patterns an approximate idea of 
the required parameter values needed to place the dominant roots of the 
system in a prescribed position. The parameter plane technique becomes 
a very useful tool in this type of situation, however, since the pre- 
scribed root location automatically specifies the required value of zeta 
and also of . A cursory inspection of equations (2.3-1) through (2.3-4) 
will reveal that the following curves can be drawn in the parameter plane ^ 
which by convention has alpha on the abcissa and beta on the ordinate: 

1. Constant zeta curves, in which zeta is held fixed at a 

specified value, and is allowed to vary through a pre- 

determined range. 

2. Constant undamped frequency or "natural frequency" ( iOyj ) 

curves in which is held at a specified fixed value, 

and zeta is allowed to vary from 0 to 1. 

3. Constant ^ — cOy? curves, in which zeta and cUy^ are both 
varied such that their product is always constant at pre- 
selected curve values. 

4. Constant sigma ( c' ) curves, which represent lines paral- 
lel to the imaginary axis in the complex frequency plane. 

5. Constant "damped frequency" ( oU ) which represent lines 
parallel to the real axis in the complex frequency plane, 
and are constructed by holding cO fixed, and varying sigma 
from zero through a pre-selected negative value range. 

6. Constant real root ( cr +jO) curves, which represent 
specified fixed points on the negative real axis of the 
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complex frequency plane. (Thus it is seen that certain 
points on the s-plane are transformed into lines within 
the parameter plane. This factor plays a predominant role 
in singular line theory). 

Not all of these curves need be used by any means. The combination 
of any two of the first five curves positively identifies a complex root 
in the left half of the complex frequency plane, while the sixth type of 
curve locates the real roots on the negative sigma axis of the s-plane. 

Figure 2-5 is one type of parameter plane representation, in which Kc( “ 
alpha and = beta. It shows the constant zeta, constant U)yf and real 
root curves. One can see that for dominant roots at S = -16.5 + j33 the 
required values of alpha and beta will be at the intersection of the J « 

0.5 and - 33 curves. In this example, they intersect at alpha * 

0.0005 and beta = 0.036. The third root is real and can be determined by 
the real root curve which coincides with these values of alpha and beta. 

It is seen that the real root is approximately at S = -80. These values 
are indicated by point A on figure 2-5. 

The curves for this example are deceivingly simple, and therefore 
easy to interpret. Additional care must be employed in higher order equa- 
tions, where all the root locations are not so easily defined. In order 
to ensure the desired dominant characteristics, one must certify the loca- 
tion of all existent roots. In exceedingly complex cases, this might neces- 
sitate a second, enlarged parameter plane plot centered at the point in 
question. This, however, is a greatly simplified task compared to single 
parameter techniques. 

2.4 Paramater Plane, Non-Linear Case. 

Equations (2.2-6) and (2.2-7) have been developed to handle two variables 
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when they are individual coefficients of a polynomial. Equations (2.3-1) 
and (2.3-2), or similar forms, apply to two variables embedded linearly 
in the polynomial equation. The following development is concerned with 
non-linear arrangements of two variables within the coefficients of a 
given polynomial. Consider first the case where the variables appear in 
product form. This arrangement is termed the alpha - beta product case, 
because the coefficients of the polynomial can be represented as 
°K ^ K “‘P 

Substituting this into (2.1-5) and (2.1-6), one can write: 

cL ^ (~l) b ^ -h 8 ^ (' 0 tUy^ ^ K~ I ( ^ ^ 

Kr’O <<=0 

■h K ^ K-i ^ 

+ £ U^., (0 = O 

K=0 ^ 

(X £ h^, U ( i) + ^ f-') V 

t' K ' 

■h oi (3 £ ('O' dj, Wy, W (S) 

K —0 

+ i. (-i) ^ e U ^ { V =0 

Krg 

Employing notational substitutions similar to those in the preceding sec- 
tion, we may write: 

^'= ‘^k n) 

B,-- £ (-()'' b, w/ Wk ( V 

K = 0 

Cl— (-1) Cf^ UJ^ ^K-/ 

± ('O'" Uk ( n 
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(V 

= #=. <" '^*^- 

E;,= i e^ CJ^^ u^ (P 

thereby reducing the preceding equations to the form 

«• S, ( 5 ,u;„). 4- (3 C, H; “^>7 j 4- D,(i,i^„) E, = o 

oc 4-o(.(S D2.(j,W^„) +E2.('?,<)=0 

The functional subscripts have been included up to this point to emphasize 
the applicable dependent variables. It will now be necessary to elim- 
inate this descriptive subscript notation in order to achieve better 
comprehension of mathematical manipulations. Henceforth the dependent 
variables will be regarded as zeta and unless otherwise specified, 

and their inherent inclusion will be assumed in the remainder of this 
section. Accordingly, the above equations may be rewritten as 



c?c 


P C, 


+- oi (S 


0, 


+ 

il 


0 


(X. -1- 




■f OL p 




- Fz- 


0 



The solution of these simultaneous equations for either alpha or beta is 
possible provided their Jacobian is not zero, where the Jacobian is de- 
noted: 



J - (S, 4 f3 - fc, 4 ^ D|)( 
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Presuming this to be the case, the solution for beta is found as follows 



CB;, + PDJL + pc, + E, J = O 

C-B, - p Dj|lc< B 2 + ^0 

LC| D 2 -Cj,D| 3 j3 Lc, B?. -C;;^ B, -I- E| D, 1 (3 E, Bx~E^S,-0 

which may be written in determinant notation as 

Acp P A ^ — 0 



where 



A c D — 


c, 




^ C B 


c, 






D. 


Dz 




B, 


Bz 



and so forth. Employing the quadratic formula, beta is found to be 



n _ ( ^Ci 3 ^ ^ C^cB ^(^co ^^e') 

I - — 

(2.4.1) 

Proceeding in a similar fashion, it is possible to solve for alpha by 
the following process: 

i c^-^- ^ -»-C^/ + o^D|)P + E,~] ~ 0 

C- C, - o' f'l j L Bj. -t" ( c_2, +-0'- Dz) P Ezl ~ O 

from which: 

(B,Oz - BzD,)oC^ +(B, q-BzC, A)'=^ £,C^-£zC, =0 
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or in determinant notation: 



OC -h p) oL "h /\ 0 

The desired solution for alpha then becomes 



oC 



(^BC ^ Y ^ ^ l2C ^ 60 ^ ^c) 

2 - ^ BO 



(2.4-2) 



While the above solutions can be used to obtain alpha and beta, it is 
advantageous to simplify them for adaptation to computerized programming. 
Analyzing the terms under the radical signs, it is possible to show their 
equivalence as follows: 

^ ^ ^ £C ^ IE 0^ ^ 

Noting that ^ BC * squaring both sides yields 

which, upon expansion, becomes 

^ ^ BC^ ^BO ^ 

ec ^ ^ ^ ^ BC ~ ^ ^ CD ^ B B 

Transposing and collecting terms 

Dividing through by 4, and expanding determinants, 

(B,C^-B^C,)( JE, P,)- f B, B^P,)(BX2 .-^xC,) 

-h (C, -C2^0,)( = 0 



B, Dji 1 1 B2_ Cl E ^ B, Df ^ \ B, ^ 

-hB,C^ 1)2^2. Pi^jL Px^t ^ 
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9 

Thus, we may define the terms under both radicals as , where 

and the solutions for alpha and beta become 

0( — 

^ ~ o A 

^ i5D 

o — f-Aec A^p) ± R 

^ 

2. A ^ 

of the four possible combinations for alpha and beta, it is necessary to 
determine which of these satisfy the original simultaneous equations. 

The following results are obtained: 



D^i B,<x. ^ C, p -f- D, oc (3 4- F, J ^ 0 



— r, L Ro. ^ oC (3 4- E^2,J — ^ 

(B,P, t- (C, Dz -C^D. ) /3 l=,Px-B^D, - 0 



In determinant form this may be written: 



Substituting the above relations for alpha and beta 

(±ec _ (-^3c ^go) ± P ^A_ =0 

Z 

Therefore, the radical coefficients for one term, either alpha or beta, 
must be the negative of the radical coefficient for the other term. This 



permits the following solutions for alpha and beta: 



^ K — 



gp +■ Afcc)-+- ('-)) R 



Z A 



BP 



~ ^ Be) R 



; k = l,Z (2.4-3) 



; k=rZ (2.4-4) 
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SINGULAR POINTS 



Inspection of (2.4-3 and (2.4-4) reveals that limitations are en- 
countered whenever or ^cd equal zero. For these cases, it is 

necessary to resort to the original set of simultaneous equations for 
the complementary solution as follows: 

Should be zero, the two solutions for beta are obtained from 

(2.4-4) and then substituted into the original simultaneous equations 
where the solution for alpha is obtained by the following relationship: 









B, 



Pk Da 



Bo 



j = (2.4-5) 



When equals zero, the analogous approach yields two solutions 

for alpha from (2.4-3) which are likewise substituted into the original 
equations, yielding the following values of beta: 



P K — 



^ ^ C, Of 



(2.4-6) 

OC ^ C 2 . -^0;^ 



SINGULAR LINES 

The solution of a set of V] simultaneous equations in 7} unknowns 
is not always obtainable. When a solution does indeed exist, it some- 
times possesses unusual characteristics in that it is not a unique solu- 
tion. These cases have been categorized as singular lines, and their ap- 
plications to parameter plane theory are presented here for completeness. 
The following development is included for the purpose of showing the neces- 
sary and sufficient conditions for singular lines. 

Any set of simultaneous equations may be represented in the following 
matrix notation: 

[ f\l LXl = il3j 
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where J is the coefficient matrix of order 7ii X Yl 
cx;] is the variable matrix of order T) X I 
C B H is the constant matrix of order m X I 
and Qi B] is the augmented matrix. 

Any set of simultaneous equations is said to be consistent , (i.e., a 
solution is guaranteed to exist) when the rank of the coefficient matrix 
and the rank of the augmented matrix are equal. A unique solution may be 
found only when the rank of the matrix equals the number of unknowns, 
Otherwise the matrix is said to be singular . These three statements are 
very subtle, because their simplicity belies their importance. As an 
illustration of their application, consider the linear second order 
equations which were derived in the preceding section for coefficients of 
the form: 

Re { F( S)} ^ B, oC -f- C, P -h p, = 0 ( 2 . 4 - 7 ) 

IrnfFC?)} =, e^oC -h Cl ^ 4 - = 0 ( 2 . 4 - 8 ) 

from which: 

(X ~ ~ , — X » o _ 0| — y 

S.q-|3jC, “ H ’ ^ H 

Thus for any characteristic equation of the form 
7) ^ 

± r =0 

where the coefficients are linear combinations of alpha and beta, it is 
possible to construct a planar representation as a function of the three 
terms, X, Y and 2 > with alpha and beta as the variable parameters. 

The inconsistency arises when the denominator determinant ^ is zero, 
and when at least one of the numerator determinants is not zero. Since 
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division of a finite number by zero yields an infinite (or undefined) 
quotient, it is readily seen that the two simultaneous equations yield 
two parallel lines, which by definition are everywhere equidistant and 
therefore never intersect. Thus an inconsistent set of simultaneous 
equations has no solution. A unique solution is obtained whenever the H. 
term, or denominator coefficient is not zero. This is the most general 
type of solution, denoting the intersection of two lines at a finite 
point in the alpha - beta plane. Should one or both of the numerator 
determinants X and Y be zero, it merely places the unique value on one 
or both axes respectively. The latter case of course is just a solution 
at the origin. The third classification for simultaneous equation solu- 
tions is the singular case which exists if and only if all three deter- 
minants X, V > Z are identically zero. Such a case can exist 

when the number of unknowns in a set of equations exceeds its rank. Thus 
one or more of the equations is said to be linearly dependent. Geometri- 
cally this represents two collinear lines which by definition intersect 
everywhere, thereby giving an infinite number of solutions. Extensive 
value can be gained by the use of singular lines for design purposes, 
especially in self adaptive systems. It therefore becomes highly desir- 
able to be able to recognize the existence of such lines when depicting a 
system characteristic equation in the parameter plane. The theory of 
linear dependence stipulates that two linearly dependent equations differ 
only by a constant. Thus we may say that for (2.4-7) and (2.4-8) to be 
linearly dependent, the following conditions must hold: 

= k 6, 

C;. = K C, 

D;. k D, 
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Because of the exceptionally involved calculations using the Chebyshev 
polynomial method or the binomial expansion method, the trigometric solu- 
tions (2.3-5) and (2.3-6) are amenable to singular line synthesis develop- 
ments. Employing those equations, (2.4-7) and (2.4-8) can be broken down 
into the following equivalences: 



k e, = 


n 

k ^ 

l-O 




cci id 


= e,= 


1 

^ t t sin L & 


k c, = 


K i 

l~0 


C', 


Cor t ^ 


~ — 


V i 

^ sin lS 

i'-O ^ ^ 


K Di 


n 

K ^ 
(-0 




cos 16 




0 1 sin i & 

i -0 



Using this criterion, it is possible to define the required coefficient 
values necessary to give singular line solutions, or to determine if 
any such combinations are indeed possible for a given polynomial order. 

Various extensions can be made to include all forms of non-linear 
polynomials. Such applications merely increase the number of trigono- 
metric equivalence equations, and involve conic sections rather than 
planar surfaces, provided the independent variables are limited to second 
order. For higher orders, no geometrical interpretation is possible; how- 
ever the fundamental theory remains unchanged. 

SINGULAR SOLUTIONS 

While most equations have explicit solutions, certain differential 
equations cannot be explicitly defined by specifying the values of the 
arbitrary constants. In Laplace or Fourier transform theory, this corres- 
ponds to non-parametric equations. The solutions to equations written in 
this form are called singular solutions [8]. Certain forthcoming studies 



39 



can be facilitated by expressing solutions for alpha and beta in non- 



parametric form. Up to this point, the simultaneous equations have been 
written as: 

X, CoC, (3^ S; Wn ) = 0 

Xj_(oC,(3) ^ 

with their associated solutions for alpha and beta written as 

(X- = R ( i 

(3 = 3 C S , ) 

The following theory, therefore, is applicable to reducing the above equa- 
tions into the form: 
y (oc>p, D = 0 
with associated solutions 
oC C ( (3 ) % ) 

(3 = D C OC; 

and the alternate form 

H ( Ot; (3^ ) = O 

with the solutions 

OC = F ( (3, 0J„) 

(3 F ( <X , 

It will now be shown that for all integers j and K » 
i (-1)" n) =0 (2.4-7) 

whenever the relations (2.1-5) and (2.1-6) can be employed. That this is 
the case may be proved as follows: Given any polynomial of the form 

i fo" Uk n) =0 

K = 0 

it is possible to multiply both sides by L( j-j ( ^) *, 

^ (-1)*' "-K ^ lx (0 = 0 (2.4-8) 

h=o 
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Given also that 



£ (-1)" U^.,fi)=0 

K = o 

Multiply this by U j ( ^) to obtain: 

i i-'f Wjc'n (0 = 0 

K -C 

Subtracting (2.4-8) from this yields: 

1 (-;>*' Cwj (0 U^.i(V -Uj., (l) u^(0]=0 

K=-0 

Recalling (2.2-3) we may rewrite the above in th^ following form: 

1 (-1)" [ U^-.j ( i)j ^ 0 (2.4-7) 

KzrO 

With the aid of this relationship , it is possible to generate additional 
equations which will enable (2.1-5) and (2.1-6) to be written in non- 
parametric form. 

To derive non-parametric equations for constant zeta loci, it is neces- 
sary to eliminate uJyf from the generalized simultaneous equations which 
then become 

y i) = 0 

The necessary steps involve Sylvester's Dialytic Method of Elimination 
[9] for two polynomial equations of different order. An alternate approach, 
restricted to equations of equal order, can be used with Bezout's Method 
of Elimination but is not presented here [10]. Consider the following two 
polynomial equations of order 77) and Y) where Tn may not necessarily 
be equal to 7? , and where the coefficients of omega are functions of 

alpha, beta, and zeta. 
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Yz (^m) — 

= £. (^^) (^) i) ^y\ 

K-<3 




K 



=: O 



It will be necessary to form a square coefficient matrix out of the 
combined equations to enable the determinant to be taken. This deter- 
minant is then set equal to zero in order to obtain the non-parametric 
solution for all cases but the trival one which occurs when simply 

reduces to zero. The square matrix is of order rrj x n and is formed 
by writing C COri) Yz ^ follows: 



In a similar fashion each term in Yz(^n) multiplied by 



cation on the next page. 

Equations Yj (CU^) and y^(^n) were not constrained to be of equal 
order, so that the preceding steps might show the method to be appli- 
cable for unequal cases. Practical use of this procedure for the elimina- 
tion of to in context with the parameter plane, however, results in two 





This corresponds to multiplying each term in V I; "* 

Y? -3 

^ respectively. 




respectively. The matrix equivalent is shown for amplifi- 
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equations of equal order when the following substitutions are made for 
and q y : 

f, = (-0%, u, (V 

9k=H)'^o,^ 



The non-parametric equation for constant is derived by tactics 

analogous to those used to generate the constant zeta equation. Again, 
a set of equations is obtained for the purpose of taking the determinant 
of the coefficient matrix which, by definition, must be a square matrix 
if the determinant is to exist. Since the variable matrix is 

only of the first power, it is possible to generate 7)"/ additional 
equations by use of (2.4-7) which was previously derived. Letting j 
range from zero to 7) I in (2.4-7), the following set of simultaneous 
equations is obtained: 
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a;/ 
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u,., CS)^- 0 


i- 

K^o 


Hf 






tr-2 (V ^ ^ 


K=o 
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K^C 
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i-\f 


\ 
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Ti 


U (0 ^ 
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H~o 


fi)' 






^ >f -( yi-x; (0~ ^ 


k 








IX ^.(y,.,) (J)= O 



Writing these equations in expanded form; 

(- l)^ ^0 ^0 (0 f-l) u, ( §) ^ ^ ^ -f- f-i) 

('lya, LUy,^U -,(0 ^ (-l)o.,iifn Uo(V ^ ■*" 

(~l) ^Yi (~0 ^o({) +* *' ^rt-X^ 9~0 



(-l)^ Ao U -(rf~i) (~0 (^n ^ --(n-Z)^ H- • • • 

. c . + (-i)''a^ U, (§) = 0 
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These equations may be simplified by substituting the following 



relations: 

Uo (0 = c 

Thus we may write: 

O (-1) a, cu„ •• • +H)"a„ w"Ur,(S) = 0 

~(-l)° U,(l) -t- O ■ ■ ‘ +(-lf 



I 

I 

— (-lY CL o Ufi~, (i) — f-/) • • • -f- C”/) U I 0 



Which is seen to be a set of n-1 equations with n-1 unknowns. In 
matrix notation these become: 



(-/) a, uj„' + (-1) -t- ■■■ 

|'^0^O;L^n ~(~0 + (- 1 ')'^ + ••• 

’Mf,' + 



U, 

U, 

U3 

Uh 



= 0 

(2.4-12) 



rt -2 n-a 



j(~l) 

L 



('-oV.w; 



u„-, 

L J 
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The above developments have shown how a set of two order linear- 

ly independent equations with four variables can be reduced to one equa- 
tion of order n-1 with only three variables. This is accomplished by 
solving the coefficient matrices of (2.4-10) and (2.4-12) which elimin- 
ate zeta and omega respectively. The application of these methods will be 
indicated in the following article where the mathematical concepts of the 
frequency response curves are presented. 



2.5 Frequency Response Applications of Parameter Plane Theory 

Frequency response curves can be portrayed in many ways. Some of the 
more useful ones will be considered in sectipns 3 and 4. The various 
representations of frequency response curves can be derived from the 
following fundamental theory relating magnitude to frequency. Consider 
again the transfer function of a mathematical model which represents the 
system to be studied. This can be reduced to a ratio of two polynomials 
by partitioning, or other suitable means to arrive at: 



T(s) = 



gp ■f-g.s 

+b, S +b2.S^ ^ ^ h' 

j-0 



ivrs) 

DCS) 



(2.5-1) 



To find the absolute magnitude of this transfer function, the following 



is noted: 

M ^ s 



T 



4 

= [TCJw)^ [TC-Jw^ 



(Although the notation standards of sectipn 2 clearly state the distinc 
tion between LU and lO , it is important to note that natural 
angular frequency is no longer being used.) The following steps will 
illustrate this transition: 
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T (. j uj) = 



+ a,jto' -h -I- a^(jujy 



a-n 0 '“-') 



t>0 b, jw- *- 

Rewriting this to separate the real and imaginary parts, we have: 

a, f- ii uj)^ -h • • ^ a, jco a2 • " 



T (ju.) 






T (ji> 0 ) = 



r - a 2 ^ LV^ -ha,^ - •• J -i- j 

[^ho - bx ^ ‘ * 'J ^ -h^ 4> 



In abbreviated notation, this becomes: 

T(l'jj) = lVeyg» -r- N oJJ I juj) 

C) Q ^ f) C J 0 oct .■! 

Preceding in a similar manner we find the following: 



T (-juO 



_ [a. -«i 



+ o^u;‘< - ■ ■ - J -j [a.w 



E b„ ~ bj. • ' ’ J ~J H “ bj bi' 



■3 



Employing abbreviated notation again, we may write: 

^ fVevev) ( 

T (-Ji^’) = ^ 

eve n (ji^) — C vaJ ( J 

Thus the square of the absolute magnitude is found to be: 

e v'e n ( J ^') N oJ J (J J L ^ S v s n ( J ~ hJ 



= 






M = 



[D 9 ve^- 0 ^>'’ b’oJJ Oeve n (jbo) - D odJ (j^)] 



c J Nocid ( J uj ) ^ 

Osvf^n ( J J — P 0 d cl 



(2.5-2) 







[ b. - b^_ + h^uo‘' - hg u; '- 4 • • • 





47 



This can be written in the following suiranation notation: 



%rx 






\jj 



M 



K = 0 
: ev^eo 






^Yl-I 

£ (-0 
K o(^d 






Xyy) 



^ (-0 ^ d, 



Ji-0 
a ev€r) 



aJ 



+- 



Xm - 1 Jt - / js 

^ (-/) ^ 



x-\ 

X oJji 



V , k 

t (- 1 ^ ^ 

« ^ K=0 



zn 



ZK 



UU 



m 

^=0 



(- 0 ^ d,j, 



Let 




- ^%X 

Substituting these new variables above, a more concise notation may 



be formalized for easier illustration of pertinent relationships. 



'n , K K 

^ ^ (-0 ^ 

(V| _ K = 0 



#0 f-')' 



(2.5-3) 



In general, the coefficients of (2.5-1) are non-linear combina- 
tions of at least two variable parameters, designated alpha and beta. 

The linear combinations were discussed in article 2.3 and the associated 
form of any polynomial coefficient was found to be 
h ^ OL -h ^ 

For the non-linear case, the simplest combination was seen to be the 
alpha-beta product combination. This subject was treated in article 2.4 
where the resultant coefficients of a given polynomial were broken down in 
the following manner: 
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ct + (3 + 0(.(3 

Theoretically, there is no limit to the degree of difficulty which exists 
in the non-linear case. A few examples of the complexity involved in 
non-linear parameter combinations are presented here for completeness. 

oL -h (3 cl oC (3 -H 



CX ^ t o( + 


Ck r- + 


J 


K -f- >• V 


-f 








b ^ oc +• 






■*- +■ -f^ 4- ^^(3^ 


+- 




-h 


'■ K ^ j K P ^ 






4- 




-h 




> 





Clearly, it can be seen that such combinations are exceedingly unwieldy. 
Attempts to handle all but the simplest arrangements are extremely labor! 
ous. Computer programming is still in the embrionic stage for parameter 
plane applications. Attempts to utilize this tool have been extended 
only to the alpha-beta non-linear case. This restriction is acute when 
considering the frequency response curves. A cursory review of the 
mathematical developments at the beginning of this article will show that 
the polynomial coefficients of (2.5-3) are obtained by squaring the poly- 
nomial coefficients of (2.5-1). Thus, even the simple linear case of 
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article 2.3 results in non-linear combinations of second order for fre- 
quency response studies. To handle this problem, it becomes necessary 
to depart from the procedures which were employed in previous develop- 
ments. By assigning a value to either alpha or beta, the task is reduced 
to merely a quadratic formula solution for the remaining unknown variable. 
This can be illustrated in the following manner: Given a non-linear co- 

efficient combination of the alpha - beta product form: 
aK= oc 4- p + d ^ o( p + 

Assign a value to alpha, and denote the value as A: 

This may be rearranged for clarity, 

^ E (2.5-4) 

where the first term on the right hand side denotes a constant value, and 
the second bracketed term represents a coeficient of the independent vari- 
able beta. In a similar manner, a fixed value may be assigned to beta. 

For illustration purposes, let it be denoted as B. This gives: 

a = h ^ oL -h B 4- d ^ o( B 

which is likewise rewritten as: 

a [ c^ e + ej + [t>if + oc ( 2 . 5 - 5 ) 

Again, the first bracketed term represents a constant while the second 
denotes the coefficient of the independent variable alpha. The computer 
program for frequency response curves, program PARAMS, is based upon this 
solution procedure. Consider first the case where magnitude and omega are 
specified. By assigning incremental values to either alpha or beta, the 
corresponding solution may be obtained for the other. This is the method 
employed in the graphical representations which have both alpha and beta 
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on the axes. For the graphs where magnitude is an axis variable, it is 
incremented while a corresponding value of beta is obtained as a function 
of alpha and omega. The same philosophy is utilized when omega is an 
axis variable. For the case where both magnitude and omega are on the 
axes, the values of alpha and beta are pre-selecte4 and the familiar Bode 
plot results. All of these graphical arrangements are treated with con- 
siderable detail in section 3. Their introduction at this point was 
simply to point out that the solutions are based on (2.5-3) and either 
(2.5-4) or (2.5-5). Since (2.5-3) is the fundamental gain versus fre- 
quency equation, it can be anticipated that all frequency response para- 
meter plots will contain the same basic information as a series of con- 
ventional Bode frequency response curves. When two variable parameters 
alpha and beta are utilized together with magnitude and omega, the problem 
becomes one of four dimensional space. When still another polynomial co- 
efficient variable parameter is considered, an additional dimension is 
added. In spite of this, it is still possible to derive one set of para- 
meter curves from another, since each is just a different projection from 
the same multi-dimensional model. The value of each separate representa- 
tion lies in the ability to depict certain important characteristics which 
would not be obvious otherwise. Additionally, certain configurations 
facilitate the study of a dominant parameter, or the effect of certain 
parameters on different performance requirements. Certain mathematical 
considerations must be pointed out at this time concerning the inter- 
dependence of these graphical plots. Clearly, the construction of PARAM-7, 
the Bode plot, offers no problem. In this case, the desired values of 
alpha, beta, and omega are inserted into (2.5-3) and the corresponding 
magnitude squared value is then computed. For the remaining programs 
(PARAM-1) through (PARAM-6) , a rearrangement of (2.5-3) is necessary since 
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magnitude is no longer the unknown variable. As an illustration, con- 
sider (2,5-4) where alpha has been specified. Beta then becomes the un- 
known value to be solved. For added clarity it is better to revert back 
to (2.5-2) in order to point out why a second order system is always en- 
countered in the linear. case, as well as the non-linear alpha - beta pro- 
duct case. Since beta is the unknown variable, pre-selected values of 
magnitude, omega, and alpha must be presumed. Performing these substitu- 
tions, (2.5-2) may be written as 

H ^ 

where, using the notation of (2.5-4) 

f K [ t. fl + ej - E e^J uj^ + ■ ■ • 

^ ^ E '~h 1^0 ~ b- <^2. ^ J ^ b- [ b~ ci fjj Uf ^ — 

[ t>d fl + e.] - I hzfi b ■ ■ > 

1 ' ■ ■ 

I 

I 

I 

= [c, b- difl] - [ £j. + fl] i- b du, fl'J LO^ 



Therefore, one can easily see that by squaring the right hand bracketed 
values and combining like terms, 



M 



Sb/6o 4 -swi3<p + Sib ez p 



+ SPSI f 3 + SoSZ ( 3 * 
where the terms SNBO, SNBI, SNB2, SDBO, SDBl, and SDB2 are so chosen be- 
cause they are the variables used in the computer program to be dis- 
cussed later. Since the numerical value of has been specified, this 

equation can be cross -multiplied. Transposing the right hand side of the 
resultant equation, one can write: 
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SNBC -M^-SDBO +\_SNBl -M^SPiSlJfi +LSA/fi2 -M^SD<3i]/3^=0 

This is clearly a quadratic in beta, and employment of the quadratic equa- 
tion yields the desired results. Imaginary numerical values have no 
physical significance for graphical plotting and are therefore discarded. 
However both the plus and minus signs in front of the square root term 
of the quadratic formula must be considered. Thus it is possible to 
have two values of beta for a specified combination of magnitude, omega, 
and alpha. This is an important consideration in graphical plotting 
techniques which can lead to spurious lines if not properly accounted 
for. Effectively the problem requires calculating all the prescribed 
graph points using first one sign of the square root term, and plotting 
them. Then the computational process is repeated with the alternate 
sign of the square root term. This illustration has shown how a solu- 
tion would be obtained for beta when alpha, magnitude, and omega are 
known. Precisely the same procedure would be followed to obtain a solu- 
tion for alpha when the values of beta, magnitude, and omega are given; 
only (2.5-5) would be used in lieu of (2.5-4) 

RESONANT PEAK MAGNITUDE 

Performance specifications of certain systems include stringent re- 
quirements on the presence of a stipulated resonant peak magnitude. Con- 
versely, other design studies are concerned with the guaranteed absence 
of a resonant peak magnitude above a pre -determined level. Either of 
these requirements results in the study of frequency response curves with 
particular emphasis on the peak magnitude and the frequency at which it 
occurs. To facilitate such studies, it would be advantageous to construct 
resonant peak magnitude curves for the system ss a function of the 
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polynomial variable coefficients, alpha and beta. This involves the 
utilization of the non-parametric theory presented at the end of article 

2.4. 



Normally, one thinks of only one maximum point in the entire fre- 
quency spectrum because of dominant root characteristics of low order 
polynomials. Certain high order systems, however, can be represented by 
characteristic equations which have more than one pair of dominant com- 
plex roots. Such arrangements will yield more than one maximum point on 
the frequency response curve. Only one of these points can be consider- 
ed as the global maximum, while all other peaks must be relegated to the 
role of local maxima. To avoid ambiguities, the notation M|o(u^)will be 
restricted to the global maximum point only. The corresponding value 
of omega will be denoted by p . Since it is desirable to refer to 
local maximum points as well, the term M shall be construed to be 

applicable to all resonance points, both local and global. The respec- 
tive angular frequencies at which these occur will be labelled as UJ^ . 
Thus it is possible to have any number of M^(lo) and U) y. values for a 
given transfer function. Only one of these can be the /V|| 5 ^iand OJp , 
however. 

Recalling (2.5-3) where 6^ and are both functions of the vari- 

able parameters alpha and beta, it is possible to state this equation in 
the following functional notation: 



M ^ M 



f, P; 

fj. ^ 



(2.5-6) 
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Assuming at least one pair of dominant complex conjugate roots, it is 
possible to find the following relation from (2.5-6): 






3L 

r 



f I ( Ci ) (^ > 



^ uj) 



UJ-’jJyr 



Co<., (5, Wy) 



Ilj = UJp 



fc (®; 



(2.5-7) 



(2.5-8) 



Fundamental theory of equations dictates the following correlation be- 
tween (2.5-6), (2.5-7) and (2.5-8): 

1. A relation in the form of (2.5-7) exists whenever the 
slope of (2.5-6) is equal to zero and its radius of curvature is nega- 
tive. Expressing this statement in mathematical terms: 



(j [m ^ (u^)] 
M ^ ((/')] 



= 0 
= 0 



d ix> 

2. Only that value of 
mizes (2.5-6) satisfies (2.5-8). 



(2.5-9) 



(2.5-10) 

UJp in (2.5-7) which globally maxi- 
Thus it is possible to state: 



M 




7T/ ^ X i rm u/y\ 




UJ = 






(2.5-11) 



The preceding theory is adequate in that it gives the necessary and suf- 
ficient conditions for the existence of Satisfying the above 

criteria is fairly straight forward whenever the values of alpha and beta 
have been selected or prescribed. The dilemma arises when seeking the 
algebraic relationship for Mp (l^) as a function of alpha or beta. 
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Although three equations are available, the inequality in the second de- 
rivative criterion of (2.5-10) imposes a severe limitation on the solu- 
tion technique. The most straight forward approach would be to utilize 
(2.5-7). From these it is possible to obtain (2.5-11). Such a method 
of attack is only useful when both alpha and beta are given. This would 
require a two dimensional scan for all possible values of alpha and beta. 
Various computerized programs exist in optimal control theory for the 
determination of (2.5-7). These, however, are based on apriori knowledge 
or random search techniques using Lagrange multipliers. Not only are 
these methods time consuming, they are not conducive to highly sensitive 
frequency response curves where resonant peaks are in close proximity to 
one another. An alternative solution employs the non-parametric theory 
which results in (2.4-10). This procedure is applicable only when the 
two equations involved can be set equal to zero. Therefore it is neces- 
sary to investigate the locus of inflection points of the magnitude versus 
frequency curve. From this information, it is possible to determine how 
many local maximum points exist. It can be shoxm that for T1 inflection 
points, there can be no more than Tl maxima. In practice, however, only 
V - I will be found for most systems. This procedure has the important 
restriction of eliminating omega. Thus it is impossible to employ (2.5-6) 
in order to find the value of M at each inflection point. At best, 

this result makes it possible to construct the locus of inflection points 
in the parameter plane when alpha and beta are the axis variables. It 
then becomes necessary to employ some other means of determining which side 
of this locus represents the minima, or maxima. Full utilization of the 
non-parametric theory can produce this desired result. Inspection of the 
developments leading to the coefficient matrix of (2.4-12) will show that 
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the combined application of (2.4-10) and (2.4-12) provides the necessary 
information for (2.5-6) as follows: 



1. The results of (2.4-10) yield a solution of the form: 



X 



f 7 P ^ 



(2.5-12) 



2. The results of (2.4-12) provide the following alternate 



solutions : 



y, ~ i 



(2.5-13) 



or 



>2 - f, ( 



(2.5-14) 



depending on which coefficient variable is eliminated. From these re- 
lations, it is possible to employ (2.5-6) directly by use of (2.5-12) 
and (2.5-13), or (2.5-12) and (2.5-14). In this manner, it would be 
possible to employ the optimal control theory previously mentioned. 
However the more direct approach would be the utilization of the gradi- 
ent technique which is presently employed in many optimal control pro- 
blems. Alternately, either of these pairs may be solved simultaneously 
to provide a locus of inflection points as a function of the two re- 
maining variables, alpha and omega, or beta and omega. This latter ap- 
proach would be most desirable for frequency response parameter plots 
employing these variables as the graphical axes. 

The preceding theory has shown several ways of handling the peak 
resonant magnitude problem. Inherent in each of these, however, is a 
vast array of mathematical manipulations which is limited by solution 
techniques presently available. The following development will indicate 
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this more clearly. Consider again the fundamental relationship of (2.5-3) 
upon which these derivations are based. The following result is obtained 
for (2.5-9): 



^ tA (X) 
i * 



^ f-l) 

U = o 



^ r-') 

K=0 



rl 7 



(2.5-15) 



Let g^= numerator coefficient of (2.5-15) with the sign embedded, and 
let f) « = denominator coefficient. Thus (2.5-15) can be rewritten as: 



^ A! (X.) 
c) y 



r>5 -f-yi-i 

X 

C -0 



/I 

£ (-ir 



The equivalent form of (2.5-10) then becomes: 






L I, 



(2.5-16) 



I'l -T If - 1 _ 777 / /p /_/ -7 



2/-/ ■ 






Thus from (2.5-16) it is noted that the order of the denominator poly- 
nomial is raised to a power of four times the order of the denominator in 
the original transfer function. In like manner, it is evident that the 
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numerator polynomial of (2.5-16) is also of very high order. Therefore, 
all but the simplest of transfer functions would require computer assist- 
ance to solve the determinants of (2.5-15) and (2.5-16). Although many 
programs exist for explicit determinant solutions, these cannot be used 
when one of the terms is an unknown variable. Until a generalized com- 
puter program is developed to handle these cases, the construction of peak 
resonant curves can only ba achieved by a laborious calculation of second- 
ary curves. Such procedures do exist, whereby tangent lines are drawn 
from a series of curves. To be highly accurate, however, an excessive 
number of these are required for sensitive systems. 

PHASE CALCULATIONS 

The preceding discussions in this article have been limited to magni- 
tude studies in the frequency domain. An equally large field of interest 
is that of phase angle characteristics. The utility of the magnitude 
versus frequency plot lies not only in analyzing or predicting the fre- 
quency response of a system, but also in the ability to check its stabil- 
ity as well. The correlation between magnitude and phase curves provides 
this verification by means of phase margin and gain margin calculations. 

It therefore becomes desirable to derive the phase calculations as a func- 
tion of the three variables: alpha, beta and omega. 

Initial developments in this article showed how (2.5-1) could be alter- 
ed to produce 



^ ev e n ( j ^ odd ( J 

T(jtL) ^ 

ev(5ri + Dodd ( -1^) 
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Considering the numerator and denominator separately, either can be re- 
garded as the real and imaginary components of a vector with a given magni- 
tude and angle. Magnitude calculations have already been developed and 
are of no interest here. Phase angle calculations, however, are known 
to involve the arc-tangent of the vector components as follows: 



Phase P = 



Im -j 


> T C jt^)} 




Re 


[ T(j w)' 


1 J 



(2.5-17) 



Two solution procedures are therefore possible. One is to solve independ- 
ently for the angle of the numerator vector and subtract the denominator 
vector angle from it. Thus: 



Nad 0“^) 


-1 

— tcL r\ 


Pcold 


eVeh ( j 


p€\/e^ 





(2.5-18) 



The alternate approach involves only one trigonometric conversion. Clear- 
ing the imaginary term from the denominator of T~ Cjtt^) , 
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Therefore, substituting this into (2.5-17), 



e 



-f 

tan 



fs/odJ (juj) De^/er}(J(^) Nev'en ( jt^) OoJiJ 

N even (j>^) pQveYi(}u))’-Nodd[\{Aj) Oodd(j(^ 



(2.5-19) 



Either (2.5-18) or (2.5-19) may be used for the prescribed phase calcu- 
lation. For the Bode plot (PARAM-7) , both have the same degree of diffi- 
culty in computer applications, although (2.5-18) is preferred on long 
hand calculations. Extensions to phase curves in the remaining para- 
meter plane frequency response graphs (PRAAM-1 through PARAM-6) have not 
been attempted. These would be possible without the use of the quadratic 
formula if (2.5-18) were used. Selection of (2.5-19) requires the usage 
of the quadratic formula with both its square root coefficients as il- 
lustrated in the magnitude portion of this article. 
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LWCOMPENSATL'D THIRD JRbDR SYSTEM 



CHARflCTDRlSTI C EQUATION 
•+ GO S~ SOCS ^ K—O 



F I & u c 2 - / 




/-iCHLE^I.OOE+Ol UNJT5 INCH. 
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’'LhUIS. RUOT LOCUS FOR 5^'i.3+60S^-i2+500S+K-0 
CftIN,(l<J URRIE.S FROn 0 TO 100,000 

FIGURE 2-2 
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THIRD ORDER system WITH TACHO METER 

P, H D P C C B / £ I' m V V EE £ P B /? f K 

CHfinPCTLRl 5 ru EOUHTION 
S’ + ("6/:' ■/ (soo + KKi)S + K = O- 

■ r ' I e II p £ 1-3 



64 




. rHLE=i.OOE+Oi UNITS INCH. 

V-3CRLE=i.00E+0i UNITS INCH. 

GLfiUlS, ROOT LOCUS FOR S**3+60S^*2+5005+K=0 
FEEDBRCK COMPENSRTION ROOT RELOCRTION ZONES 

FIGURE 2-4 
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X-SCflLE=1.00E-03 UNITS INCH. 

'■ SCflLE=1.00E-02 UNITS INCH. 

i-LRUIS, PRRRMETER PLRNE CURUES 

(60+100, DOOR) S**2+ (500+100, OOOB) 5+100,000=0 

FIGURE 2-5 
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3. Computerized Graphical Construction Techniques For Frequency 

Response Studies. 

Graphical analysis has always held a virtuous distinction in all 
facets of industry and business. Since the Second World War, however, 
the advances of technology promoted the engineering profession beyond 
the existent graphical methods which were derived for the single purpose 
of varying only one parameter at a time. In spite of this, the ancient 
adage "one picture is worth a thousand words" has withstood the test of 
time, as witnessed by many practicing engineers today. Many pictorial 
representations of single variable systems are familiar; however with the 
advent of multiple parameter variations, graphical representation in two 
dimension's becomes more complex. Considerable experience and moderate 
imagination are necessary to visualize three dimensional configurations 
when presented in two dimensional planes, in order to adequately make an 
intelligent intrepretation for meaningful results. 

The recent developments employing algebraic methods for variation 
of two parameters have resulted in the successful use of digital com- 
puter facilities to perform the required manipulations fpr graphical 
representation of mathematical equations. The preceding section was 
devoted to a fundamental development of the underlying mathematical theory 
employed in two parameter design studies. This section will include a 
broad presentation of certain graphical methods which are most likely to 
be employed in two parameter electrical control problems, with primary 
application to frequency response studies. This section is also con- 
cerned with a detailed inspection and description of the computer pro- 
grams utilized in formulating those graphs. The discussion is not direct- 
ed toward programming, but rather an effort is made to clarify certain 
computer oriented problems so that engineers who desire to utilize any 
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program contained herein may do so with minimum difficulty and yet gain 
maximum benefit from its comprehensive construction. Illustrative ex- 
amples of each of these programs appear in sections 4 and 5 where several 
case studies are investigated in depth. 

3.1 Selection of Graphical Representations. 

For any given set of six variables, the study of their respective 
inter-relationship by means of two dimensional projections results in a 
myriad of graphical portrayals. The magnitude of the problem can best 
be seen in the following development: Assume that an arbitrary number 

of variable parameters are to be considered. If T1 is the number of 
design variables, it is possible to show that the number of possible 
graphical representations is 71 -factorial ( 7? f ). All graphs reveal 
essentially the same information content when one axis variable is switch 
ed from the ordinate to the abcissa, and vica versa. Therefore, redun- 
dancy is eliminated by allowing the abcissa and ordinate variables to be 
interchangeable. This reduces the number of combinations by a factor of 
two, thereby making the revised number of graphical representations equal 
to % ( 7? ! ). Although a considerable reduction in the size of the ar- 

ray has been effected by interchanging the axis variables, it is highly 
desirable to reduce this size further, since a six variable array means 
360 graphs, a five variable array requires 60 representations, and so 
forth. To further reduce this compilation, it is often possible to inter- 
change any two parameters alpha ( ) and beta ( P ). Such would be 

the case in the parameter plane technique where they are, by definition, 
arbitrary variables of the characteristic equation. Even this reduction 
leaves a formidable task of evaluating over 180 graphs for a system in- 
corporating six variables, and any system with five variables requires 
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an array of 33 different graphical representations. For this reason 
alone, the present state of the art is limited to at most four variables 
for practical design techniques employing graphical displays. Any system 
with four variables requires basically seven different graphical presenta- 
tions for comprehensive evaluation. Seldom are all of these required for 
a satisfactory design, although barring the elements of time and expense 
they would most likely be necessary for optimal results. Effectively the 
problem at hand involves the inspection of four independent variables. 

This can be regarded as a four dimensional space which is not physically 
realizable. By assigning a constant value to one of these four variables, 
the problem is reduced to a three dimensional sub -space representation 
which does have a physical shape. It*s meaning, however is restricted to 
that particular value of the fourth variable. To present this three di- 
mensional sub-space on a two dimensional plane requires that a large 
number of curves be projected onto this plane, where each curve represents 
a different value of the third variable. The following discussion is 
based on this concept, and the resultant graphs should be interpreted ac- 
cordingly. 

Most specifications for the dynamic performance of linear systems are 
written in terms of the time domain and usually include such terms as peak 
overshoot, transient duration, and steady state error. It therefore be- 
comes necessary for the control engineer to bridge the time domain with 
the more utilitarian frequency domain. Correlation is accomplished by 
means of the Fourier integral when rigorous methods are called for. This 
treatment is avoided in practice however, because of the inherent arduous 
task involved in all but the simplest of integrations. It can be shown 
that by suitable approximation, many linear problems fall into the cate- 
gory of second order systems. For engineers not experienced in partitioning 
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methods and root locus techniques, this statement would seem to be easily 
contested; however, such an analogy actually exists in most cases. For 
systems which cannot be resolved into an equivalent second order approxi- 
mation, the absence of dominant complex roots is evidenced by the shape of 
the closed loop frequency response curve. Since considerable information 
exists in present day literature on the generalized second order system, 
much time and labor is saved by manipulating a given problem into this 
form whenever possible, and extracting the equivalent damping and band- 
width specifications from pre-computed curves and tables. 

A cursory review of the terms used in control engineering applica- 
tions will reveal that some of the more common design criteria are steady 
state magnitude ( M ), phase angle designated theta (0 ), natural 

angular frequency designated omega ( bUy, ), and damping ratio designa- 
ted zeta ( ^ )• Some problems are more amenable to the rectangular 

coordinates of the complex frequency plane. These are designated as the 
damping constant sigma ( cr ) and the damped angular frequency omega 
( bU ) . It is again emphasized that two forms of angular frequency are 
designated by the word omega. Much of the existing literature uses this 
term very loosely, and extreme care must be exercised not to mistake un- 
damped, or natural angular frequency ( ^y) ) for damped angular frequency 

( UJ ). The relation between these two variables is found in article 2.1, 
and it should be noted that these terms are synonomous only when the damp- 
ing ratio ( ^ ) equals zero. 

Section 2 developed the notational concept of alpha ( ) and beta 

( P ) as representing two arbitrary variables embedded within the co- 
efficients of the characteristic equation. The preceding paragraph has 
shown that eight variables are frequently encountered in the resultant 
electrical engineering control problems. As previously stressed, the 



70 



simultaneous usage of all of these is not feasible. Fortunately, however, 
one may break the control problem down into two basic categories: stabil- 

ity and frequency response. These two facets are interrelated, however, 
and for that reason some attention must be given to system stability 
when performing frequency response studies. 

The parameter plane stability problem may effectively be considered 
as the study of the relation between the characteristic equation poly- 
nomial and the stability performance variables ^ and , or O' and 

LO . When the two parameters alpha and beta are embedded within the 
characteristic equation, the scope of the problem increases to a study of 
at least four variables. Thus one might be inclined to consider seven 
possible graphical arrangements. Unlike the frequency response studies, 
however, the stability problem has been attacked from a slightly differ- 
ent viewpoint. Only alpha and beta have been plotted along the axes, and 
no attempt has been made to develop the necessary equations and programs 
to portray a different set of graph coordinates. Thus only one basic 
graphical presentation is used for parameter plane stability studies. The 
example in section 2 illustrated this method of stability study by employ- 
ing the following alpha - beta parameter plane curves: constant zeta ( f ) 

curves; constant omega ( ) curves; constant zeta - omega ( 

curves; constant sigma { CT' ) curves; constant omega { lP ) curves; and 
real root ( CT + J C ) curves. Although the constant zeta - omega ( 
curves are identical with the constant sigma ( <7^ ) curves, the alge 

braic methods required for each are different. The names assigned to each 
have therefore been kept separate to avoid confusion in the respective 
graphical construction techniques. Additionally, the real root ( CJ' + j ) 
curves represent constant sigma points in the complex frequency plane which 
are transformed into lines on the parameter plane. As such, they are not 
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to be confused with the constant sigma ( ) or constant zeta - omega 

( - u) ) curves. Although partial programs exist for these curves, a 

program entitled PARAM-D has been written for the primary purpose of 
consolidating these programs into one package for added flexability and 
versatility. This permits any combination of the above curves to be pro- 
cessed on one graph. 

Complex engineering problems often involve a stability study of systems 
with more than two arbitrary characteristic equation coefficients. Pro- 
gram PARAM-D has been extended to include a third variable which is here- 
in defined as gamma ( Y )• This scheme will allow stability analysis 
studies to be visualized in three dimensions, rather than two. When con- 
fronted with this type of problem, it then becomes essentially a four 
variable problem. Therefore it is necessary to hold one of these four 
independent variables constant, and generate a series of curves using 
the remaining three. This is accomplished by assigning a constant value 
to one of the stability performance variables: ^ ^ uJy^ ^ cT+jO. 

The remaining three variables are those in the characteristic equation co- 
efficients, namely alpha, beta and gamma. With these, the program at- 
temps to develop a three dimensional surface for the fixed value of the 
specified stability performance variable. The x, y and z axis variables 
are alpha, beta and gamma respectively. Graphical presentations are re- 
stricted to two dimensions, however. The dilemma therefore is that this 
three dimensional surface must be projected onto a two dimensional plane, 
which is herein defined by the x and y axes. The only way that this pro- 
jection can be accomplished is to construct a large number of third vari- 
able (gamma) curves on the graph. Each curve then represents a cross 
section of the three dimensional surface. The physical shape of the sur- 
face can then be visualized providing a sufficiently large number of 



72 



these curves are thus constructed. 



The frequency response study has been viewed in a manner which is 
slightly different than the above stability program approach. It was 
initially decided to make a thorough analysis of the four variables: oL , 

^ » M » and U) . As such, alpha and beta can be regarded as two 

variables embedded in the characteristic equation coefficients, while 
magnitude and omega are the two frequency response performance criteria 
of interest. Treating these as four independent variables results in 
seven basic graphical combinations. These have been titled PARAM-1 
through PARAM-7. Additionally, phase is investigated in PARAM-7, where 
its correlation with magnitude is most easily made. All possible graphi- 
cal presentations of interest have been utilized, and a summary of these 
is presented in figure 3-1. They have been implemented in one consoli- 
dated computer program PARAMS which is included in appendix I. A brief 
summary of these graph characteristics is presented below, while a more 
thorough discussion is reserved for article 3.2. 

PARAM-1 is a graph of alpha versus beta, with magnitude held at a 
fixed value while drawing constant omega curves. Interchanging alpha and 
beta produces absolutely no alteration of curves, other than a ninety de- 
gree re-orientation. PARAM-2 is a graph of alpha versus beta again, with 
omega held at a fixed value this time, while constant magnitude curves are 
drawn. As in PARAM-1, no difference other than a ninety degree re-orienta- 
tion is obtained when the axis variables are interchanged. PARAM-3 places 
alpha on the ordinate, and magnitude on the abcissa. Omega is held at a 
fixed value throughout the computational process, while constant beta 
curves are generated. Now interchanging alpha and beta in PARAM-3 will 
result in an entirely different set of curves, since the sensitivity of 
these two variables is normally different. PARAM-4 is a graph of alpha 
on the ordinate versus omega on the abcissa. Beta is held constant while 
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magnitude curves are generated. Again, different curves will be obtained 
when alpha and beta are interchanged, since holding one characteristic 
equation parameter constant will produce an entirely different effect if 
another parameter is held at a fixed value instead. PARAM-5 generates 
constant omega curves with beta held at a fixed value, and graphs these 
curves against magnitude on the abcissa, and alpha on the ordinate. As 
with PARAM-4, the PARAM-5 curves are quite different when the roles of 
alpha and beta are interchanged. PARAM-6 is a graph of alpha on the ord- 
inate and omega on the abcissa. Magnitude is held at a fixed value through- 
out the calculations while constant beta curves are generated on the alpha- 
omega plane. This is the same arrangement as PARAM-3 with the roles of 
magnitude and omega being reversed. Interchanging the alpha and the beta 
in PARAM-6 will also cause a sharp alteration in the curve configurations. 
The familiar Bode plot is herein reproduced as PARAM-7. Needing no re- 
view, it is simply a logarithmic plot of magnitude on the ordinate, and 
omega along the abcissa. With beta held fixed at some constant value, 
constant alpha curves are generated. Conversely, when alpha is held 
fixed at some constant value, the beta curves are generated and plotted 
on the magnitude - omega plane. One problem exists with the logarithmic 
plot, however. In order to draw a useful logarithmic grid, it was neces- 
sary to shift the origin to the lower right hand corner of the graph paper. 
In this manner, the decade span could cover a total of 15 inches rather 
than a maximum of nine. The net effect has been to require the abcissa 
to be plotted along the right hand edge of the graph while the ordinate 
is plotted along the bottom where the abcissa otherwise would be, thereby 
coinciding with the X-axis. Thus the ordinate originates in the bottom 
right hand corner and increases toward the lower left hand corner, or 
original origin. The same process is required for PARAM-4 and PARAM 6 
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whenever specifying the logarithmic plots. Detailed coverage of the 
available options will be found in article 3.2 which is devoted to 
program PARAMS. 

The programs PARAM-D and PARAMS utilize only four variables, al- 
though their adaptation to others has been made. The method of solution, 
however, has been with two variables along the axes, while the remaining 
two have been employed as constant curves, and in the case of program 
PARAMS one variable has been ehld fixed at a constant value. In program 
PARAM-D, the axis variables have been selected as alpha and beta only; 
no attempt has been made to develop any other projection whereby other 
variables might replace alpha and/or beta as axis coordinates. Program 
PARAMS has attempted to focus attention on several two dimensional pro- 
jections in the hopes that peculiarities not evident in one planar 
presentation might become conspicuous in another mode. Due to the ex- 
tensive additional projections with the incorporation of one more vari- 
able, it was decided to refrain from including phase ( © ) as a para- 

meter in all but the Bode diagram, PARAM-7, where its utility and value 
have already been demonstrated. 

The following article will develop the basic capabilities of program 
PARAMS with regards to its application in frequency response studies. 

3.2 Digital Computer Program PARAMS. 

The preceding article and earlier example in section 2 presented in 
some detail the capabilities of program PARAM-D for stability studies. 
Program PARAMS has been developed to study some pertinent factors in fre- 
quency response which might prove useful in systems and controls design. 

For each of the graphical representations entitled PARAM-1 through 
PARAM-7, a separate computer program was developed in FORTRAN II language 
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for use on the Computer Data Corporation (C.D.C.) Model 1604 digital 
computer [11]. Each of these programs was compiled in a similar fashion, 
and employed the algebraic methods derived by Doctors G. J. Thaler and 
A. G. Thompson [12]. As previously discussed in this section, each graphi- 
cal configuration has been developed in order that the controls design 
engineer might gain additional tools to facilitate studies in analysis 
and synthesis. For the study of one particular problem, it is often 
necessary to analyze the systems behavior from many different angles, and 
then select a smaller number of characteristics for further intensive 
study. Unfortunately, every real life problem Involves a different set 
of constraints, and for that reason there can be no set pattern for suc- 
cessful design application in all cases. Many situations reduce to trial 
and error techniques on preliminary studies in order to obtain the best 
compromise when a number of possible solutions exist. The ability to ac- 
complish this objective using many separate computer programs is indeed 
possible but cumbersome. Instead of this method, one consolidated com- 
puter program titled PARAMS has been written in FORTRAN IV language for 
use on the International Business Machine (IBM) model 360 computer. The 
entire computer program, together with explanatory flowgraphs and illustra- 
tive printouts, has been Included as Appendix I. Details of its use, to- 
gether with the advantages and disadvantages, will now be presented. 

Versatility is indeed a virtue, provided it does not gain domination 
over simplicity. While considerable apprehension might be the first re- 
action upon reviewing Appendix I, a systemized pattern will become evident 
upon further inspection. With this program, complexity does not enter into 
the computational procedures, but rather a considerable amount of detail 
has been included in the printout format for the sake of clarity when re- 
viewing the material contained within the graphs. At most, the computational 



76 



process is intricate in the initializing stage wherein the various axis 
maximum and minimum limits, together with the curve and fixed constant 
values are redesignated for the ensuing common computational and graphi- 
cal construction sections. This can be verified by noting that the numera- 
tor and denominator coefficient calculation sections are common to PARAM-1 
through PARAM-6. It will further be seen that the gain computation for 
PARAM-7 employs different formulae but precisely the same logic, and for 
that reason is incorporated in the same numerator and denominator coef- 
ficient calculation sections as used for PARAM-1 through PARAM-6. Only the 
phase calculations for PARAM-7 require a different process for generation 
of the coefficients. For this reason, this computational section is by- 
passed for all other processes. 

The logical commencement point for design studies is to define a 
range of parameter values which guarantees all roots to lie in the left 
half plane, for clearly the question of stability must remain foremost in 
the designer's mind. Program PARAM-D has been devised to meet this demand, 
and its capabilities to provide a stability perspective have already been 
delineated. In determining the various possible frequency response 
characteristics for a given system, several graphical representations 
will be invariably studied. Most likely, however, program PARAM-7 will 
be employed in the final analysis as an overall check on system specifica- 
tion criteria. The remaining programs, PARAM-1 through PARAM-6 are most 
informative in the initial and intermediate stages of design study. It 
would therefore seem that the inclusion of all frequency response para- 
meter plane programs into one package would be most useful. Different 
design studies will normally emphasize certain performance criteria more 
than others. Thus seldom would the same sequence of graphical presenta- 
tions be called for in two consecutive studies. In fact, PARAM-7 might be 
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omitted for the initial coarse adjustments, but yet included in the final 
stages of refined improvements. In this manner each system study can be 
tailored to meet certain specifications and suit specific needs. Also, 
certain systems will be more complex than others. Therefore, more than 
one mathematical model might be considered for simultaneous inspection. 

The program is constructed to handle any number of such studies at one 
time. Each of these normally requires a complete set of data cards, and 
a brief explanation of their usage will follow. Additionally, an abbrevi- 
ated description section appears in the program preface, found in Appendix 
I. The total number of these data card sets must be indicated on the 
first card following the main program. This desired number of runs is 
designated by the variable NRUN. The program simply processes the re- 
maining data cards by commencing a ’'DO LOOP" with the variable IRUN going 
from 1 to NRUN. The first card of each data set indicates which graphi- 
cal presentations (PARAM-1 through PARAM-7) are desired. Again, there is 
no limit to the number of these. The program simply commences another 
"DO LOOP" at statement number 26, with IPARAM going from 1 to KPARAM. 

This cycles each of the graphical constructions in the order called for. 

In this manner, a great multiplicity of variations exist for the simultan- 
eous study of different preoblems. At this point it is pertinent to in- 
dicate a technique by which additional finesse may be utilized. Although 
most studies involve extensive rearrangement of a data set from one system 
to another, some require only a slight modification. In these latter cases, 
it is not necessary to submit an entire additional set of cards for each 
study. By first specifying the desired graphical plots of one system, it 
is possible to consecutively list the desired sequence for another system. 
Then a simple insertion of the desired data modification immediately follow- 
ing statement number 26 will alter the original data in the specified manner. 
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As an example, suppose one desired to first construct PARAM-5 with alpha 
on the ordinate, and then with beta on the ordinate. Simply list KPARAM 
as 2, and then indicate 5 twice for two consecutive plots of PARAM-5. 

Insert the following statement on an additional card behind statement 26: 

**1? (IPARAM. GE. 2) MODORD = 2*'. Thus the required change in the input data 
is effected at the desired time by one card rather than a complete set of 
data cards. 

Although no standard names have been assigned to these different pro- 
grams the title of each graphical representation (PARAM-1 through PARAM-7) 
has been retained to eliminate confusion in terminology for those persons 
previously exposed to the partial programs which are already in existence. 
The result has been a rather complex programming sequence for printout of 
graphical information and preliminary initializing computations, but as 
has already been stated, the basic manipulations are common to all pro- 
grams, and are centrally located in the program from statement number 350 
through statement 426. A cursory inspection of these 76 statement numbers 
as well as the remainder of the program will reveal that its construction 
is modular in form, and each section is a separate entity unto itself with 
adequate illustrative sub-titles. The overall program has been so ar- 
ranged that many sections can be removed if it is desired to employ their 
use elsewhere. Conversely, additional extensions are easily facilitated 
by this manner of program fabrication. 

Careful consideration has been used in this program design to retain 
flexibility for use on multiple problem types. Although the program is 
slightly larger in size than any of its component programs which were 
previously mentioned, it is smaller than the combination of any two of 
them. While intensive effort to streamline programs is not of particular 
significance, the repetition of computer processing is needless and 



79 



inefficient. Also, the preparation of data cards for a large number of 
programs serves only to increase the probability of error with consequent- 
ial waste in time, expense and labor. More important is the fact that 
computer programming is still a dynamic field which often necessitates 
program modifications to meet the requirements of more sophisticated equip- 
ment. Thus one might be motivated to use a single computer program in 
lieu of eight separate ones, and to this extent consolidation is warranted. 

To ensure adequate clarity of results, an elaborate printout of in- 
formation has been established. This information can be basically cata- 
gorized into three groups, each of which is controlled by a single input 
data variable IWRITE. The first reproduces all input data for verifica- 
tion, and is effected whenever the value of IWRITE is any non-negative 
number. The second category gives detailed information about how each 
graph will be constructed. Included in this section are the graph titles, 
the range of values for the ordinate and abcissa variables, the curves to 
be generated, and the parameter that is to be held constant throughout the 
computational process of that particular program. This explanatory sec- 
tion will be printed for all positive values of IWRITE. Finally, the com- 
puted results for each curve are printed out at pre-selected intervals. 

The value assigned to IWRITE determines at what intervals the computation- 
al results will be printed, and the range of permissible values is from 
one to 900. Thus an assigned value of one would result in the reproduction 
of every calculated curve value, whereas a value of ten would give a print- 
ed value for every tenth calculated point. Approximately one third of the 
program is devoted to printing out the above information for feedback pur- 
poses. If desired, these sections can be partially or entirely bypassed 
by the appropriate selection of the input data variable IWRITE. Any nega- 
tive value assigned to this variable will provide complete calculation. 
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but without any printout whatever. 

It has already been pointed out that if the roles of alpha and beta 
are allowed to be interchanged, a much smaller number of graphs are re- 
quired, In order to accomplish this objective, some criterion must exist 
for the choice of which variable (alpha or beta) will be plotted on the 
ordinate in the case of PARAM-1 through PARAM-6, and which variable will 
be held at a fixed constant value while the other generates the desired 
curves in the case of PARAM-7. The input data variable titled MODOEU) has 
been designated for this purpose, and so by this selection the type of 
curves is controlled. With a value of one for MODORD, alpha curves will 
be generated in the case of PARAM-7, and alpha will be plotted on the 
ordinate for the rest of the programs. Conversely beta will be used in 
lieu of alpha for those functions when the value of MODORD is set equal 
to two. This flexibility for PARAM-1 and PARAM-2 might be unnecessary, 
because as indicated earlier, a ninety degree reorientation of the graph 
is the only effect of interchanging alpha and beta when these variables 
are on the axes. However, the present day graphical plotting facilities 
are generally constructed in such a manner that the abcissa, or X-axis, 
is more severely restricted in length than the ordinate, or Y-axis. The 
facilities available at the inception of program PARAMS were in fact such 
that the ordinate length could be approximately twice as great as the 
length of the abcissa. It was envisioned that in some cases an expansion 
of one variable might be more desirable than another. Extending the ap- 
plicability of MODORD to PARAM-1 and PARAM-2 makes this provision possible, 
which otherwise would require modification of the main program to achieve 
such results. Furthermore, it is advantageous to retain a common orienta- 
tion of graphical construction when viewing several projections simultan- 
eously. 
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The value of magnitude often occurs in several different forms, and 
thus it provides a source of confusion if the terminology is not clearly 
specified. The input data variable IMAG serves to eliminate this poten- 
tial error by distinctly categorizing magnitude into one of three forms. 

The program is constructed in such a manner that the form of magnitude as 
specified by IMAG is maintained consistently for all graph plots. Al- 
though this is strictly a matter of convenience, it was found that various 
approaches to design require magnitude to be handled in different forms; 
hence this inclusion. Thus for magnitude to be interpreted as decibels, 
a value of zero would be assigned to IMAG, and the curve values of PARAM-2 
and PARAM-4 would be listed in decibels, while the abcissa values of PARAM- 
3 and PARAM-5 would likewise be interpreted accordingly. If actual magni- 
tude values were to be used instead, IMAG would necessarily have to be 
changed to one. Assigning a value of two to IMAG would dictate all values 
to be handled as magnitude squared. In all facets of life it seems that 
most conventions are fraught with exceptions, and this program has not 
escaped that dilemma. Since all Bode plots associate magnitude with deci- 
bels, it was decided to maintain this practice; consequently the selection 
of IMAG has no effect on the form of PARAM-7 graphs. They automatically 
plot magnitude along the ordinate in the form of decibels. Conventional 
practice has also dictated that omega be logarithmically plotted along the 
abcissa for all Bode diagrams. Other frequency response graphs which also 
plot omega on one axis have not acquired any standardized form, and for that 
reason provisions have been made to select a square linear grid or a logari- 
thmic grid as deemed appropriate. The input data variable IGRAPH provides 
this option for the various possible combinations of PARAM-4, PARAM-6 
and PARAM-7. Additionally, this variable is used to prevent any graphs from 
being generated by assigning to it a value of minus two. This provision is 
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useful only for efficient use of available hardware when desired graph 
scales are initially unknown. 

In order that maximum versatility might be gained, the memory stor- 
age requirements of this program have been expanded to present day capa- 
bilities. A review of the program DIMENSION statement just preceding the 
input data read in section will show that the bulk of the storage capa- 
city is consumed by the two dimensional variables S through W. 

These are employed for the various optional arrangements of PARAM-7 and 
may be combined into one common storage array for all other calculations 
when employing this program in smaller computer facilities. In this 
manner, the memory core may be reduced from 500K bytes to less than 200K 
bytes. This is accomplished by inserting the following statement after 
the dimension statement: EQUIVAIANCE (SjTjUjVjW) . 

The additional variable, phase, has been incorporated in PARAM-7 
with several variations possible. The input data variable IPHASE is 
used to govern the relation between the magnitude and phase plots. 

No phase curves will be drawn when IPHASE is set equal to a negative num- 
ber. Setting IPHASE equal to zero will enable phase curves to be super- 
imposed upon the magnitude curves for ease in investigating stability 
by means of phase and gain margins. Assigning any positive number to 
IPHASE will provide separate phase curves. To allow for various super- 
positions, and various separate phase curve configurations as well, it 
was decided to employ an additional variable IQUAD which allowed the phase 
graphs to be plotted in the following two manners: When IQUAD has an as- 
signed value of zero, phase curves will be automatically scaled at 45 
degrees per inch, with 180 degrees on the magnitude axis. When the vari- 
able IQUAD is positive, however, the phase curves will be plotted with a 
scale equal to the value of IQUAD, starting with 0 degrees at the top of 
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the graph. In this manner, it is possible to cover any desired range of 
phase angles. 

Basic conformal mapping theory indicates that a polar plot can be 
generated from the information contained in the magnitude and phase plots. 
Particular importance might be attached to the polar plot representation 
of a system when the Nyquist criterion is preferred over the Bode stabil- 
ity analysis method. Although it is a duplicative process and therefore 
not strictly required for information content, its inclusion has been 
assumed appropriate for the sake of completeness. The basic construc- 
tion of the polar plane and polar curves are contained within the program 
statement numbers 537 and 454 under the heading "PARAM-7 Polar Plot Curves”. 
Its usage is governed by the rather appropriate input data variable IPOLAR, 
and a polar plot will be drawn when its value is one. Conversely, no 
polar plot will be drawn for any other value of IPOLAR. 

An additional provision has been incorporated in program PARAMS 
for extended usage of the Bode plot. Certain advanced feedback design 
techniques call for complex systems to be studied from a cumulative stand- 
point. Such design measures require vector addition of component Bode 
diagrams and as such, the linear addition of each magnitude plot and each 
phase plot will yield erroneous results. Instead, the vector quantities 
of each design segment are summed in complex notation, and the final re- 
sult is then plotted as a single magnitude curve, and also a single phase 
curve when the latter is requested. The program subsection designated to 
accomplish this process is titled ”PARAM - 7 Superposition Section”, and 
is contained by the program statement numbers 520 and 529. It is con- 
trolled by the input data variable NPAR7 , where the magnitude of this 
variable dictates how many component terms are desired to be added toget- 
her. It is necessary to emphasize at this point that this section performs 
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a vector superposition of feedback terms, and not merely a scalar super- 
position. Thus whenever NPAR7 is greater than one, this section will 
add NPAR7 minus one terms vectorially to the first term before graphing 
the results. In order that these components may be properly calculated, 
provisions have been made to insert the component equation terms in the 
same coefficient format as normally submitted for the conventional cal- 
culations. These additional terms, however, are appended to the initial 
group of input data cards. Therefore NPAR7 minus one groups of coeffici- 
ents must be submitted directly behind the title cards of the rest of the 
program. 

A brief summary of the basic operational procedures which are common 
to each of the programs will now be undertaken. The schematic corollary 
in appendix I will parallel the following outline and should be consulted 
together with the main program for more detail. 

The first six statements in the "Input Data Read in Section" pro- 
vide ^in proper FORTRAN language, the necessary format for most of the 
read and write statements of the program. Their use is specified in the 
"Data Card Information Section" which appears at the beginning of the 
program. The remainder of the input data read in section is devoted to 
reading all the 19, 21, or 23 data cards as outlined in the preceding in- 
formation section. The next step is to print out the input data when so 
requested for verification. This procedure is accomplished by the section 
bounded by the statement numbers seven through 23, inclusive. Following 
the input data verification section, the program initializes the various 
parameters and computes the necessary controlling variables from such in- 
put data as axes maximum and minimum values. Within this section, con- 
tained by statements 24 to 54, the desired value of NPARAM is selected and 
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this value guides the rest of the program sequence by assigning the proper 
axis variables and their limits, the correct number and type of curves to 
be generated etc. When the logarithmic plot is requested for PARAM - 4 
or PARAM - 6, and whenever the Bode plot is to be generated, the section 
titled "Log Grid Graphical Construction" is executed. Extending from 
statement number 55 to 58, this sub-section utilizes the information 
generated in the initializing section to determine the ordinate scale 
and origin shift if applicable. Then it calculates the required number 
of decades by rounding up the maximum limit of omega to the nearest power 
of ten, and rounding down the minimum omega value to its nearest power of 
ten. Utilizing these factors and the axes lengths, a logarithmic grid is 
constructed in such a manner that the omega values are properly labelled 
along the right hand edge of the graph paper, while the ordinate values 
are indicated across the bottom in such a manner that their magnitudes are 
correct but the signs reversed. Thus the ordinate scale must be correct- 
ed to read increasing from right to left. Although this arrangement is 
less than optimum, it has proved to be much better than some alternatives 
which were previously attempted. 

The most prolific section of the program ranges from statement number 
59 to 200 under the sub heading "Graphical Information Printout Section". 
Its utility lies in the desire to establish positive correlation between 
printed data and graphical presentations when several similar graphical 
analyses are computed simultaneously. Following the extensive information 
printout, the main computational process commences by initiating calcula- 
tions for the first curve with statement number 202. Next, the coefficient 
preceding the quadratic formula square root is assigned by the variable 
ZSIGN. Following this designation, the abcissa is incremented 900 times 
from the minimum specified value to the maximum value. Additionally the 



86 



curve value is selected, based upon the value of the variable in statement 
number 202. Then the appropriate parameter which is to be held at a 
fixed constant value is assigned its pre-selected value. The ensuing 
calculations commencing with statement number 350 are in conformance with 
the theory presented in section 2. The magnitude protion of the Bode 
program requires the coefficients designated SUMN for the numerator and 
SUMD for the denominator, while the remaining programs utilize the 
numerator variables SNB2, SNBl, and SNBO in addition to the denominator 
variables SDB2, SDBl, and SDBO. Such requirements are generated under 
the sub headings ''Numerator Coefficients Calculations" and "Denominator 
Coefficient Calculations". Following these sections, the necessary phase 
coefficients are computed from statement number 386 through 400, and the 
additional Bode curve ordinate value calculations are continued through 
statement number 410. The section titled "Quadratic Formula Solution of 
Dependent Variable, PARAM-1 through PARAM-6" is self explanatory and 
simply computes the ordinate value cwnmensurate with the assigned fixed 
value of the parameter which is kept constant, together with the specified 
value of the curve which is being generated, and the corresponding point 
on the abcissa. The computed results are then printed in accordance with 
the provisions contained between statement number 426 and 448. Having 
computed this ordinate value, the abcissa which is the independent vari- 
able is again incremented and a corresponding ordinate, or dependent vari- 
able point is calculated. After reaching the maximum permissible abcissa 
value, the process is again repeated for the alternate solution of the 
quadratic formula. When both possible ranges of values have been scanned, 
the next curve is specified, and the entire procedure repeated. Upon ter- 
minating the curve requirements, the program then reverts to statement 
number 26 for the next cycle in the run. When all requested graphical 
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construction cycles have been accomplished, the next run is processed 
in like manner. The overall program is terminated when all desired runs 
have been computed and graphed. 

The graphical plotting procedures will be discussed separately be- 
cause their intricate operations tend to destroy the continuity of the 
preceding discussion on the overall program operation. Basically, the 
graphical plotting subroutine is executed whenever the following condi- 
tions are engaged; 1) upon initially setting up the logarithmic plot 
when specified in accordance with the above dictates; 2) when two or more 
computed points fall within the prescribed graph range and the next point 
lies outside the ordinate graph limits; 3) when the end of each scan of 
the abcissa has been reached - once for each value of the quadratic equa- 
tion square root coefficient; 4) at the end of each set of curves for one 
graph, as accomplished by the graphical termination section found at the 
end of the program beginning with statement number 999. 

The versatile arrangement necessitates several provisions for the 
many possible combinations of the Bode plot. Since the quadratic formula 
is not employed in these calculations, only one incrementing scan of the 
abcissa is performed. Additionally, the computed magnitude values have 
been programmed in such a manner that when the ordinate graph limits are 
exceeded, the external points are fictitiously redefined to lie on the 
appropriate ordinate graph limits. This procedure was necessitated by the 
additional phase calculations which remained within the graph limits 
throughout the entire range of the prescribed omega values. Finally, a 
separate provision had to be made for the construction of the phase curves, 
since an additional logarithmic graph was required when phase curves were 
plotted separately. This procedure has been accordingly programmed in the 
section titled "PARAM - 7 Phase Curves" found between statement numbers 529 
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and 537. The section titled "PARAM - 7 Polar Plot Curves’" has already 
been discussed and is a separate entity unto itself, utilizing the pre- 
computed values of magnitude and phase, and plotting same in polar coordi- 
nates. This process is accomplished between statement numbers 537 and 999 
in the program. 

The preceding development has summarily shown that this program con- 
tains broad applications, and is well suited for additional modifica- 
tions to meet particular needs. Blanket provisions for another graphical 
configuration has been purposely included in the program sequence labelled 
NPARAM- 8, and should future expansion be deemed appropriate, the required 
modification will merely entail filling the appropriate spaces with the 
desired initializing variables and applicable computational equations. 

To facilitate such extensions and also to enable maximum comprehension 
of the program operations, two logic flow diagrams have been included in 
Appendix I. The first is a generalized sequentail listing of logical 
operations which are applicable to all graphical configurations. The 
second is a detailed flow chart for PARAM - 6, so chosen because it in- 
cludes most of the alternatives presented in this discussion. 

A final comment is appropriate concerning the application of computer 
program PARAMS. As indicated earlier in this article, considerable atten- 
tion has been given to versatility. The associated price of computational 
time and memory storage has consequently been paid. The procedures to be 
used in reducing memory storage for PARAM-1 through PARAM-6 has already 
been indicated. Compilation time is roughly 1.5 minutes. An additional 
30 seconds is consumed for each printout of preliminary graphical informa- 
tion. The actual calculation of 16 curves for each graph requires approxi- 
mately 3 minutes. Thus a considerable amount of computer time is consumed 
when a complete set of curves is requested for all possible graphical 
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representations. Various methods exist to reduce this time requirement. 

By submitted an object deck, compilation time is deleted. This imposes 
the restriction of submitting a complete data card set for each run, 
however. Eliminating the input data printout, and the preliminary graphi- 
cal information printouts will also reduce the time requirements some- 
what. This, however, is not recommended because the elaborate informa- 
tion presentation provides an exceedingly valuable check on the validity 
of the graphical results. Finally, the most time is consumed by a highly 
sensitive incrementation of the abcissa variable. It will be recalled 
that the range of the abcissa variable is divided into 900 evenly spaced 
steps. For each of these abcissa values, a corresponding ordinate value 
is computed. Although many situations lend themselves to a larger in- 
crementation, it was decided to maintain the fine step sizes for highly 
sensitive system curves. The maximum sensitivity of 900 points per graph 
line is the restriction placed upon the program by the CALCOMP plotter in 
use at the Naval Postgraduate School computer facility. Reduction of this 
incrementation sensitivity is therefore another means by which computation- 
al time can be further reduced. The associated loss of information con- 
tent must be considered in such cases. As will be pointed out in the 
next section, not all PARAMS frequency response curves would be required 
or desired in every design study. The judicious choice of the most valu- 
able plots for a given set of constraints will normally eliminate approxi- 
mately half of the program possibilities. Only by experience, insight and 
intelligent design procedures can a good trade-off be achieved between the 
tools available and the data required. 
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4. 



Tachometer and Acceleration Feedback Compensation Studies. 

The mathematical developments of section 2 have provided the basis 



for constructing the computer programs which were discussed in Section 3. 
This section will apply these methods to a system study employing tacho- 
meter and acceleration feedback as compensation techniques. The frequency 
response parameter curves will be analyzed according to the respective 
utility of each as a design tool. Meticulous detail is employed in the 
initial discussion of these curves. This is to provide a clear under- 
standing of the limitations and advantages of each respective graphical 
representation . 

4.1 Third order 

Although an entire library of control engineering studies has been 
created, there has been no attempt to distinguish a set of example prob- 
lems as superlative for generalized case studies. Therefore a simple 
third order system has been considered for an initial application of the 
parameter plane frequency response curves. The illustrative example which 
was introduced in section 2 will now be explored in depth. It will be re- 
called that Figure 2-3 presented the schematic diagram of the third order 
system with compensation included. The stability curves for the system 
where mapped into the parameter plane shown in Figure 2-5. In that pre- 
sentation alpha, the abcissa variable, represented acceleration feedback 
( K CL )• Similarly beta, the ordinate variable, indicated tachometer 
feedback ( ). Reference to the characteristic equation of the sys- 

tem shown in Figure 2-3 shows that the steady state velocity error is 




500 + K Kt 
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Inspection of shows that the smallest allowable value of beta 

will minimize the steady state velocity error. It is presumed that this 
is one of the predominant performance criteria. Furthermore, one can 
arbitrarily stipulate the dominant complex roots to have a damping ratio 
( ^ ) of 0.5. With these system stability requirements, one can see 

that optimum performance will occur at point A in Figure 2-5. Thus the 
optimal value of beta is 0.036, while the optimal value of alpha is 
0.00047. 

A tacit assumption up to this point has been that frequency perform- 
ance constraints are of minor importance. Most design studies do not 
conform to this type of ommission; however stability is of interest in 
all systems. It has therefore been assumed that the two stability per- 
formance criteria above are desirable but not mandatory. This permits 
frequency response specifications to be considered with some initial 
reference point. Figure 4-1 shows an approximate correlation between 
frequency response, transient performance, and root locations in the com- 
plex frequency plane. Additionally the separate root loci with alpha and 
beta varied in Figures 4-2 and 4-3 respectively show the approximate vari- 
ation of root locations as a function of one parameter only. The corres- 
ponding effect in the frequency domain is represented by the Bode plots 
in Figures 4-4 through 4-7. The graphical relationships presented in 
Figures 4-1 through 4-7 are certainly not new to the design engineer. 
Rather, they are the fundamental tools used in most system studies today. 
As such, they are included primarily to establish a firm foundation upon 
which parameter plane frequency response curves can be built. In this 
manner, a better understanding can be gained concerning the utility of 
the latter curves which are unfamiliar to most engineers. 
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PARAM-1 



Complete analysis of the information contained in PARAM-1 plots 
will now be endeavored. Figure 4-8 shows the relationships between alpha, 
beta and omega which exist for a constant magnitude of 12 decibels. It 
is clear from the legend at the bottom of the graph that alpha is the 
abcissa variable, while beta values are indicated on the ordinate. Each 
curve in the alpha - beta plane represents a fixed magnitude of 12 deci- 
bels, but a different value of omega. Thus the curve of W =10 rad/sec 
barely appears in the lower right hand corner. As omega increases, the 
curves can be seen to shift toward higher values of beta and lower 
values of alpha. It can also be seen that the area enclosed by each 
curve decreases as omega increases to = 70 rad/sec. The meaning of 
each of these 14 curves can be explained as follows: For the third order 

system in question, any point on a specified omega curve depicts the 
locus of alpha and beta values which will guarantee a magnitude of 12 
decibels at that specific frequency. This by itself does not portray 
much information, because it says nothing about the area enclosed by the 
curves, nor the area outside the locus of these curves. In most situations, 
it is possible to glean considerably more information out of these curves. 
Consider for a moment that value of alpha and beta at which two omega 
curves intersect. Point B of Figure 4-8 shows one intersection of the 25 
rad/sec and UJ = 30 rad/sec curves. For this value of alpha and beta, a 
magnitude of 12 decibels is guaranteed to exist at both values of omega 
indicated. This observation is important, for it means that the magni- 
tude is either greater than 12 decibels, or less than that value for the 
intermediate range of omega. Furthermore, it is possible to assume that 
the magnitude cannot be both greater than and less than 12 decibels for 
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these intervening values of omega. The reason for this is simply that 
the magnitude curve is continuous, and must pass through 12 decibels in 
order to change from a high value to a lower value, or vica versa. If a 
sufficiently large number of omega curves are drawn in the region of 
interest, it can be seen that, at Point B, only two omega curves do actual 
ly intersect. Therefore, this range of omega values represents the 12 
decibel bandwidth of a system which satisfies these values of alpha and 
beta. With this insight, it can be concluded that a larger bandwidth is 
obtained for a point slightly below point B« Thus at a value of alpha 
and beta indicated by point C, interpolation of the omega curves shows 
an approximate bandwidth from CO = 24 to (jJ =31. Conversely, preced- 
ing in the opposite direction, the bandwidth is reduced to the point D 
at which only one frequency curve satisfies the specified value of alpha 
and beta. No other frequency curves will intersect this point, for it 
lies on the tangent locus of all frequency curves evaluated at 12 deci- 
bels magnitude. Thus this tangent locus defines all the values of alpha 
and beta which guarantee a peak resonant magnitude of 12 decibels. The 
value of the resonant peak frequency is not specified by the tangent locus 
The reason for this is because omega is an implicit function of the tan- 
gent locus equation. Article 2.5 discussed various methods to generate 
the required non-parametric form of this p Lio) equation, and the 
limitations encountered in each of those approaches. For those reasons, 
the tangent locus curve was not constructed. Instead, a sufficiently 
large number of omega curves were constructed to show the approximate loca 
tion of the tangent locus. Although the resonant peak frequency is not 
directly obtainable from the tangent locus, it is clearly seen that such 
information is readily found by the intersection of a particular omega 
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curve with this tangent line. 

The discussion thus far has presumed that all values inside the 
tangent locus represent peak resonant points greater than 12 decibels. 

This is intuitively obvious, but no proof has been made of this fact. 
Figures 4-9 through 4-12 have been included to verify this presumption, 
and also to illustrate additional phenomena. Figure 4-9 clearly shows 
that when the magnitude is reduced to 6 decibels, a greater area is en- 
closed by each of the omega curves. This results in a tangent line rep- 
resenting the locus of peak resonant magnitude equal to 6 decibels. 
Comparison of this tangent locus with the one in Figure 4-8 will reveal 
that the 12 decibel tangent locus falls within the area bounded by the 
omega curves of Figure 4-9. The points B and D of Figure 4-8 have been 
reproduced on the constant 6 decibel graph for comparison. It can be 
noted that point B represents a 12 decibel bandwidth of 5 rad/sec and a 
6 decibel bandwidth of approximately 13 rad/sec. Additionally, point D 
depicts a resonant peak magnitude of 12 decibels at 28 rad/sec with a 6 
decibel bandwidth of 13 rad/sec. Similar information can be obtained 
for other values of alpha and beta throughout the parameter plane. By 
similar analogy, it is possible to determine what range of alpha - beta 
values will guarantee a peak resonance not to exceed a pre-selected value. 
Thus all values between the M ^ - 12 and the M p =s 6 curves 

will guarantee all magnitudes to be less than 12 decibels, with a 6 
decibel bandwidth of variable size depending on the choice of alpha and 
beta. A closer inspection of Figure 4-9 will reveal that the maximum 6 
decibel bandwidth attainable is approximately 17 rad/sec for the situation 
just described. This maximum bandwidth occurs approximately alpha equal 
to 0.0001 and beta equal to 0.018. At this point the resonant peak magni- 
tude is slightly less than 12 decibels, with a resonant angular frequency 
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of approximately 40 rad/sec. The cut-off frequencies occur at CU - 32 
rad/sec and 60 =49 rad/sec. All of the above information has been 
derived from the evaluation of two PARAM-1 graphical plots. It will be 
recalled that the PARAM-7 curves of Figures 4-4 through 4-7 merely por- 
trayed Bode curves as a function of one parameter only. These PARAM-1 
curves denote a portion of many Bode curves as a function of two system 
parameters. The value of these two is restricted to only two magnitudes, 
but certainly the analysis of any two such graphs provides a certain in- 
sight as to the system behavior for other magnitudes. Interpolation can 
be used for the intermediate magnitude range. Such a procedure is not 
recommended without some caution for extensions of magnitude studies out- 
side the graphed limits. The consideration of Figures 4-10 through 4-12 
will illustrate this point more clearly. For a constant magnitude of 3 
decibels. Figure 4-10 shows a dramatic increase in the applicable range 
of alpha - beta values. The graph scales of all PARAM-1 plots were kept 
the same simply to illustrate the large increase in area encompassed by 
each constant frequency curve as the value of magnitude is decreased. 

As might be anticipated, a much larger 3 decibel bandwidth exists for 
the points, B, C and D shown on Figure 4-8. The bandwidth of these three 
points do not vary greatly, however. Each is approximately 20 rad/sec. 

This illustrates the sensitivity of beta, the tachometer feedback in the 
area under discussion. For a given value of alpha equal to 0.008, a 
slight change in beta greatly affects the size of the peak reasonant 
magnitude. Conversely, an equal change in alpha reflects a vast difference 
in omega, the location of this peak resonant magnitude, plus a slight 
change in its magnitude value. Although the sensitivity of the alpha vari- 
able might not seem apparent by looking at the curves, it must be remembered 
that there is a different graph scale for each of the axis variables. 
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Figures 4-11 and 4-12 show most clearly why the analysis of merely 
two PARAM-1 frequency response representations might be misleading. 
Consider point B of Figure 4-8 which is reproduced on the zero decibel 
magnitude graph of Figure 4-11. The constant omega curves of this plot 
show that no two curves intersect at point B. Clearly, it is seen that 
all values of omega less than 40 rad/sec encompass this choice of alpha 
and beta. Therefore a cut-off frequency of approximately 39 rad/sec is 
obtained for the zero decibel bandwidth of that system. Similarly, 
from Figure 4-12, the -3 decibel bandwidth is about 43 rad/sec. Using 
only these latter two PARAM-1 plots, one cannot discern any resonant peak 
for the system frequency response associated with these values of alpha 
and beta. It is only through the apriori information portrayed in any 
one of the Figures 4-8 through 4-10 by which we are able to make such a 
statement that a resonant magnitude actually exists. Only by the combined 
interpolation of all the Figures 4-8 through 4-10, can we state with a 
fair degree of confidence that the value of this resonant magnitude is 
approximately 15 decibels. This of course could be checked by construct- 
ing another PARAM-1 curve, or perhaps a Bode plot. Neither of these are 
necessary for initial coarse approximations. It will soon be shown that 
the combined usage of PARAM-1 and PARAM-2 provides a much better method. 
The reason for the inclusion of so many PARAM-1 curves was simply to il- 
lustrate how a valuable insight can be obtained for system frequency re- 
sponse studies when considering two parameter variations simultaneously. 
Within the limitations just indicated, it is possible to determine how 
two parameter variations affect bandwidth and peak resonance from a care- 
ful inspection of several PARAM-1 graphs. 

Before preceding to PARAM-2 graphical configurations, it will be 
advantageous to review Figures 4-8 through 4-12 in a three dimensional 
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concept. Consider for a moment a curved conical surface similar in 
shape to a bull’s horn. Also, let this curved cone be elliptical in 
shape, rather than round. Next assume that the base of this surface is 
in the lower right hand corner, nearest to the observer. The smaller end 
of this surface would therefore be in the upper left hand corner of the 
observation sector, at a certain distance away from the observer. Now 
superimpose the omega axis on the centerline of this conical surface in 
such a manner that the low values of omega are at the base of this sur- 
face, or nearest to the observation point. The higher values of omega 
would consequently lie further away from this observation point, in- 
creasing in value toward the curved conical tip. With this three dimen- 
sional surface concept in mind, look at Figure 4-8. The curved omega 
axis originates in the lower right hand corner and extends through the 
alpha - beta plane, curving upward and slightly to the right. Figure 4-9 
depicts this curved conical surface as larger in size, (i.e,, closer to 
the observer), with the base of this surface shifted slightly to the left, 
thereby bringing it more into the line of sight. Bringing this curved 
conical surface still closer to the observation point makes it appear 
larger, as seen in Figure 4-10. The zero decibel PARAM-1 graph. Figure 
4-11, portrays an even larger shape with the base shifted even more to 
the left. Finally Figure 4-12 can be regarded observing this curved sur- 
face with its base very close, similar in appearance to passing through 
a curved tunnel under a river. 

The preceding illustrative example can be considered as elementary, 
but is fundamental to the understanding of the various graphical rep- 
resentations to follow. No attempt will be made to correlate the il- 
lustration to the PARAM-7 curves of Figures 4-4 through 4-7, since those 
are essentially single parameter variations of the series of two parameter 
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surfaces just described. Better use of the illustrative example is made 
by correlating PARAM-1 with PARAM-2. 

PARAM-2. 

Consider now a cross sectional cut of the curved conical surface 
described above. PARAM-2 effectively maps constant magnitude curves onto 
the alpha - beta parameter plane for a fixed value of angular frequency. 

In this manner, it is possible to show how a particular omega value, 
such as to = 10 rad/sec will vary as a function of magnitude in the 
alpha - beta plane. Consider again the U) = 10 rad/sec curve in Figures 
4-8 through 4-11. The associated PAE^-2 plot is shown in Figure 4-13. 
Since the tU = 10 rad/sec curve represented a cross section of the curved 
conical surface near its base, it is possible to observe how this base 
is rotated toward the left as magnitude is decreased. Further inspec- 
tion of Figure 4-13 shows how the curve increases in size with a reduc- 
tion of magnitude. The -3 decibel, CJ = 10 rad/sec curve is in fact 
shifted completely off the graph of Figure 4-13. This effectively shows 
that no alpha - beta value within the prescribed range will have a magni- 
tude of -3 decibels at an angular frequency of 10 rad/sec. No information 
is conveyed about any other angular frequency in Figure 4-13. To obtain 
such information would necessitate a different PARAM-2 plot, or a compari- 
son with one or more PARAM-1 curves. This does not preclude the correla- 
tion of other frequency response curves. They can be consulted for pert- 
inent relationships, and will be discussed later. Figure 4-14 is a graphi- 
cal representation of constant magnitude curves for a fixed angular fre- 
quency of 20 rad/sec. Inspection of these curves shows that they lie very 
closely to the centroid of the curved conical surface discussed in the il- 
lustrative example of this article. This can be verified by their 
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concentric shape. Due to the graphical construction technique, these 
curves do not extend to the bottom of the graph. This can be explained 
by the fact that these curves are almost vertical at the lower graph 
edge. Even though the abcissa variable is incremented 900 times during 
each curve computation, this is not sensitive enough to allow the detec- 
tion of such a steep slope. If an infinitessimally small step size were 
possible, these curve values would be detected. This limitation of the 
graphical plotting procedure must be kept in mind. The mathematical 
theory of section 2 indicated that two quadratic formula solutions were 
necessary for the complete solution of any curve generated by PARAM-1 
through PARAM-6. The computer technique discussion of section 3 dis- 
cussed how each graph curve was constructed first for one sign of the 
square root coefficient in the quadratic formula, and then for the other. 
Subsequent graphs will show how these two curve solutions do not always 
connect. The reason for this again lies in the finite incrementation of 
the abcissa variable. Returning once again to Figure 4-14, it is possible 
to note certain characteristics regarding the angular frequency of 20 rad/ 
sec as a function of the two parameters alpha and beta. The intersection 
of the peak values of each magnitude curve defines a vertical line in the 
parameter plane. Therefore it is possible to conclude that variations in 
magnitude for an angular frequency of 20 rad/sec are highly sensitive to 
beta values, and quite insensitive to changes in alpha. This says nothing 
about resonance magnitudes. Rather, each curve guarantees a particular 
magnitude value for a given angular frequency. Also, each curve defines 
an alpha - beta region within which the magnitude will always be greater 
than the boundary value for that value of omega. As a further analysis of 
the system response for point B of Figure 4-8, it can be seen that such a 
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system will have a magnitude of 6 decibels at a frequency of 20 rad/sec. 
Figures 4-15 and 4-16 show considerably more information than the previous 
two PARAM-2 graphs, because the center of these eliptical curves is indi- 
cated. This center corresponds to the centerline of the three dimensional 
curved conical surface which has been repeatedly referenced up to this 
point. It should be observed that the center of these ellipses in Figures 
4-15 and 4-16 denote the maximum magnitudes attainable at their respec- 
tive angular frequencies. Thus by analogy, the centerline of the three 
dimensional surfaces in Figures 4-8 through 4-16 serves a dual purpose. 

Not only does it specify a particular angular frequency, but it denotes 
the locus of the highest peak resonant magnitudes for any alpha - beta 
combination of the system? This is an important conclusion, because it 
indicates whether a desired peak resonant magnitude value can ever be 
reached for the system under consideration. The comparison of PARAM-1 
and PARAM-2 graphs further shows how the conical surfaces expand about 
this centerline with decreasing magnitude. From these, it is possible 
to visualize the locus of constant magnitude conical surfaces, each 
surrounding the centerline just discussed. It is therefore seen that a 
reduction of tachometer gain coupled with an increase in acceleration 
gain will shift the resonance peak to a higher frequency. Point B of 
Figure 4-8 has also been reproduced on Figures 4-15 and 4-16. These 
show that a system with the indicated alpha - beta values will have a 
nagnitude of 12 decibels at 6 l) = 30 rad/sec, and -1 decibels at an angular 
frequency of 40 rad/sec. At this point, enough information has been ob- 
tained to obtain a fairly accurate Bode plot for the parameter values given 
at point B. Further comparison of the 15 decibel magnitude curve in 
Figures 4-14 through 4-16 show that point B does in fact have a magnitude 
equal to 15 decibels. This confirms the assumption of this value based on 
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interpolation of the PARAM-1 curves alone. Since enough information has 
been depicted in the PARAM-1 and PARAM-2 curves to portray a Bode plot 
of point B, it is obvious that similar Bode plots could be obtained for 
any other desired combination of alpha-beta points. This is exceedingly 
valuable, because it allows the design engineer to construct any number 
of coarse Bode plots by evaluating the data contained in Figures 4-8 
through 4-16. Therefore it is concluded that the combined utilization 
of PARAM-1 and PARAM-2 is a fundamental starting point for all two para- 
meter frequency response design studies. 

The preceding discussion concerning the use of PARAM-1 and PARAM-2 
was very comprehensive because of the following three key factors; First, 
a thorough understanding of new concepts can only be gained by a slow and 
detailed presentation of the uses and limitations of such innovations. 
Secondly, the combined application of PARAM-1 and PARAM-2 plots provides 
the broadest information spectrum for frequency response design studies. 
Finally, a thorough understanding of these two parameter plots is pre- 
requisite to visualizing the following frequency response projections 
which are essentially a function of one system parameter. The remainder 
of this article will cover the various roles played by PARAM-3 through 
PARAM-6. 

PARAM-3 . 

Section 3 developed the major graphical presentations which were 
considered to be of most use to design engineer. In that discussion, 
it was indicated that two basically different configurations were avail- 
able in PARAM-3. Figure 3-1 showed in schematic form these differences. 
Attention will first be focused on the M0D0RD=1 variety as shown in 
Figure 4-17. Before discussing the information content of this graph, 
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an apparent mistake in the graphical plotting procedure is noteworthy. 

This set of curves was purposely offset to illustrate a limitation en- 
countered by careless selection of axis variable ranges. Graph scales 
are limited to one significant figure. As such the automatic scaling of 
the program PARAMS does not always produce the values which might be ex- 
pected. The requested range of alpha values was from -0.001 to 0.007. 
Since the ordinate length was specified as ten inches, the program comput- 
ed a grid scale of 0.0008 units per inch for alpha. Accordingly this 
scale was used; however, it placed the origin one inch above the minimum 
value in accordance with program plotting limitations. Thus the requested 
range of alpha values was plotted properly, even though the choice of 
ordinate values was improperly chosen. The preceding was for information 
purposes only; PARAM-3 plots will now be studied. 

Figure 4-17 portrays a set of constant beta curves on the magnitude 
- alpha plane for a given fixed angular frequency value of laJ = 30 rad/ 
sec. This graph can best be understood by comparing it with the PARAM-2 
plot of the same constant omega, namely Figure 4-15. Effectively, the 
PARAM-3, M0D0RD=1 curves take the information contained in the horizontal 
cross-sections of PARAM-2 and map this data onto the magnitude - alpha 
plane. This can be most easily visualized by taking the base of Figure 
4-15 and comparing it with the beta = 0.0 curve in Figure 4-17. With 
these two viewpoints, one readily sees another facet of PARAM-2 plots. 
Consider for a moment a topographical map of land terrain. Such maps 
represent constant altitude as a continuous curve. Each such curve is 
always a closed contour, with either peaks or valleys contained within 
the curve boundaries. PARAM-2 is such a topographical plot, with magni- 
tude equal to elevation. This concept is borne out by the magnitude- 
alpha plane of Figure 4-17. The constant beta curves of this latter 
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graph can be regarded as vertical slices into the topographical hill- 
side. As such, one sees which cross sections have actually been sliced 
out of the PARAM-2 magnitude contours. The complete range of PARAM-2 
alpha and beta values has been spanned by the ordinate range and constant 
beta curves in Figure 4-17. Accordingly one might conclude that a good 
indication of the system’s magnitude behavior has been conveyed by the 
PARAM-3, M0D0RD=1 graph for = 30 rad/sec. This is most emphatically 
not true.* Analysis of Figure 4-17 indicates that the greatest magnitude 
occurs at beta * 0.0 and that the magnitude contours decrease as beta ap- 
proaches 0.05. Careful review of Figure 4-15 will disclose a serious 
fallacy in such a conclusion. Figure 4-17 shows no indication of the 
magnitude peak resonance at beta = 0.004. Without this apriori knowledge, 
a completely erroneous interpretation would have been made by analyzing 
the curves in PARAM-3, M0D0RD=1. This example stresses once again the 
importance of viewing system frequency response as a function of the two 
system parameters. 

Consideration will now be given to the alternate PARAM-3 curve desig- 
nated by M0D0RD=2. This graph presents constant alpha curves in the beta 
magnitude plane for a fixed omega. As such, these curves simply slice 
vertical cross-sectional areas out of the magnitude contours of PARAM-2. 
Figure 4-18 shows such a set of curves for a constant angular frequency of 
30 rad/sec. Looking at Figure 4-15 again, all curve values lie within a 
fairly small range of alpha values. Accordingly, one would only expect 
graph curves for alpha less than 0.0021 to lie in the beta - magnitude 
plane for magnitudes greater than or equal to -3 decibels. Furthermore 
one would expect a series of magnitude peaks to occur at beta = .004 with 
the maximum amgnitude occuring on the alpha * 0.0005 curve. These factors 
are verified in Figure 4-18 where the alpha = 0.0005 curve has maximum 
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magnitude greatly in excess of 21 decibels. All other curves represent 
lower magnitude ranges as was anticipated. In this set of alpha curves, 
the maximum magnitude curve was intentionally constructed to indicate 
the large magnitude at peak resonance. This curve might have been over- 
looked, were it not for the apriori knowledge gained from Figure 4-15. 

The preceding discussion of PARAM-3 curves showed serious limita- 
tions in their usage when the PARAM-2 curve of the same constant omega is 
not utilized. When used in conjunction with PARAM-2 graphs, however, the 
PARAM-3 curves do provide some useful information. They are particularly 
useful for the area enclosed by the largest value of magnitude in the 
PARAM-2 plots. This area tells nothing about the peak resonant magnitude, 
other than the frequency at which it occurs. By constructing one PARAM-3 
curve along the bisector of these magnitude curves, one quickly gains the 
value of peak magnitude and the magnitude sensitivity about this resonant 
point. Should the precise peak magnitude location not be so obvious on 
other PARAM-2 plots, a series of curves in the PARAM-3 graph will clearly 
indicate its location. Due to the elliptical shape of the PARAM-2 magni- 
tude curves for this system, the PARAM-3, M0D0RD=2 configuration would 
be better suited than the PARAI'1-3, M0D0RD=1 graph. The reason for this 
is simply that the magnitude is more sensitive to alpha than beta. There- 
fore a much smaller range of alpha values need be considered than would be 
required for beta curves. 

Up to this point, only one value of omega has been discussed. The 
PARAM-3 curves for other angular frequencies can be visualized by studying 
Figures 4-13, 4-14, and 4-16. The PARAM-2 plot of Figure 4-13 shows a 
very gentle magnitude slope along constant alpha and constant beta lines. 
Thus for constant alpha lines, one would expect the magnitude to decrease 
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slowly with an increase in beta. The magnitude would decrease more 
rapidly for larger constant alpha curve values than they would for smaller 
curve values. Conversely, for constant beta lines, it is seen that magni- 
tude will increase slowly as alpha is increased. Large constant beta 
curve values would show a slower increase in magnitude than small beta 
curve values. Figures 4-19 and 4-20 are included for verification of 
these trends. 

Inspection of Figure 4-14 shows that the PARAM-3, M0D0R=1 graph will 
have a maximum magnitude at alpha = 0.0018 for all beta curves. The maxi- 
mum value of magnitude for any \alue of alpha will be obtained on the beta 
= 0 curve. All other beta curves will lie inside this magnitude boundary. 
The shpae of these curves will be similar to those shown on Figure 4-18. 
Since the magnitude is less sensitive at 20 rad/sec than at 30 rad/sec, the 
curves will be less pointed. The PARAM-3, MODORD-2 curves will not be of 
the shape just discussed. This is because the vertical cross-sections of 
the PARAM-2 plot show a decrease in magnitude as beta increases. All alpha 
curves will lie inside the constant alpha curve equal to 0.0018, for this 
curve represents a ridge of maximum magnitude. The constant 0-) = 40 rad/ 
sec magnitude curves of Figure 4-16 show that the peak resonance occurs on 
the beta axis at an alpha value of approximately 0.012. The sensitivity 
of magnitude is greatest in this graph. Accordingly, the PARAM-3, M0D0RD=1 
curves will be much more pointed than those shown in Figure 4-17. Also 
the center of these magnitude peaks will be shifted down to the abcissa. 

The PARAM-3, M0D0RD=2 curves will be likewise sharper than those shown in 
Figure 4-18, but this time the centerline of the magnitude peaks will lie 
on the beta = 0.012 line. 
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It has been shown that for this system the PARAM-3 curves are of 
limited use. Detailed correlation between PARAM-2 and PARAM-3 plots was 
made to show the inter-relationships which exist. The reason for discus- 
sing PARAM-3 graphs was to show how such information can be picked off of 
the PARAM-2 plots, thereby eliminating the requirement for another graphi- 
cal presentation. A more important reason for this detailed analysis was 
to point out the correlations applicable to this system under study. Many 
systems will undoubtedly be more complex than this. Such cases might re- 
quire the use of these curves. Therefore it is essential to obtain a firm 
understanding of these graphical advantages and limitations. 



PARAM-4. 

The graphical concepts of PARAM-4 will now be discussed for the third 
order system under consideration. Again, two basically different plots are 
obtained for this program. The PAElAM-4, M0D0RD=1 configuration graphs 
constant magnitude curves onto the alpha - omega plane for a fixed value 
of beta. Figure 4-21 shows one such plot for a given fixed beta value of 
0.036. At first glance, these curves seem very strange. Once again, how- 
ever, recourse to the two parameter frequency response curves of PARAM-1 
will show a striking relationship. Recall for a moment the three dimen- 
sional concept that was introduced in the study of PARAM-1 graphs. Each 
three dimensional curved elliptical cone shaped surface was a function of 
alpha, beta, and omega for a given value of magnitude. As magnitude was 
decreased, the conical surface base was observed to rotate into the obser- 
vation line of sight, while the overall size expanded. The simultaneous 
rotation and expansion caused the base portion to expand more than the 
smaller end. Furthermore, recall that omega increased in value from the 
base toward the small end. The PARAM-4, M0D0RD=1 graph effectively slices 
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a horizontal cross section out of this PARAM-1 three dimensional sur- 
face at the prescribed value of beta. To see this more clearly, Figures 
4-8 through 4-12 will be reviewed in order. Consider first the constant 
decibel magnitude PARAM-1 graph of Figure 4-8. By cutting a perpendicu- 
lar cross section along the beta = 0.036 line, it is seen that only one 
omega curve is severed namely the CO = 65 rad/sec curve. Note, however, 
that omega is a continuous elliptical surface; only a small number of 
cross sections are depicted in Figure 4-8. By interpolation, it can be 
seen that a perpendicular cross section of PARAM-1 will intersect the 
range of omega values between approximately 62 through 67. Furthermore 
note that no other omega curves have been severed in this manner. This 
intersection occurs at a range of alpha values between -0.002 and -0.004. 
Now turn to Figure 4-21 and inspect the 12 decibel magnitude curves re- 
produced there. By a similar analysis of Figures 4-9 through 4-12, it is 
possible to establish a positive correlation between the PARAM-1 and PARAM- 
4, M0D0RD=1 graphs. Next consider the PARAM-4, M0D0RD=2 plot shown in 
Figure 4-22. This configuration graphs constant magnitude curves onto 
the beta - omega plane for a fixed value of alpha. By methods analogous 
to those presented above, it is possible to slice a vertical cross section 
of the PARAM-1 graph at the prescribed value of alpha. The information 
thus obtained is seen to appear as PARAM-4, M0D0RD=2. The PARAM-4 curves 
are effectively topographical maps with magnitude equivalent to elevation. 
Their utility lies in the ability to correlate magnitude, frequency and one 
system parameter. They are restricted by the second system parameter which 
has been held at a constant value. As such they would only be suited for 
supplying additional frequency response information after an initial choice 
of one variable has been made. Providing a sufficient number of magnitude 
curves have been constructed on these graphs, they become in essence a 
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locus of Bode plots. When used in this manner, they become valuable 
tools. The constant alpha and beta values chosen for figures 4-21 and 
4-22 were so chosen because they coincided with the values used in 
figures 4-4 through 4-7. Thus a comparison of these PARAM-4 curves with 
the familiar Bode magnitude plots can easily be made. The PARAM-4 
curves can therefore be regarded as most useful for coarse Bode plot 
studies, since they portray magnitude-frequency information for a much 
larger number of alpha or beta values than would be feasible in a con- 
ventional Bode plot study. In this manner a design engineer would most 
likely construct enough PARAM-1 and PARAM-2 plots to make a judicuous 
first approximation for each of the system parameters. Using these, 
he could then construct both PARAM-4 plots and check his system perfor- 
mance with a fair degree of accuracy. After suitable corrections, he 
would then use PARAM-7 for a final refined system frequency response 
performance check. 

PARAM-5 

The discussions of PARAM-3 and PARAM-4 have considered cross section- 
al cuts of the PARAM-2 and PARAM-1 planes respectively. A slightly dif- 
ferent concept is portrayed in PARAM-5 and PARAM-6. Consider first the 
PARAM-5, M0D0RD=1 graph shown in figure 4-23. This shows a series of 
omega curves plotted upon the alpha magnitude plane for a fixed value 
of beta = 0.025. Each constant magnitude depicts a horizontal cross 
section of the PARAM-1 graph corresponding to that constant magnitude. 

The location of the horiztonal cross-section is fixed by the constant 
beta = 0.025 for the PARAM-5 plot. This is most clearly shown by the 
following comparison. Arbitrarily choose the constant magnitude line of 
6 decibels. This automatically selectes the PARAM-1 plot shown in 
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figure 4-9. Since beta is fixed at a value of 0.025, this is the 
horizontal line of figure 4-9 to be considered. Positive correlation 
between the omega curves crossing this beta =* 0.025 line in figure 4.9 
is thereby made with those omega curves crossing the magnitude line of 
6 decibels in figure 4-23. Effectively the PARAM-5, MODORD graph shows 
the series of PARAM-1 plots for any given constant beta value. Merely 
studying the different magnitude lines, one can see how a constant beta 
cross section would look on different PARAM-1 plots. The PARAM-5, 
M0D0RD=2 graphs show this same behavior, with the roles of alpha and 
beta interchanged. 

A canparison can also be made with a series of PARAM-2 plots, since 
they represent another three dimensional surface projected into the 
alpha beta parameter plane. In that case, however, each omega curve in 
the PARAM-5 graph represents a different PARAM-2 plot. Consider for a 
moment the ui=20 rad/sec curve of figure 4-23. This curve depicts the 
topological relief of a cross section of the PARAM-2 magnitude contour 
graph when that cross section occurs at beta = 0.025. Similarly, dif- 
ferent omega curves of figure 4-23 represent such a cross section in 
other PARAM-2 plots. 

The correlation between PARAM-5 and PARAM-1 as well as PARAM-2 pro- 
vides a rather important conclusion. One normally wants to get a good 
indication of system behavior by simultaneously varying alpha and beta. 
It is not necessary to construct a large number of PARAM-1 plots to do 
this. Neither is it necessary to construct a large number of PARAM-2 
plots. By constructing only a few PARAM-1 graphs, as well as two or 
at most three PARAM-2 graphs, one can gain valuable insight into the 
system behavior in the voids by studying the PARAM-5 plots. Thus for an 
initial system design study, one might be inclined to use the following 
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six graphs: two PARAM-1 plots, two PARAM-2 plots, and two PARAM-5 

plots. The constant magnitude values of the PARAM-1 plots would most 
likely be chosen from the bandwidth and resonant peak criteria. The 
constant omega values for the two FARAM-2 plots would normally be dictat- 
ed by the cut-off frequencies. Alternatly, the peak resonant frequency 
might be used as a governing choice. By choosing an intermediate value 
of alpha and beta, the PARAM-5 curves would indicate how other PARAM-1 
and PARAM-2 graphs would have appeared. For illustrative purposes, 
suppose a problem called for a -3 decibel bandwidth of 60 rad/sec, with a 
peak resonance magnitude of no more than 6 decibels around 40 rad/sec. 
Construction of figures 4-9, 4-12, and 4-16, one would be able to con- 
clude the following information; 

1) The tangent locus of the omega curves in figure 4-9 define 
the limiting values of alpha and beta as dictated by the 
resonance magnitude constraint. 

2) The intersection of the tangent locus and the constant 
UJ - ^0 rad/sec curve will determine which values of 
alpha and beta will give this desired resonant frequency. 

3) Looking at figure 4-12, one sees that the alpha and beta 
values determined above will not provide a -3 decibel band- 
width of 60 rad/sec. Instead, it can be seen by interpola- 
tion that the bandwidth will be approximately 57 rad/sec. 

4) Figure 4-16 shows the magnitude relationship at the de- 
sired resonant frequency. Thus by studying this graph, 
one sees that a resonant magnitude of approximately 7 
decibels will result from a 60 rad/sec bandwidth. 

The design example above has shown a coramon problem which exists for 
most engineering studies. The desired performance measures cannot be 
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attained; instead some compromise must be made. This conclusion was 
derived by simply observing three frequency response parameter plane 
graphs. The advantage of PARAM-5 graphs lies in the ability to corre- 
late performance sensitivities with parameter variations. Even though 
the curves shown in figure 4-23 do not coincide with the beta most likely 
to be selected, they do illustrate the behavior of magnitude and omega 
at a neighboring beta value. Thus it is possible to project that in- 
formation onto one of the PARAM-1, PARAM-2 graphs and extrapolate the 
data thus obtained. Additionally, one can do the same with the PARAM-5 
curves shown in figure 4-24. This figure shows the magnitude and omega 
variations as a function of beta for a constant alpha also near the 
point of interest. In effect the curves in figure 4-24 reflect the be- 
havior about the peak resonant point of l^>=35 rad/sec. It is seen that 
a variation in beta will increase the resonant magnitude for smaller beta 
values, with a corresponding decrease in magnitude for higher beta values. 
Figure 4-24 shows another important factor. At beta = 0.036 one sees that 
a flat frequency response of 0 decibels will be obtained for a bandwidth 
of approximately 35 rad/sec. This can be verified in figures 4-11 and 
4-25. 

Figures 4-25 and 4-26 have been included to show more vividly how the 
interpolation of PARAM-5 plots might be made. Consider the PARAM-5 graphs 
as topological portrayals of omega. Figure 4-25 shows how these omega 
contours change from those of Figure 4-23 as beta increases. Note that 
all omega curves are reduced in magnitude when beta is increased. Note 
further that a reduction in the number of omega curves will result in a 
loss of important information. The eight curves shown in figure 4-25 
might lead the observer to conclude that magnitude decreases continually 
as alpha is decreased. Only by careful inspection of the = 40 rad/sec 
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curve can one see that this is not the case. By inspecting the complete 
array of curves in figure 4-23, however, one sees a resonance trend for 
negative alpha values. In this figure, resonance occurs at approximately 
55 rad/sec. Using this information, another PARAM-2 graph can be visual- 
ized for constant = 55 rad/sec. The curves would be similar in shape 
to those sho^Tn in figure 4-15 and 4-16, but smaller in size. The center 
of these elliptical curves would occur at beta = 0,036 and alpha ~ -0.003. 
The analogous procedure could also be employed using either figures 4-24 
or 4-26. 

PARAM-6 

The discussions of PARAM-3 through PARAM-5 have shown how various 
graphical configurations result from taking cross sectional cuts in 
the PARAM-1 and PARAM-2 frequency response parameter plane curves. The 
curves shown in PARAM-6 depict another such cross sectional projection. 
Consider first the curves shown in figure 4-27. They correspond to 
constant beta curves mapped onto the alpha-omega plane for a given pre- 
selected magnitude value of 12 decibels. Accordingly, this graph corres- 
ponds to the PARAM-1 plot shown in figure 4-8. The full utilization of 
the entire set of graphs will now become clear. Effectively the PARAM-6 
graphs show the third variable of the curved elliptical surface which 
was projected onto the parameter plane in PARAM-1. In this manner, the 
third variable, omega, now becomes one of the axis variables while cross 
sectional cuts of the three dimensional alpha-beta-omega surface are ob- 
tained by slicing down along constant beta lines or constant alpha lines. 
The latter cross sections are pictured in figure 4-28. Thus when the 
corresponding PARAM-1 and PARAM-6 graphs are viewed simultaneously, one 
can visualize the three dimensional alpha-beta-omega surface in its true 
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form. The illustrative example of an elliptically shape Dull's horn is 
applicable in all but one facet which will now be corrected. It was 
originally stipulated that the omega axis passed down the conical surface, 
and was rotated and shifted slightly as magnitude decreased. This stipula- 
tion was made to simplify the visualization of the surface. Inspection of 
the PARAM-6 curves will show a slight modification of this illustration. 

The omega axis does not pass down the surface, but instead projects per- 
pendicularly through the origin of the PARAM-1 plots. In this manner it 
is seen that the surface has shifted slightly, and been rotated in the 
three dimensional alpha-beta-omega space.' Now the illustration is cor- 
rect in its entirety. Analyzing the corresponding PARAM-1 and PARAM-6 
graphs gives a complete description of the portrayed surface in three 
dimensions. Inspection of any one of the PARAM-2 plots reveals the 
change in shape of this three dimensional surface as the fourth dimen- 
sional parameter, magnitude, is varied. Figures 4-29 through 4-32 are 
included for amplification and verification of the above results. 
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5. 



CONCLUSIONS 



Section 4 presented an exhaustive study of the application of fre- 
quency response parameter plane curves to one particular system. The 
selection of the third order system was based on the ease of interpre- 
tation of each frequency response graph discussed. Extensions to systems 
of higher order will undoubtedly change the graphical configurations to 
some extent. The degree of this variation was smaller than expected for 
a similar fourth order system. Figure 5-1 shows the particular fourth 
order tachometer and acceleration feedback system which was studied. 

Figure 5-2 presents the system stability performance curves in the alpha- 
beta parameter plane. Figures 5-3 through 5-12 depict the associated 
frequency response curves. Comparison with the third order system discuss- 
ed in section 4 shows strikingly similar results. 

Several concluding comments are noteworthy concerning the use PARAM-1 
through PARAM-6. The alpha - beta parameter plane provides more informa- 
tion than any other representation for system studies involving the simul- 
taneous variation of two parameters. Several different PARAM-1 plots 
are considered to be most useful for initial system studies, when they 
are accompanied by their corresponding PARAM-6 graphs. In this manner, 
a three dimensional space of alpha, beta, a’^d omega can be visualized for 
each of the constant magnitudes used. The additional incorporation of one 
or two PARAM-2 graphs will describe the effect of a variation in the fourth 
dimension. Stated in other words, the information in PARAM-2 shows how the 
three dimensional PARAM-1, PARAM-6 surface varies in magnitude. PARAM-3 
and PARAM-5 curves yield a more sensitive portrayal of the magnitude varia- 
tions; however these are restricted to a single parameter variation. They 
are useful for accurately presenting a maximum peak resonance magnitude 
when this information is not available from the peak resonant tangent 
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locus of PARAM-1. PARAM-4 effectively shows a series of Bode curves 
as a function of one system parameter. As such they would be advantag- 
eous for intermediate studies, when the final adjustments are being 
contemplated. Finally, the familiar Bode plot representations of PARAM-7 
show the most accurate frequency response, and are therefore most useful 
in the terminal stages of design. 
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i<-SCRLE=1.00E-03 UNITS INCH. 

V-3CRLE=1.00E-Q2 UNITS INCH. 

GLRUIS, PRRRMETER PLRNE CURUES, 4TH ORDER T/R FB 
S**^+160Ss^*3+ (6. 5E3+5E6R) 3**2+ [5Ea+5E6B) S+5E6=0 

FIGURE 5-2 
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>^-SCRLE=1.00E-03 UNITS INCH. 
y-SCRLE=5.00E-03 UNITS INCH. 

GLRUIS. PflRflM 1, 4TH ORDER T/R F.B.,OMEGR CURUES 
ORDINRTE=BETR, RBCISSR=RLPHR, CONMRGNI TUDE= 6DB 



FIGURE 5-3 
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X-SCRLE=1.00E-03 UNITS INCH. 
y-SCRLE=5.00E-03 UNITS INCH. 



GLRUIS. PflRRM 2, ^TH ORDER T/R F-B-, MRG CURUES 
ORDINRTE=BETR. RBCISSR=RLPHR, C0N0MEGR=30DB 



FIGURE 5-4 
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V5CFiLE=3.00E+00 UNITS INCH. 
^-3CHLE=1.00E-03 UNITS INCH. 



TlRUIS, PflRRM 3, 4TH ORDER T/fl F.B.. BETR CURUES 
""RDINRTE=RLRHR, RBCISSR=MRG , C0N0MEGR=30RRD/SEC 



FIGURE 5-5 
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X-SCRLE=3.00E+00 UNITS INCH. 

V-SCRLE=5.00E-03 UNITS INCH. 

GLRUIS, PRRRM 3, ^TH ORDER T/fi F-B.,RLPHR CURUES 
ORDINRTE-BETfi , RBCISSfi=MRG . C0N0MEGfi=30RRD/SEC 



FIGURE 5-6 
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s 




/-„'nLE=:.00E--01 UNITS INCH. 

•'-3ChLE=1.00E-03 units INCH. 

F.B. , MflG CIIRUEb 
flBCISSn=OMEGR, C0NB£Tfl=0. 0^6 



PRRRM dTH ORDER T/R 



'REINRTE=RLPHR, 



FIGURE 5-7 
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X-SCRLE=1.00E-*-Ql UNITS INCH. 



V-SCRLE=5.00E-03 UNITS INCH. 

GLRUIS, PRRRM d, 4TH ORDER T/R F-B., MRG CURUES 
ORDINRTE=BETR , RBCISSR=OMEGR, G0NRLPHR=0. 0008? 



FIGURE 5-8 
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X-:."RLE=3.00E-^00 UNITS INCH. 

•'-:"W-E=1.00E-03 UNITS INCH. 

oi.ril. IS, PnRRM 5, 4TH ORDER T/R F.B.,OMEGR CURUES 
DRDIRRTE^RLPHR, RBCISSR=MRG , G0NBETR=0. 0^6 



FIGURE 5-9 
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X-SCRLE=3.00E+00 UNITS INCH. 

V-SCRLE=5.00E-03 UNITS INCH. 

GLRUIS, PRRRM 5, dTH ORDER T/R F.B.,OMEGR CURUES 
ORDINRTE=BETR , RBCISSR=MRG . C0NRLRHR=0. 00087 

FIGURE 5-10 
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V-SCflLE=1.00E-01 UNITS INCH. 

■'-SCRLE=1.00E-03 UNITS INCH. 

GLRUIS, PRRRM 6, 4TH ORDER T/R F-B-, BETR CURUE 
ORDIHRTE^Beffl, RBCISSR=OMEGR, CONORGNI TUDE=6DB 



FIGURE 5-11 
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X-SCRLE=1.00E+01 UNITS INCH. 

V-SCfiLE=1.00E-03 UNITS INCH. 

GLRUIS. PRRRM 6, 4TH ORDER T/R F.B.,RLPHR CURUES 
0RDINRTE=RLPHR, RBCISSR=OMEGR. G0NMRGNITUDE=6DB 

FICUSE 5-12 
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APPENDIX I 

PROGRAM PARAMS ♦♦♦*♦♦♦♦♦* 

C THIS PROGRAM CONVERTS A GIVEN TRANSFER FUNCTION TO A SET OF FREQUENCY 
CDOMAIN REPRESENTATIONS FOR ANALYSIS AND DESIGN STUDIES. THIS TRANSFER 
CFUNCTIUN IS READ AS A RATIO OF TWO COMPLEX FREQUENCY POLYNOMIALS. WITH 
CTWO VARIABLE PARAMETERS ALPHA AND BETA EMBEDDED IN THE COEFFICIENTS. 
CTHE PROGRAM THEN TRANSFORMS THIS TRANSFER FUNCTION INTO THE ANGULAR 
CFREOUENCY DOMAIN, YIELDING AN EQUATION IN TERMS OF FOUR VARIABLES: 
CALPHA, BETA, MAGNITUDE, AND OMEGA. SEVEN TWO-DIMENSIONAL GRAPHICAL 
CREPRESENTATIONS ARE GENERATED BY FIXING ONE VARIABLE AS A CONSTANT; 
CASSIGNING DESIRED CURVE VALUES TO ANOTHER; AND LETTING THE REMAINING 
CTWO VARIABLES BE THE ORDINATE AND ABCISSA VALUES. A THIRD COEFFICIENT 
CVARIABLE PARAMETER, GAMMA, MAY ALSO BE EMBEDDED INTO THE POLYNOMIAL. 

CIN THIS CASE, THE ABOVE PROCESS IS REPEATED FOR EACH VAIUF OF THF. THIRD 

CVARIABLE. DATA CARD INFORMATION FOLLOWS 

card information SECTI 0N************+********** 

CCARD FORMAT SYI 



8110 



NORN 

NORD 

MODORD 



NALPHA 

N8ETA 

NMAG 

NOMEGA 

NGAMMA 



8110 IMAG 



LX 

LY 

CONALP 

CONBET 



USAGE OF SYMBOL 
NUMBER OF DESIRED RUNS. EACH MUST HAVE A COMPLETE 
SET OF DATA CARDS LISTED BELOW. NOT EVERY SYMBOL 
IS USED IN ALL CASES. LEAVE UNUSED SYMBOLS BLANK, 
on NOT OMIT ANY BLANK CARDSj UNLESS SO SPECIFIED. 
INSERT ADDITIONAL CARDS BEHIND EACH FULL CARO. 
NUMBER OF GRAPHICAL PRESENTATIONS TO BE PROGRAMMED 
NUMERICAL TITLES (1 TO 7) OF EACH PARAM TO BE PUN; 
MUST AGREE IN NUMBER WITH KPARAM ABOVE. 

ORDER OF NUMERATOR POLYNOMIAL (LIMITED TO 49) 

ORDER OF DENOMINATOR POLYNOMIAL (LIMITED TO 49) 
Y-AXIS VARIABLE, AND CURVE VARIABLE CONTROL 
MODORD=l: ALPHA IS ORDINATE VARIABLE (PARAM l-6) 
M0D0RD=2: BETA IS ORDINATE VARIABLE (PARAM 1-6) 
MCDORDsl: ALPHA CURVES WILL BE PLOTTED (PARAM 7) 
M0DURD=2; BETA CURVES WILL BE PLOTTED (PARAM 7) 
NUMBER OF ALPHA CURVES DESIRED (LIMITED TO 16) 
NUMBER OF BETA CURVES DESIRED (LIMITED TO 16) 
NUMBER OF MAGNITUDE CURVES DESIRED (LIMITED TP 16) 
NUMBER OF omega CURVES DESIRED (LIMITED TO 16) 
POSITIVE INTEGER DENOTES NUMBER OF THIRD PARAMETER 
VALUES TO BE EMBEDDED INTO EACH PARAM RijN. INSERT 
COEFFICIENTS FOR THESE GAMMA VALUES ON SEPARATE 
CAROS DIRECTLY BEHIND REST OF COEFFICIENTS; SEF«19 
NEGATIVE INTEGER DENOTES NUMBER OF PAPAM-7 VECTOR 
ADDITIONS FOR GENERALIZED BLOCK DIAGRAM STUDIES. 
THE (NEGATIVE) VALUE OF NGAMMA DETERMINES HOW MANY 
BODE CALCULATIONS WILL BE MADE BEFORE PLOTTING THF 
VECTOR RESULTANT AS ONE MAGNITUDE, AND ONE PHASE. 
INSERT ADDITIONAL SETS OF COEFFICIENTS AT FNO . 
MAGNITUDE VALUE FORM CONTROL 

IMAG=0: ALL VALUES IN DECIBELS (PARAM 7 THIS FORM) 
IMAG=1 MEANS ALL VALUES IN ACTUAL MAGNITUDE 
IMAG=2 MEANS ALL VALUES IN MAGNITUDE SQUARED 
PRINTOUT CONTROL 
IWRITE=-1 means no printout 

IWRITE=0 means only INPUT DATA WILL BE PRINTED 
IWRITE=1 TO 900 YIELDS INPUT DATA, GRAPHICAL 

information, and computational results, the number 

OF COMPUTATIONS BETWEEN EACH PRINTOUT IS DICTATED 
BY THE VALUE OF IWRITE. IF IWRITE IS GREATER THAN 
900, ONLY INPUT DATA G GRAPH INFO. WILL BF PLOTTED 
FREQUENCY SCALING CONTROL 
IGRAPH=-2: NO GRAPHS WILL BE PLOTTED. 

IGPAPH=-1 :ONLY PARAM 7 WILL BE PLOTTED ON LOG GRID 
IGRAPH=0; PARAMS 6£7 WILL BE PLOTTED ON LOG GRID. 
IGRAPH=1 :PAPAMS 4,6, G7 WILL BE PLOTTED ON LOG GRID 
IGRAPH=2: PARAMS 4G7 WILL PE PLOTTED ON LOG GRID. 
PARAM 7 PHASE PLOT CONTROL 
IPHASE=-1 MEANS NO PHASE CURVES 

IPHASE^O: phase CURVES SUPERIMPOSED ON MAGNITUDE 
IPHASE=1: PHASE CURVES PLOTTED ON SEPARATE GRAPH 
PARAM 7 PHASE PLOT SCALING CONTROL 

I0UAD=0 YIELDS PHASE CURVES FROM 0 TO -180 DEGREES 

I0UAD=1 YIELDS PHASE CURVES FROM-90 TO 4-90 DEGREES 

PARAM 7 POLAR PLOT CONTROL 

IP0LAR=0 MEANS NO POLAR PLOT 

IPOLAR=l GIVES POLAR PLOT ON SEPARATE GRAPH 

LENGTH. OF X-AXIS (LIMITED TO 9 INCHES) 

length of Y-axis (limited to 15 inches) 

VALUE OF ALPHA WHEN USED AS A CONSTANT 
VALUE OF BETA WHEN USED AS A CONSTANT 
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c 


CUNMag 


MAGNI 


c 


CUNOMt 


VAIUE 




•♦*DATA 


CARD r 


CCARO FORMAT 


Symbol 


( 


C 6 8F10.3 


ALPHmn 


MIN. ' 


C 


ALPHMX 


MAX. ' 


C 


BFTAMN 


•MIN. ' 


C 


RE TAMX 


MAX. ’ 


C 


XMAGMN 


MIN. ' 


C 


XMAGMX 


MAX. ' 


C 


OMEGMN 


MIN. ( 


C 


CMEGMX 


MA X. 1 


C 




NIUE : 



A CONSTANT; SEE I^AG 
A CONSTANT (NOT ZERO) 



USAGE OF SYMBOL 



8E10. 3 
Pt 10. 1 
8E10.3 
8E 10. 3 
3E 10. 3 
8E 10. 3 
6E10. 3 
8E10.3 
8E10.3 
HEIO. 3 
PE10.3 
8E10. 3 



WHEN USED AS an axis VARIABLE 
VALUE DP ALPHA WHEN USEO AS AN AXIS VARIABLE 
VALUE HE BETA WHEN USEO AS AN AXIS VARIABLE 
MAX. VALUE OE BETA WHEN USEO AS AN AXIS VARIABLE 
MIN. mag. value when used AS AXIS VAR.*, SEE IMAG 
MAX. MAG. VALUE WHEN USED AS AXIS VAR.; SEE IMAG 

omega value when used as axis variable 
omega value when used as axis vaoiarle 

_!IE log GRin IS REQUESTED, ( SEE IGRAPHI, OMEGMN 
will auto, changed to nearest power OE 10, AND 
PPEGMX CHANGED LIKEWISE. THE NUM8EF OF DECADES IS 
automatically COMPUTED FROM THESE NEW VALUES. 
DESIRED VALUES FOR ALPHA CURVES (SEE NALPHA) 
DESIRED VALUES FOR BETA CURVES (SE^ NRETA) 

DESIRED VALUES «=0R MAGNITUDE CURVES! SEE MMAGEIMAG) 
DESIRED POSITIVE OMEGA CURVE VALUES (SEE MOMFGA) 
NIImpraTOP CONSTANT COEFFICIENTS IN ASCENDING ORDER 
NUMERATOR ALPHA COEFFICIENTS IN ASCENDING ORDER 
NUMERATOR BETA COEFFICIENTS IN ASCENDING ORDER 
NUM. alpha-beta coefficients IN ASCENDING ORDER 
denominator CONST. COEFFICIENTS IN ASCENDING ORDER 
denominator alpha COEFFICIENTS IN ASCENDING ORDER 
denominator beta coefficients in ascending ORDER 
DtNOM. alpha-beta COEFFICIENTS IN ASCENDING ORDER 
OMIT CARDS * IR THRU « 27 IF NGAmma LESS THAN ONE. 
VALUES OF THE THIRD PARAMETER; MUST MATCH NGAMMA 
NUP. gamma-const. coefficients IN ASCENDING ORDER 

NUM. oAMma-ALPHA coefficients in ASCENDING ORDER 

gamma-beta coefficients in ASCENDING ORDER 

GAmma-ALPHA-BET A COEFFS. IN ASCENDING ORDER 

GAMMA-C0N»^T. COEFFICIENTS IN ASCENDING ORDER 
gamma-alpha COEFFICIENTS IM ASCENDING OPDFR 

gamma-beta coefficients in ascending opdfr 

_ . gamma-alpha-beta COEFFS. IN ASCENDING OPDFR 
FIRST LINE graph TITLE (COLUMNS ITHRijAB) 

SLCOND LINF OF GRAPH TITLE (COLUMNS ITHRIIAB) 
the number of titles fop »ARAM 7 MUST MATCH THF 

number of graphs requested. ( SFF IRHASF C IPOLAR) 
1ST 1/2 PHASE TITLE;0MIT UNLESS NPAR Am= 7R I PHASF=1 
2NC 1/2 PHASE TITLF;nMlT UNLESS NPAR A«= 7S I PH AS F= 1 
1ST 1/2 POLAR TITLF;OMIT UNLESS NP AR AM = 7r. I POL AR= 1 
2ND 1/2 POLAR TITLE;OMIT UNLESS NPARAM = 7P.I P3LAR=1 

■.t*t***************^******** ******** ********^^***************:***P«^^** 

REALmSIT ITLE( 12) , JTI TLE( 12) ,KTI TLE( 12) ,J0KFR(9)/1H , SHPAP AM- , 7HFijP 
ether, 5MALPHA,4HBETA .SHOMEGA, BHMAGN I TUO , 5HG A mm a , BHS EP AR ATE / 

PEAL DM/1M1/,LA( 66)/4H 0.1, 4H C.2,AH 0,3, AH O.A,AH 0.‘^,AH O.ft.AH 0 
E.7,AM 0.8 , ah 0,9, ah 1.C,AH 1.0, ah 2,0, AH 3.0, AM A.O,AH G.O,AH 6.0, 
CAM 7.0, ah B.O.AH 9.0,AH1C.O,AH10.0,AH20.0 ,AH30.0,AHAO.O, AH50.0, 
r. AM60.0, AH70.0, AM8C. 0, AM90.0,AHIOO. ,AH100.,AH200.,AH300. , AHAOO. , 
EAH500. ,AH600. ,AH700. , AHBOO . , AH900 . , AH 1000 , AH 1000 , AH2000 , AH3000 , 
EAHAOOO,AM5000,AH60CO,AH7000,AH8000,AH9000,1H ,1H1,1H7,1H3,1HA,1H5, 
£1H6, 1H7,hb,1h9,2H10,2H11,2H12,2H13,2H1A,?h16,7H16/ 

DIMENSION JPAHAM (E ) , A( 50) ,B( 50) ,D( 50) ,F( 60) ,HA( 50) ,HAT( 50 ) ,HBT( 50 ) 
f. ,HCT(50),HD(50) ,HDT( 50) ,HET( 50) ,HF T ( 50 ) , OL DMA ( 50 ) , OLDHA T ( 50 ) . OL DHB 
ET( 50) ,OLDMCT( 50) ,0L0HD(50) ,C1L0HDT( 50 ) , OL DHET ( 50 ) ,01 OHFTI 50 ) ,GAMhA( 
C50),GAMHAT(50) ,GAmhbT( 50 ) ,GAMHCT (50 ) ,GAMHD( 50) ,GAMHOT( 50 ) , GAMHFT( 5 
£0) ,GAMHFT( 50) , IC ( 18 ) ,Z ALPHA! 16),Z8ETA(18) ,Z GAMMA! 1 6 ) , 7 MAGSO ( 1 6 ) , 70 



0. 3 
0.3 
0. 3 
0. 3 
0.3 
0 , ‘ 



8E10. 3 

6AR 

6AB 



6A8 

6AB 

PA3 



Z ALPHA 

ZBETA 

ZmaG 

Z OMEGA 

HA 

HAT 

HBT 

HCT 

HD 

HOT 

HET 

HFT 

• * ♦ ♦ 

ZGAmma 

gamha 

gamhat 

GAMhrt 

GAmhCT 

Gamhd 

GAMhdt 

GAMhE T 

GAMMFT 

(TITLE 

(TITLE 



JT I TLE 
JT ITLE 
KTITLl 
RT ITLE 



NUM. 
NUM. 
DEN. 
DEN. 
DE^ . 
DEN. 



,900),W(16,900),X(900), 



£MEGA( 16) ,S( 16,900) ,T( 16,900) ,U( 16,900) ,V 
£Y( 900) 

CALL CANCcL(2) 

EQUIVALENCE! S,T,U,V,W) 

1 FORMAT(PIIO) 

2 FORMAT! 10X.8I lO) 

3 F0RmaT(6AB) 

A FORMAT! lOX ,6A8) 

5 FORMAT! BE10.3) 

6 FORMAT! lOX ,1P lOE 12. 3 ) 

***********************\t\tP[,jj data READ IN ************************** 
READ! 5,1) VPUM 
DO 1000 IRUN=1,NRUN 

PEAD!5,l) KP ARAM, ( JPARAM(I), 1 = 1, KP ARAM) 

READ! 5, 1 )NORN,NORD,MOOORn,NALPHA,NBETA,NMAG,NOMFGA,NGAMMA 
PEA0!5,1) IMAG , (WRITE , I GR A PH , I PH ASE , I QUAD , I POL AR , L X , L Y 
READ! 5,5 )CONALP,CP^'8ET,CONMAG,CONOMF 
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REAO( 5i 5) ALPHMN, ALPHHX,8ETAMN»BETAMX|XMAGMN»XMAGMX,0MEGMN, OMEGMX 
REAO(5,5) ( ZALPHA( I ) ,I*1,NALPHA) 

REAO(5,5) (ZBETA(I),t=l,NBETA) 

REAO(5«5) ( ZGAMMA( I i » iNGAMMA) 

READ(5,5) (ZMAGS0( I) ,1=1 ,N»1AG) 

REAO( 5, 5) (ZOMEGA( I ) , I =1 , NOMEGA ) 

LOE=NDRN+l 
MOE=NDRO+l 

REAO(5,5) (HA (I), 1=1, LOE) 

READ(5,5) (HAT(I), 1=1, LOE) 

REAO(5,5) (HBT(I), 1=1, LOE) 

READ(5,5) (HCT(l), 1=1, LOE) 

REAO(5,5) (HO ( I), 1=1, MOE) 

REAO(5,5) (HOT(I), I = l,*iOE) 

REAO(5,5) (HET(I), 1=1, MOE) 



READ(5,5) (HFT(I), 1=1, MOE) 

F( (NGAMMA ) ,GE.( 1) ) READ(5,5) (ZGAMMA( 



! ZGAMMA! 


) ,1 = 


, NGAMMA) 


(GAMHA ! 


) , 




,LOE) 


(GAMHAT! 


), 




,LOF) 


(GAMHBT! 


) , 




,LOE) 


(GAMHCT! 


) , 


X 


,LOE) 


(GAmhD !1 


) , 


a 


,MOE) 


(GAMHDT! 


) , 


= 


,MOE) 


(GAMHET! 


I) , 


-s 


,MOE) 


(GAMHFT! 


) , 


= 


,MOE) 



I- 1 l . I 1 J I 

IF( ( IWRITE).LT.(O) ) GO TO 2A 

NPUT data print out SFCT inN*****^*********** 

WRITE (6, 7) IRUN,NORN,NORD,KPARAM,( JOKER(2) ,JPARAM( I),I=1,KPARAM) 

7 FORMAT( IHl , lOX, • ♦♦♦♦♦ THE INPUT DATA FOP RUN*, I?,' IS ♦♦♦♦♦*,//, 
610X, ’ORDER OF NUMERA TOR : • . I 2 , / / , lOX . • ORDER OF DENOMI N ATOR : • , I 2 , / / , 
610X,*THE FOLLOWING', 12, • GRAPHICAL CONFIGURATIONS WILL RE COMPUTED 
G:*,//,20X,8(A8,I1,3X),/) 

IF( (NGAMMA) ,GE. ( 1 ) ) WRITE(6,8) NGAMMA, I ZGAMMA (I ), I = 1 , NGAMMA ) 

8 FQRMAT( lOX, 12, • GAMMA, OR THIRD VAR. PARAMETER COEFFICIENTS WILL P 
6E USED. • ,//,10X, 'THE GAMMA VALUES AR E : • , 1 P8 E 12 . 3 , / , 10 X, 1 P8E 1 2 . R , / ) 

WRITE(6, 10) 

10 FORMAT! iOX, 'THE COEFFICIENTS FOR THE NUMERATOR ARE:*,/) 

DO 11 1=1, LOE 

J=I-1 

IF(NGAMMA .LE.O) WRITE(6,12) HA (1 ) , H AT ( I ) , HRT ( 1 ) , HCT ( I ) , J 

12 F0RMAT(20X, • ( *,1PE10.3, * + *.lPE10.3, **ALPHA 4- • , 1 P E 10 . 3 , * *B ETA + * 

G, 1PE10.3, **ALPHA*8ETA ) ♦ S**( *, 12 ,*)*,/ ) 

IF(NGAMMA .GT.O) WRITE(6, 9) HA ( I ) , GAMHA ( I ) , HAT ( I ) ,G AMHAT ( I ) , HBT ( 1 
G) ,GAMHBT( I ) ,HCT( I ) , GAMHC T ( I ) , J 

9 FORMAT! 15X, M(*,1PE10.3,* + *,IPE10.3, **GAMMA) + ( * , 1 P E 10 . 3, * ♦ • . IPF 1 0 . 
G3, **GAMMA)* ALPHA .... * ,//,20X, * *■ ( * , 1 PE 10 . 3 , * ♦ * , 1 PEIO . 3 , **GAMMA)*BE 
GTA4-( *,1PE10.3,*4-*, IPEIO. 3, **GAMMA) ♦ALPHA*BETA ) *♦ ( * ,12, * ) * , / ) 

11 CONTINUE 

WRITE(6, 13) 

13 FORMAT! /lOX, * THE COEFFICIENTS FOR THE DENOMINATOR ARP:*,/) 

DO 1<, 1 = 1, MOE 

J=I-1 

IFINGAMMA .LE.O) WRITE(6,12) HO ( I ) ,HDT ( I ) , HET ( I ) , HFT ( I ) , J 

GT.O) WRITE!6, 9) HO ! I ) ,G AMHD ( I ) , HOT ( I ) , G AMHDT ( I ) , HET ( T 



MODORO 



NALPHA 



NRFTA 



IFINGAMMA . . . 

G) ,GAMHET! I ) ,HFT( I ) ,GAMHFT( I ) , J 
lA CONTINUE 
WRITE(6,15) 

15 F0RMAT!15X,' NORN NORD 

GG NOMEGA NGAMMA *,/) 

WRITE!6,2 ) NORN, NORD, MODORD, NALPHA, NBETA.NMAG, NOMEGA, NGAMMA 
WRITEI6, 16) 

16 FORMAT! /, 17X, * IMAG IWRITE IGRAPH IPHASE IQUAD 

GLAR LXAxtS LVAXIS »,/) 

WRITE(6,2) IMAG , I WR I TE , I GR APH , I PH ASE , I QUAD, I POL AR , LX , L Y 
WRITE! 6, 17) 

17 FORMAT! / ,12X, * CCNALPHA CONBETA CONMAG CONOMEGA 

WRITE (6, 6 )CONAL P,CONRET, CONMAG, CONOME 

WRITEI6, 18) 

18 FORMAT! /,12X, • ALPHMIN ALPHMAX BETAMIN PETAMAX 

GMIN XMAGMAX 0«EGAMIN OMEGAMAX *,/) 



, /) 



XMAG 



TE! 6, 6) ALPHMN ,ALPHMX, BET A MN, BET AM x, XMAGMN , XMAGMX , OMEGMN , OMEGMX 
TE(6,19) JOKER(A) ,!ZALPHA!I ) ,1 = 1, NALPHA) 

WR TE(6,19) JOKER! 5) , ( ZBETA ! I ) , 1=1 , NBETA ) 

WR TE(6,19) J0KEP!7),(ZMAGS0(1),I=1,NMAG) 

WR TE(6,19) J0KER!6),(Z0MEGA(I ), 1=1, NOMEGA) 

19 FORMAT! /,10X,A8, 'CURVE V ALUE S: * , 1 P 8 E 12 . 3 , /3 IX , 1 P8 E 1 2 . 3 ) 

IF! IWRITE) 20,20,22 

20 WRITE!6,21) 

21 FORMAT! /,10X, 'THERE WILL BE NO PRINTOUT OF COMPUTED RESULTS') 
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GO TO 2k 

22 w«ITE(6,23)IW«I TF 

23 FURMAT( /, lOX, 'CJMptjTcO RESULTS WILL HE PRINTED EVERY', 13,* STEPS*) 

c**#**^***o***»****initializing computational set tings **♦♦*♦ *<'**♦*♦♦♦♦** 

2k IF( MOUORO.GF.2) GO TO 25 
oummyi=«etawn 

nUMMY2*BETAMX 

betamn=alphmn 

betamx=alphmx 

alphmn=pum«yi 

ALPH«X=OUM«Y2 

25 DD 1000 IPARAM=l .KPARAM 
NPARAm=jPARAM(IPARAM) 

mgamma=ngamma+i 

00 1000 IGAwma=1 ,mgAm'^A 
GAMMA = ZGAMMA( IGA**MA) 

IF(NGAMma.GE.I) Gr TO 517 

26 XS= ( ALPHMX-ALPHMN)/LX 
IF((ABSIXS)).LE.(1.0E-I0n XS=l.OE-lO 
YS=(BETAMX-rtfcTAMN)/LY 
IFnABSfYSn.LE.I l.OE-10)) YS=1.0F-10 

1 x = -betamn/ys 

I Y=-ALPH“N/XS 
M0D=1 
JCURVE=C 
IPAR7=1 

RCAD(5,3)( ITITLEI I ),I=1,12) 

IE(NPARAm.nE.7)G0 to 3A 
IM I“AG-l)27,2P,2‘i 

27 DPwax = xmaG‘^X 
CBM IN=X'^AG'«N 
GO TO ?C 

28 CRMAX=20*ALnG10 ( XMAGmx) 

CrtM in = 20*ALOG10( XMAG'-N) 

GO TO 30 

24 CBMAX=10*ALnGl0( XMAGMX) 

0BMIN=IC*AL0G1C( XMAG^N) 

30 XS=(Ommax-OBMIN)/LX 
I Y=OBMAX/XS 

or 31 1=1,16 
S( I , I )=0.0 
T( I ,1)=0.0 
UI I , 1)=0.0 

v( 1 , 1 )=o.o 

31 WI I ,l)=C.O 

IF( IPHASr.EO.O) GO TO 32 
READ! 5,3)( JTITLFI 1 ), 1 = 1,12) 

32 IF( ( IPOL AR ) .LE . ( 0) ) GJ TO 33 
READ(5,3)(KTITLF(I),I=1,12) 

33 IF( MQOORO-2) A5,AA,4A 

^k GO TO (35,A0,A1,46,A7,51 ,5R,59) ,NPARAM 

35 NL I MlT=NOMbGA 

36 IF( I«AG-l)37,3ri,3S 

37 OUMMYI=10.0**(CONmag/10.0) 

IE(NPARAM-l )59,59,A3 

38 OUMmyI=CONMAG*CONMAG 
IF( NPARAM- 1)59,59,^3 

39 OUMMYI=CCNmag 
IF(NPARAM-1)59,59,43 

^0 NLImIT=NMAG 
k\ OMEGA=CrNOME 

IF(NPARAm-2)59,59,4? 
k2 XS= ( XMAGMX-XMAGMN) /LX 
I Y=-XMAGMN/XS 
h 3 IF( M000R0-2)'tA,45,A5 
kk NLIm1T=N9ETA 

IF(nPARAM-6)59,53,55 

45 NLIMIT = N«LPHA 
IF(NPARAM-6)59,53,55 

46 NL1“IT=NMAG 
XS=(OMEGMX-OMFGMN) /LX 
I Y=-CMEGMN/XS 

47 IF(Ma0ORO-2)48,49,49 

48 CONST = CON!3ET 

I F( NPARAM-4) 52, 52, 50 

49 CONST=CnNALP 
IF(NPA9 Am-a)52, 52,50 

50 XS= ( XMAGMx-Xma&mn) /LX 

I y=-xmagmn/xs 
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NL IMIT=NOMEGA 
GO TO 5<3 

51 XS=(OMEGMX-OMEGMN»/LX 
IY=-OMEGMN/XS 

GO TO 36 

52 IF( IGRAPH)59,59.5A ^ 

53 IF( ( IGRAPH.lt. 01 .OR. (IGRAPH.GT.D) GO TO 59 
5A IY=LX-IX 

55 IF( IGRAPH.LE.-2) GO TO 59 

GRID GRAPHICAL CONST RUC T I OM* ♦*♦♦*♦♦****♦♦*♦*♦ *♦* 
IF (OMEGMN.lt. 1000. )YL0GMN=100. 

IF(OMEGMN.LT.lOO. )YL0GMN=10. 

IF(OMEGMN.LT. 10. )YL0GMN=1. 

1F(0MEGMN.LT.1. )YL0GMN=0.1 
1F(0MEGMX.GT.1.0 )YL0GMX=10.0 
IF( OMEGMX.GT. 10.0) YLOGMX =100.0 
IF( 0MEGMX.GT.100. ) YLOGMX = 1000.0 
I F( OMEGMX.GT. 1000. ) YLOGM X= 1 0000 . 0 
0UMMY3= YLOGMX/YLCGMN 
NUMDEC=ALOG10(DUMMY3) 

IX=0 

DlJMMYA=NUMOEC 
nUMMY5=LY 
YS=DUMMYA/DUMMY5 
IFINAL=LX*1 
JFINAL=LY*1 
DO 57 1A=1 

DO 56 JA= 

X(JA)=XS*( 

56 Y( JA)=YS*( JA-1) 

CALL DRAWUA.X, 

57 M0D=2 
DO 58 IB=1,NUMDEC 
DO 58 JB=1,10 

J=AL0G10( YL0GMN)*1C.+10.*( IB-1) » 10. ♦JR 
YL0G=J6 

DO 1058 KB=1,IFINAL 
X(KB)=XS*(KB-1-IY) 

1058 Y( KB)=AL0G10(YLGG)+IR-1 

,ITITLE,XS,YS ,0.IY,2,2,LX,LY,0,LAST) 



, IFINAL 
,JF1NAL 
A-l-I Y) 



,M0D,0,LA(50)»ITITLE,XS,YS ,0,1Y,2,2,LX,LY.0,LAST) 



58 CALL DRAW(KB.X,Y,M0D,0,LA( J) ,I 
IF ( ( IPAR7) .NE. ( 1 ) ) GO TO 530 



59 IF( IWR ITE.LE.O) GO TO 200 
*>0r*4t«4>«««««««4<«*4cGRAPHICAL INFORMATION PRINTOUT <;eCT I ON************** 

WRITE(6,60) IRUN, NPARAM 

60 FORMAT( IHl.lOX, ' ♦♦♦♦♦♦♦♦♦*THE PRELIMINARY DATA FOR RUN NUMQER',12, 
£», parameter NUMBER', I2f» IS ♦♦♦♦♦♦♦***•-///) 

IF(NGAMMA.GE.I) write (6,61) ZGAMMA( IGAMMA) 

61 FORMAT( 10X,'THE FOLLOWING CALCULATIONS WILL BE RASED ON THE THIRD 
^VARIABLE PARAMETER VALUE , G AMMA = • , 1 PE 1 2 . 5 , / / ) 

WRITE (6,62) ( ITITLE( I ) ,1=1,12) 

62 FORMAT! lOX, 'THE GRAPH WILL RE T I TL ED : • , // , 1 OX , 6 A8 , / / , lOX , 6 AB , / ) 

IF( NPARAM. NE. 7) GO TO 80 

WRITE(6,83) 

WRITE(6, 63)nBMIN,DRMAX,XS 

63 FORMAT! lOX, 'MAGNITUDE, THE DEPENDENT VARIABLE, WILL BE PLOTTED ON 



AT', Flo. A, 



DB/OIVISIOM' ,/) 



_ ^c, inc 

CTHE ORDINATE FROM • ,F 10 .4 , • T0',F10.A, 

IF(IPHASE) 73,64,66 

64 WRITE(6,65) 

65 FORMAT! lOX, 'PHASE CURVES WILL RE PLOTTED ON MAGNITUDE GRAPH',/) 

GO TO 68 

66 WRITE!6,67) 

67 FORMAT! lOX, 'PHASE CURVES WILL RE PLOTTED ON SEPARATE GRAPH',/) 

68 IF! lOUAD.GT.O ) GO TO 70 
WRITE!6,69) 

69 FORMAT! lOX, 'FROM ZERO DEGREES THROUGH -90 TO -130 OEGRFFS',/) 

GO TO 72 

70 WRITE!6,71) 

71 FORMAT! ioX, 'FROM -90 DEGREES THROUGH -180 TO -270 DEGREES',/) 

72 WRITE!6,62) ! JTI TLE! I ) ,1 = 1 ,12) 

73 IF! IPOLAR.LE.O) GO TO 75 
WRITE!6,74) XS 

74 FORMAT! lOX, 'A POLAR PLOT OF MAGNITUDE VS. 

CITH MAGNITUDE SCALED AT',F10.5,' DECIBELS 

WRITE!6,62) ! KT I TL E ( I ) , I = 1 , 1 2 ) 

75 IF!M0D0RD-2) 127, 121, 121 

76 FORMAT! lOX ,A8, ' , THE DEPENDENT VARIABLE, _ _ _ 

80RDINATE F ROM • , 1 PE 1 0 . 3 , ' TO ' , 1 PE 10 . 3 , ' WITH A SCALE OF ' , IPFIO. 3, / ) 

77 FORMAT! 10X,A8, THE INDEPENDENT VARIABLE, WILL BE PLOTTED ON THE 
CABCISSA, FROM* ,lPE10.3, ' TO' , 1 PEIO . 3 , ' WITH A SCALE OF ' , 1 P E 1 0 . 3 , / ) 



PHASE WILL BE PLOTTED, 
PER HALF INCH',/) 



WILL BE PLOTTED ON THE 
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7fl FORMAT! lOX, Afl, • ^.ILL ME HELD CONSTANT AT ThF VALUE OF * , 1 PE 12 . 5 , / ) 

7*5 FORMAT! 10X,A8, • wILL ME GRAPHED AS ',13,' CUPVESf CONSISTING OF TH 

CE FOLLOWING VALUES: 16( 'CURVE*, I3f ' = ',lPF12.5,/n 
8C I F( IMAG.EQ.O) WRITE(6,B3) 

IF(IMAG.EJ.I) WRITE(6,8A) 

IF( IMAG.E0.2) WRITE(6,85) 

81 IF(M000P0-2) 82, SA,1000 

82 WRITE(6,76) JOKcR ( A) , 13ET AmN. BETAMX, YS 

93 FORMAT! lOX, 'ALL MAGNITUDE COMPUTATIONS WILL PE IN DECIBELS',/) 

9A FORMAT ! lOX, 'ALL MAGNITUDE COMPUTATIONS ARE ACTUAL MAGNITUDE',/) 

95 FORMAT! lOX, 'ALL MAGNITUDE COMPUTATIONS ARE MAGNITUDE SQUARED',/) 

GO TO !86,86,121,130,130,121,200,200),NPARAM 

86 WRIT6!6,77) JOK £ R ! 5 ) , ALPHMN , AL PH MX , XS 
IFINPARAM-1) 88,88,99 
88 wRITF!6,79) JOKER ! 7) .CONMAG 

9 0 UR I TE! 6, 79) JOKER! 6) ,NOMEGA, ! I ,ZOMEGA! I ) , 1 = 1, NOMEG A) 

IF! NPARAM-5)20C, 125,200 

94 WRITE!6,76) JPK E 8 ! 5 ) , BET Amn , BE T A MX , YS 

96 GO T0!97,97,l27,lll,lll,127,200,200) .NPARAM 

97 WRITE!6,77) JOK ER ( 4 ) , A LPHMN , AL PH MX , X S 
IF ! NPARAM- 1)89,88,99 

99 WRITE(6,78) JOkFR ( 6 ) .CCNOME 
101 wRITE!6,79)J0KER!7),NMAG,!I ,ZMAGS0! I ),I=1,NMAG) 

I F! NPARAM- 4) 200, 139,200 
111 WRITE!6,78) JOK E R ( 4 ) ,C ON AL P 

113 WRITE!6,114) 

114 FORMAT! lOX, ' THE COMBINED COEFFICIENTS OF THE NUMERATOR APE',/) 

115 DO 116 1=1, LOE 
A! I ) = HA! I) «-HAT! I ) uCUNST 
B! I ) = HBT! I ) 4 -hCT! I ) *CONST 
KA=I-1 

116 WRITEI6, 117) A! I ) ,P! I ) ,KA 

117 FORMAT! 20X, IH! ,E 10.4, 3H 
=RITE!6.136) 

00 119 t=l,MOE 

0! I )=HD! I ) '•HOT! I )*CONST 
E! I ) = HET! I ) ♦^HFT! I ) *CONST 
KA=I-1 

119 WRITE! 6. 117) U!I).E!I),KA 

120 GO TO !2C0,2OO,2C0,10l, PO , 200 , 1 39 ,200 ) , N^AR AM 

121 WRI TE! 6, 79) JOKER! 5) , NB ET A , ! I , Z BE T A ! I ) , 1=1 ,NBETA) 

IF!NPARAM-6) 124, 138, 111 

124 wRITE!6,78) JOK E P ! 7 ) , CONOME 

125 WRITE!6,77) JOK ER ! 7 ) , X maGmn , XMAG MX , X S 
GO TO 200 

127 WRI TE!6,79) J0KER!4),NALPhA, ! I , /ALPHA! I ) , 1=1, N ALPHA) 
IF!NPARAM-6)124, 138, 130 

130 uRITE!6.7d) JOKER ! 5 ) ,CCNBET 

131 FORMAT! lOX, 'BETA WILL BE HELD CONSTANT AT THF VALUE r)F*,F10.4,/) 

132 WPITEI6.114) 

DO 134 1=1, LOE 

A! I ) = HA! I ) ♦MbT! I )*CCNST 

01 I ) = HAT! I ) ♦HCT! I ) ♦CONST 
KA=I-1 

134 wP ITE! 6, 135) A! I ) ,B! I ) ,KA 

135 FORMAT! 20X, IH! ,E 10.4, 3H ♦ 

WRITF!6,136) 

136 FORMAT! lOX, • T he CCMBINEO COEFFICIENTS OF THE DENOMINATOR ARF-',/) 
on 137 1=1 , mpe 

0! I ) = H0! I ) ♦•HET! I I^CDNST 
E! I ) = hQT! I )*HFT! I )^CONST 
KA=I-1 



E10.4,16H ♦ BETA ) ♦ S** ! , I ? , 1 H ) / ) 



,E10.4,17H ♦ ALPHA ) ♦ $♦ * I , I 2 , 1 H ) / ) 



138 WRITE!6,78) JUKER ! 7) ,CONmaG 

139 wRI TE!6,140) 

140 FORMAT! lOX, 'PMEGA, THE INDEPENDENT VARIABLE, WUL BE PLOTTED ON THF 
f.ABC ISSA' , / ) 

IF ! NPARAM-6) 141,143,144 

141 IF! IGPAPH) 142, 142, 144 

142 WRITE!6,146) QMEGMN,0MFGMX,XS 
GO TO 200 

l‘*3 IF! ! IGRAPH.LT.O) .CH.nGRAPH.GT.in GO TO 142 

144 wRlTE(6,l45)NlJMDEC,YLnGMN 

145 FORMAT! lOX, *WI TH A SPAN OF', 12,' DEC AOF S , ST AR T I NG AT • , FIO. 4 , // , 30X 
Q ,****<.N0TE*****',//,10X,»F0R this log plot, it WILL BE NECESSARY 
CTO Shift THF origin to the lower right hand corner, ' ,//,10X, 'WITH 
CTHE ABCISSA running up the right hand FDGE, ANO THF ORDINATE INCRF 
CASING TOWARDS THE LEFT ALONG THE BOTTOM'//) 
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I<t6 FORMAK lOXf'FROM THE INITIAL VALUE OF • , 1 PE 10 . 3 . • TO THE FINAL VALUE 
6 OF*. 1PE10.3, • WITH A SCALE OF •»1PE10.3,' RAOlANS PER INCH',/) 
C4>«**4i4>4>4<4>i>*4i4>*««4<4>4>*MAIN COMPUTATIONAL SECTION DO LOOPS**************** 

200 IF( ( IWRITE.lt. 1) .OR. ( IWRITE.GT.900n GO TO 203 
WRITEtb, 201) IRUN,NPARAM,IWRITE 

201 FORMAT! 1H1,10X, **********THE OUTPUT DATA FOR RUN NUMBER • , 1 2, • , PAP 
6AMETER NUMBER', 12,' I S********** ',////, lOX, 12 , ' COMPUTATIONS WILL 

-R INTOUT' 



,1PE10.3,/) 



68E PERFORMED BETWEEN EACH PRINTOUT',///) 

IF(NGAMMA.GT.O) WRITE(6,202) ZGAMMA! IGAMMA) 

202 FORMAT! lOX, 'THE VALUE OF GAMMA, THE THIRD VARIABLE IS: 

203 DO 520 L0M=1,NLIMIT 
JCURVE=0 

DO 519 ICOEF=l,3,2 
ZSIGN=IC0EF-2 

IF! ! IwRITE.LT.D.OR.! IWRITE.GT.900) ) GO TO 208 
WRITE!6,205 )LOM 

205 FORMAT! /, lOX, 'COMPUTATIONAL RESULTS FOR CURVE NUMBER', 13) 

IF !NPARAM-7)2C6,2CB,20B 

206 WRITE!6,207 )ZSIGN 

207 FORMAT! 1H*,55X, • USING A QUADRATIC FORMULA COEFFICIENT 0F',FA.1,/) 

208 IP=0 
YMAGPK=0.0 
0MEGPK=0.0 

DU 51b INCRMT=1,900 

GO TO ! 209, 216,253,216,212 ,254,257, 1000) ,NPARAM 

209 OMEGA=ZnMEGA!LOM) 

211 C0NST = ALPHMN + ! INCRMT-1 )*!ALPHMX-ALPHMN)/B9<?.0 

212 IF!MGOORD-2)260, 213,213 

213 DO 214 1=1, LOE 

A! I )=HA! I ) *HAT! I )*CONST 

214 B! I )=HBT! I )*HCT! I)*CONST 
DO 215 I =1,M0E 

D! I ) = HD! I) +HDT! I )*CONST 

215 E! I)=HET! I)*HFT! I)*CONST 
GO TO 263 

216 IF! IMAG-1)217,218,219 

217 DUMMY 1 = 1 0.0** !ZMAGSO!LCM)/ 10. 0) 

GO TO 220 

2 IB DUMMY l = ZMAGSO!LOM)*ZMAGSO!LnM) 

GO TO 220 

219 DUMMY1=ZMAGS0!L0M) 

220 IFrNPARAM-2)211,211,212 

253 0MEGA=C0N0ME 

254 IF! M0D0R0-2)255. 256, 1000 

255 CONST=ZBETA! LOM) 

GO TO ?60 

256 CGNST=ZALPHA! LOM) 

GO TO 213 

257 IF!MODORD-2)25Q,258,25B 

258 BETA=ZBETA!LQM) 

CONST=CONALP 

GO TO 213 

259 BETA=ZALPHA!LOM) 

CONST=CONBET 

260 DO 261 1=1, LOE 

A! I )=HA! I)*H8T! I )*CONST 

261 B! I ) = HAT! I H-HCT! I )*CONST 

00 262 1=1, MOE 

D! I )=HD! I ) *HET! I )*CONST 

262 E! I )=HDT! I )*HFT! I )*CGNST 

263 GO TO 1 350,350,275, 290,264 , 29 1 , 294 , 1 000 ) , NP AR AM 

264 OMEGA=ZOMEGA!LOM) 

27 5 ZMAG=XMAGMN4-! INCRMT-1)*! XMAGMX-XMAGmN) /B99.0 
IF! IMAG-1)276,277,278 

276 DUMMY1=10.0**!ZMAG/10.0) 

GO TO 350 

277 DUMMY1=ZMAG*ZMAG 
GO TO 350 

27B 0UMMY1=7MAG 
GO TO 350 

290 IF! IGRAPH)292,292,294 

291 IF! ! IGRAPH.GE.O) .AND.! IGRAPH.LE. 1) ) GO TO 294 

292 0MEGA=0MEGMN+INCRMT*!0MEGMX-0M£GMN)/900.0 
GO TO 350 

294 WINC=1.0*! INCRMT-1 )*NUMDEC/ 8 99.0 
0MEGA=! YLOGMN*10.0**WINC)/10.0 

C *♦♦♦♦♦♦♦♦ ♦♦♦♦♦muMERATCR coefficient calculations*********************** 

350 SNB2=0.0 
SNB1=0.0 
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SNB0=0.0 
SUMN=0.0 
N00UN=2*N0RN 
00 373 IL0QP=1,^ 

DO 373 NP=2,L0E 
DO 373 M0=1,N0PN 
MP=M0-1 
N0=NP-1 

IF( (NP+Mgj-NnOUN) 36 1,361, 373 
3bl IF( (MP+NQ )-N0RN ) 362, 362, 373 

362 IF( NP-MP ) 373, 373, 363 

363 IF(HP) 366, 366,366 
366 C=1.0 

GO TO 367 

366 C=2.0 

367 GO TO < 36d, 369,370,371) , ILOOP 

360 CSN82 = (C*B<NP+HP)*B«NP-MP) )*0MEGA**(2*N0»*<-1.0)**(MQ+1) 

SN82 = SN32 ♦ CSNP2 
GO TO 373 

369 CSNBl = (C*A<NP + MP)C!3( NP-«PJ ♦C*A(NP-MP) *R(NP+MP) )*riMFr,A** (2*MQ)*(- 

11.0) **( Mo+n 

SN81 = SNfll ♦ CSNBl 
GO TO 373 

370 CSNBO = (C*A(NP^HP)*A(NP-MP) )*0MEGA**(2*Ng)*(-1.0l**<MQ4-l) 

SN80 = SNHO ♦ CSNPC 

GO TO 373 

371 COEFN = <C*A<NP + '<P)*A(NP-HP)4C*A(NP>MP)*R(NP-MP)*BFTA4-C*A(NP-mp)*b 
KNP+'-P ) *BETA+C*H<NP+MP)*B(NP-MP) *BETA*BET A)*OMEGA-t'*( 2*N0 ) * ( -1 . 0 ) ** 
2(M0+1 ) 

SUMN = SUMN ♦ COEFN 
373 CONTINUE 

SNB2 = SN82 + B( 1)*B( 1 ) 

SNBl = SN81 ♦ 2,0»A( 1) *8( 1 ) 

SN80=SNP0 + A( 1 )*A( 1 ) 

SUMN = SlJMM ♦ A( 1)*A(1 )>2.0*Am*B(l )*0FTA^8(1 J=»B(1)*BETA'<'BETA 
C************pENnylNATQP COEFFICIENT CALCULATIONS **♦*♦♦**•*♦*♦**♦•*•***<' 
S082=0.0 
SOBl^O. 

SOBO=0. 

SUMO^O. 

N00UD=2»N'JR0 
00 385 ILOOP=l,6 
00 385 NP=2,MQE 
DO 385 mq=i,noRO 
MP=MQ-1 
NQ=NP- 1 

IF( (NP + ^'Q)-NU0lir)) 3 75,375,385 
37 5 IF( (MP♦^0) -6'ORO) 376, )76,3«5 

376 IF(nP-MP) 385,’85,377 

377 IF(mp)378,378,379 
370 C=1.0 

GO TO 3PC 

379 C=2.0 

380 GO TC (301,332,303,386), ILOOP 

381 CSDB2 = C*E(NP-*-0F)*E(NP-MP) ♦0MEGA**(2*N0)*(-1.0)*<=( MO+l) 

SDB2 = SDB2 ♦ CS0P2 

GO TO 385 

382 CSDBl = (C*J(NP + HP)<=E(NP-MP)*C*0(NP-MP)*E(NP + mP) )*rMEGA*=*( 2*NQ)*(- 

11.0) **(MQ*1) 

S091 = SD3l ♦ CSOBl 
GO TO 385 

333 CSD80 = (C*D(NP + ''P)*0(\P-MP) )*0MEGA**(2*N0)*(-1.0)**(m0+1 ) 

SDBO = S030 ♦ CSDPO 
GO TO 305 

386 COEFD = (C*l)(NP*MP)*a(NP-MP)+C*0(NP+'^P)*E(ND-MP)*8ETA^C*D(NP-MP)tF- 
1(NP>MP)>^8PTA>C*E(NP + *'P)*E(NP-MP)*8FTA*86TA) kOmPGA**! ?♦ NO )*(-!. 0)** 
2(M0 + 1 ) 

COEFD 



SUMD = SUMO 
385 CuNTINUF 

SUMO = SUMO 
S082 = SDB2 
S061 = S081 
SDBO = 5D80 



IF ( NPARAM-7) 613, 
***********#*^^'^- 

386 T0E=L0E/2.0 
PNEVEN=0. 

PDEVEN' = 0. 
PNQ00=0. 



U( 1)*D( 1 )>-2.0*n( 1)*E( 1 )*8ETA + F( 1 )*F( 1 )*8FTA*PETA 
E( 1)*E( 1 ) 

2.0*D( 1)*E( 1) 

o( i)*on ) 



7 PHASE computational SECT I0N*«*<‘**** + *«>!' 
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P00D0=0. 

J0E=L0E/2 

IF ( JOE-TOE ) 387, 3 88, 388 
J0N=(L0E4-1)/2 
MON^JON 
GO TO 389 
J0N=(L0E4^2)/2 
MON=JON-l 

IF(NORN-l) 393,389,389 
DO 390 1=2, JON 

PNODO = PNGOO*-( A( 2*1-2 ) + B( 2* I -2 ) *BE T A J *QmeGA** ( 2* I -3 ) * ( -1 . 0 ) ** ( I ) 
IF(NORN-l) 393,393,391 
DO 392 1=2, MON 

PNEVEN=PNEVEN<-< A(2*I-1 )♦ R ( 2* I - 1 ) *BE T A ) *0mEGA** ( 2* 1-2 ) * I - 1 . )**n-l) 
PNEVEN=PNEVEN+A( 1)+8(1)*BETA 
T0E=M0E/2.0 
J0E=M0E/2 

IF! JOE-TOE ) 3 9A, 395,395 
J0N=(M0E^l)/2 
MON=JON 
GO TO 396 
JON=(MOE+2)/2 
M0N=J0N-1 

IF(N0RD-1 1400,396,396 
DO 397 1=2, JON 

POOOO=PDODO<-(D( 2*1-2 ) + E( 2*1 - 2 ) *BE T A ) *OMEG A* *( 2 * I -3 ) * ( -1 . 0 ) ** ( I ) 
IF(NORO-l) 400,400, 398 
DO 399 1=2, MON 

PDEVEN=P0EVEN*-(D(2*I-1 »*E(2*I-l)*BETA)*0MEGA* + ( 2*1-2 )* ( -1 . ) ** ( I-l ) 
PDEVEN=PDcVEN+D< 1 H-E( 1 )*BETA 

TANG=( PDEVEN*PNODO-PNEVEN*PDODD) / ( PNF V FN*POE VE N+PNODO* PDOOO ) 
RAD=ATAN< TANGI 
SINE=SIN(RAO) 

IF( TANG)^02,402,401 
SINE=-SINE 

YMAGS0=ABS( SUMN/SUMO) 

YMAG=SQRT< YMAGSO ) 

YMAG0B=20.0*.4342945*AL0G( YMAG) 

YMAGAB=ABS( YMAGDB) 

IF( YMAGPK-YMAGDB 1403,404,404 

ymagpk=ymagob 
omegpk=omega 

IF ( DBM AX- YMAGDB ) 405,406,406 
YMAG0B=DBMAX 

IF < YMAGDB-OBMIN >407,408,408 
YMAG0«=DBMIN 

IF( ( XS*e.O)-YMAGAB >409,^10,4 10 
YMAGAB=8.0*XS 
GO TO 426 

♦ ♦♦♦♦♦♦**«r*OUAPRATIC FORMULA SOLUTION OF DEPENDENT V AP I ABL E , P AR A»< 
AC0=0UMMY1*SDB2-SNB2 
BC0=0UMMY1*SDB1-SNP1 
CC0=0UMMY1*SDB0-SNB0 
IF( ARS(ACO) .LT.UOE-35) GO TO 425 
OPP=8CO*BCO-4.0*ACO*CCO 
IF( DPP.lt. 0.0» GO TO 425 
DS0R=ZSIGN*S0RT( DPP) 

BETA=(DSOR-BCO)/< 2.0*AC0) 

IF ( BE TA-BETAMN >425,426, 424 
424 IF(BETA-BETAMX>426,426,425 
<•25 IF! IP-l > 500,505.507 

<•26 IF( ( IWRITE.lt. 1) .OR. (IWR1TE.GT.900> > GO TO 450 
C***««**«***<c«4c*4(«<t*C0MPUTE0 DATA PRINTOUT SECT I ON********************** 
ISTEP=ISTEP+1 
IF( IP>^28,428,<.27 
<•27 IF( ISTEP-IWRITF >450,430,430 

428 WRITEI6,429>LUM 

429 FORMAT! /,10X, • ALPHA BETA MAGNITUDE OmFGA 

■ TANGENT ; CURVE NUMBER • , I 2 , / > 



387 



388 



391 

392 

393 



394 



395 



398 

399 
40 0 



40 3 

404 

405 

406 
40 7 
408 
40 9 
410 

413 

416 



430 


IF( 


: MOOORD.I 




GO 


TO 


(433 


431 


GO 


TO 


(432 


432 


WR ] 


[TE( 


16, 6> 




GO 


TO 


448 


433 


WR] 


[TE( 


16, 6> 




GO 


TO 


448 


434 


WR 1 


[TEI 


16, 6> 




GO 


TO 


448 


435 


WR I 


[TEI 


16,6) 



34, 445 ,4 3H ,440,4 4^' , 444, 10001 , N 
C0NST,BETA,C0NMAG,0MFGA 

BETA, CONS T,C0NMAG,0MFGA 
CONST,BETA, ZMAGSQI LOM> , OMEGA 
BETA,CONST, ZMAGSO(LOM) , OMEGA 
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Rt T A, CONS TfZMAG , omega 
BETA.CDNRET »ZMAGSO( LCM) , OMEGA 
conalp,8eta,zmagsq(lom) , omega 
RETA,CON6ET,ZMAGtOMEGA 

conalp,beta,zmag,omega 

BETA»CONST,CCNMAG,nMEGA 
CONST ,8ETA,C0NMAG,CMFGA 
ZALPHAdC'^) ,C0NBET,YMAG08, omega, TANG 
CONAL P, Z BE TA( COM ) ,YMAGDB, OMEGA, TANG 
CON ST, bet A, ZMAG, OMEGA 



GO TO <,<,8 
^36 WRITE(6,6) 

GO TO 

-V37 WRITE(6,6) 

GO TO ^^8 
^38 n«ITE(6,6) 

GO TO 448 
43Q WRITt(6,ft) 

GO TO 448 

440 WRITE(6,6) 

GO TO 448 

441 WKITE(6,6) 

GO TO 448 

442 WRITE(6,6) 

GO TO 448 

‘,43 KRITE(6,6) 

GO TO 448 
4,,4 WRITE(6,6) 

GO TO 448 
445 k,RITE(6,6) 

448 ISTEP=0 

C*«****«**«*»***«*«*GRAPHI CAL ARRAY STORAGE SECT I nN**^****»* f****** 

450 IP=IP*1 

Y( IP)= BETA 

^52 GO TO ( 453, 453, 460,470, 460, 470, -VTS , 1000 ) ,NPARAM 
453 X( IP)= CONST 
GO TO 500 
460 X(IP)=ZmaG 
GO TO SCO 
470 X( IP)=iJMEGA 

IF( NP ARAM- 6 >471,472, 475 
*♦71 IF ( IGRAPH) 500,500,474 

472 IF( ( IG«APH.LT.0> .OP. (IGRAPH.GT.in GO TO 500 

474 X( IP)=LX-8ETA 
GO TO 475 

475 X( IP)=-YM4GDB 

U( LOM, 1P)=IJ(L0M, IP)^Y‘^AG 

V(LOM, IP)-( ALOGIOI CMEGA/ YLOG MN ) *LY /NUMDEC) ♦YS 
S( LCM, 1P) = S1NE*YMAGA8/( XS*TANG) 

T( L0« , IP )=SINF*YKAGAB/XS 
IC(L0«)=IP 

IF( I0UAD)476,479,478 

476 tolLOM, IP)=-4.0«RAD/1.5708 
IF(RA0)477,47S,475 

477 W(LCM, IP) = 9.0-V,(LnM, IP) 

GO TO 479 

47 8 4(LUM, IP) = XS*(LX/2.-1Y)-RAD<=XS*LX/3.1414 
IF( ( IP) .LE.( 1) )G0 TO ‘*79 
IF ( ( YMAGOB) .NE. ( CHMI N) ) GO TO 479 
IF( ( W(LOM, IP) ) .GE. ( W(LQm, IP-1) ) )G0 TO 47^ 

WILOM, 1P)=W(LCM, IP-1 ) 

V(LO*<, IP)=V(LOM, IP-1) 

4 79 Y( IP) = ( AL0G10(P‘^fGA/YICGmN)*LY/NUM0EC)*YS 

500 I F( INCR»^T.LT.900 ) GO TC 518 

IF( lP-1) 501,505,507 

501 IFIJCURVF.GT.O ) GO TO 503 
WRITE (6, 502) COM 

502 FOR'^ATI lOX, *NC GRAPH POINTS WFRF FOUND FOR CURVE NUMRFR *,12,/) 

GO TO 518 

503 wRITE(6, 504)LIIM 

504 FORMAT! lOX, 'NO FURTHER POINTS WERE FOUND FOR CURVF N).',I3,/) 

GO TO 518 

505 wRITE(6,506) 

506 FORMAT! lOX, »ONLY ONE POINT WAS FOUND, AND WILL NOT PF PLOTTED',/) 
GO TO 516 

507 J = Ln»^^50 
JCURVE=1 

IFIIGRAPH.LE.-2) GC TO 516 

IF! !NPAPAM.E0. 7) .AN'O. ! 1PAP7.lt. NPAR7) )G0 TO 513 
I F! ! ZSIGN.GT.0.0) . AND. ( JCURVE.GT.O) ) J=50 
wP1TE!6,511)L0v 

511 FORMAT!/, lOX, 'THE PREVIOUS POINTS WIIL BE PLOTTED AS CURVE ',!?,/) 
IF!NPARAM-7)^4,512,514 

512 WR ITE!6,513) YMAGPK ,0«FGPK 

513 FORMAT! / ,10X ,' THE MAXIMUM VALUE OF maGNITDOE IS*,F10.4,* LOCATED A 
CT OMEGA EQUAL T0',F10.4,/) 

J*50 

514 call draw! IP,X, Y,MC n,0,LA! J) ,IT1TLE,XS,YS,1X,IY,2,2,LX,LY,],LAST) 
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M0D=2 

516 IP=0 
5ie CONTINUE 

IF(NPARAM.EQ.7) GO TO 520 

519 CONTINUE 

520 CONTINUE 

IF( ( IPAR7.GE.NPAR7).CR.( NPARAM.NE.7) ) GO TO 529 
C*»»****gAMMA £ PARAM-7 SUPERPOSITION ; SEE INPUT DATA SYMBOL NGAMMA+*** 

517 IF ( (NGAMMA.LE .0) .OR. ( I GAMMA.GT.NGAMMA) ) GAMMA=0.0 
IF( IPAR7.GT. 1 )G0 TO 523 

DO 521 1=1, LOE 
OLDHAI n=HA( I ) 

OLDHATI I )=HAT( I ) 

OLOHBTI I )=HBT( I » 

521 OLOHCTI I )=HCT( I ) 

DO 522 1=1, MOE 
OLDHO ( I )=H0 ( I ) 

OLDHOTI n=HOT( I » 

OLDHETI n=HET( I ) 

522 OLDHFTI I ) = HFT( I ) 

523 IPAR7=IPAR7*1 
IF(NGAMMA.GT.O) GO TO 526 



REA0(5,5) (HA(I). 
READ(5,5) (HAT(I) 
REA0(5,5) (HBT(I) 
REA0(5,5) (HCT(I) 
HEA0<5,5) (HO (I) 
REA0(5,5) (HOT(I) 
R£AD(5,5) (HET(I) 
READ! 5,5) (HFT( I ) 



a, LOE ) 
1=1, LOE) 
1=1 , LOE) 
1=1, LOE) 
1=1, MOE) 
1=1 ,HOE) 
1=1 , MOE) 
1=1, MOE) 



KtAUI6,!>J l=l,"1Utl 

IF( IPAR7.lt. NPAR7) GO TO 200 
00 525 L0M=1,NLIMIT 
IP=IC(LOM) 

00 524 1=1, IP 
X( I )=-U( LOM. I ) 

524 Y( I ) = V(L0M,1 ) 

J=L0M+1 

52 5 CALL DRAW( I P , X , Y , MOO ,Q , L A ( J ) , I T I TLE , XS , YS , I X , IY,2,2,LX,LY,1,LAST) 

526 00 527 1=1, LOE 

HA (I)=OLDHA ( I )*GAMMA»GAMHA (I) 

HATH )=CLOHAT( I ) +GAMMA+GAMHAT ( I ) 

HBT( I )=nL0H8T( I ) +GAMMA»GAMH8T( I ) 

527 HCT( I )=OLOHCT( I ) +GAMMA*G AMHC T (I ) 

00 528 1=1, MOE 

HO (I)=OLDHO ( I )+GAMMA*GAMHO (I) 

HDT( I )-^OLOHOT( I )+GAMMA*GAMH0T( I ) 

HET( I )=OLOHET( I ) +GAMMA*G AMHET( I ) 

528 HFT( I )=OLOHFT( I ) +G AMMA*G AMHF T( I ) 

IF( IGAMMA.LE.NGAMMA ) GO TO 26 

529 IFMNPARAM.NE.7) .OR. ( IGRAPH.LE.-2) ) GO TO 999 
E******»»*****#*****««»*«*+p ARAM 7 phase curves*************************' 

IPAR7=2 

IF( IPHASE) 537,532,999 

530 00 531 1=1,12 

531 ITITL£( I )=JTITLE( I ) 

532 00 536 L0M=1,NLIMIT 
IO=IC(LOM) 

00 533 1=1, IP 
X( I ) = W(LOM, I ) 

533 Y( I )=V(LOM, I ) 

536 CALL ORAW( I P , X , Y , MOO , 0 , L A ( 50 ) ,JTITLE,XS,1 .,0,IY,2,?,LX,LY,1,LAST) 
IPAR7=1 

537 IF(( IPOLAR).LE.(O) ) GO TO 999 

Q^i:^>tf^i*^:t:*******************P ^ POLAR PLOT CURVES******************** 

M00=1 

00 542 1=1,8 
RAOIUS=I/2.0 
00 539 J=l,360 
R A0IAN=6. 283 18/360* J 
X( J ) = RA0IUS*C0S( RADI AN) 

539 Y( J)=RAOIUS*SIN(RAOIAN) 

CALL ORAW( J,X, Y,MO0,0,LA( 1),KTITLE,1.,1.,4,4,2,2,8,8,0,LAST) 

542 M00=2 

00 543 J=10,180,10 
X( 1)=X( J) 

X( 2)=-X( 1 ) 

Y( 1)=Y( J) 

Y( 2)=-Y( 1) 

543 CALL 0RAW(2,X,Y,2,0,LA(1) ,KTITLF,1.,1. ,4,4,2,2,R,S,0,LAST) 
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mTLE,1.0,1.0,0,0,2,2,LX,LY,l,L4ST) 



Of] 54 5 L0M = 1,NL1 ^«IT 

iP=iC(LrM) 

DO 544 J=l,IP 
X( J ) = S( LOM, J » 

544 Y( J»=T(LOM, J) 

L=LnMfl 

545 call 0PAW( ip ,X,Y, 2,0,LA( L» ,KTITLE, 1. , I . ,4,4,2, 2 t5,8,0,LAST » 
999 IF ( (^<a0) .NE. (2) ) GC TO 1000 

XI I ) = 0.0 
Yin = 0.0 
X(2» = (LX-IY»*XS 
Y(2» = 0.0 

CALL DKAW ( 2 ,X, Y,3,0,LA( 50» 

*< 00=1 

IF( I IPAR7) .NE.( I ) » GC TO 55 
1000 CONTINUE 
STOP 
tNO 

C*********SAMPlE data 
I 

4 12 4 

0 3 1 fl 

0 10 -1 I 

C.036 ■“ ‘ 

0.007 
0.0 
0.01 



0.00047 
> 0.001 
- 0.001 
0.0 
-3.0 
1.0 
45. 

100000.0 

( cfifo') 

100000. 500.0 

100000.0 



.0 

4.0 

50. 



12.0 

C.O 

0.G0047 

0.02 

3.C 

IC.O 

55. 



10 , 

C.05 

0.002 

.03 

6.0 

12.0 

60. 



-3.0 

0.003 

0,036 

0.0 

15.0 

70. 



^ 1.0 

0.004 

0.04 

12.0 

20.0 

BO. 



0,0 
0.006 
0.045 
15.0 
] 0 .0 



60. C 

lOOOCO.O 



GLAVIS, PAKAM 1, 3KD QHLFR T/A F.B., OMEGA ClJRVS 
OROINATE= ALPHA , AnC I S S A= BETA , CONMAGM TU0E = 1 20B 
GLAVIS, PARAM Z, 3K0 ORDER T/A F.P., MAG CURVES 
ORCINATE=ALPHA , ABC I SSA = 6E TA , CONCMEG A = 1 OR AD /SEC 
GLAVIS, PARAM 3RD ORDER T/A F.B., MAG CURVES 
QRDINATE = ALRHA , ABC! SSA=CmeGA, C0NBETA=0 .0 36 

?P0 ORDER T/A ^.B., BETA CURVES 

AHCISSA=CMEGA, CONMAGN I TUOE= 12 D« 
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LOGIC OUTLINE OF PROGRAM PARAMS 



1. Read number of runs submitted. 

2. Set run number equal to 1. 

3. Read input data for appropriate run. 

4. Print input data for appropriate run, if requested. 

5. Calculate graph scales and location of origin. 

6. Set NPARAM equal to the first PARAM value. 

7. Select initial values for appropriate PARAM. 

8. Print graphical information for appropriate PARAM, if requested. 

9. Construct logarithmic grid if required. 

10. Set curve number equal to 1. 

11. Select curve value for appropriate curve number. 

12. Set abcissa variable to minimum allowable value. 

13. Set number of graph points equal to 0. 

14. Compute numerator and denominator polynomial coefficients. 

15. Check for possible division by 0. 

16. Calculate ordinate value, using the quadratic formula for PARAM-1 
through PARAM- 6. 

17. Check if calculated ordinate value is within graph limits. If yes, 
continue; if no, go to statement 21. 

18. Increase number of points by 1. 

19. Print computed values. 

20. Store values necessary plotting. 

21. Increase abcissa variable by pre-computed increment. 

22. Check if abcissa variable exceeds maximum graph limit. If yes, 
continue; if no, go to statement 14. 
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23. Check if any points were found within graph limits. If none, 
so state and continue without plotting. If one, so state 

and continue without plotting. If more than one, continue and 
plot the graph points. 

24. Check if both positive and negative coefficients were used in the 
quadratic formula solution where applicable. If yes, continue; 
if no, change sign and return to statement 12. 

25. Increase curve number by 1. 

26. Check if curve number exceeds the maximum number of curves stipulated 
in input data. If yes, continue; if no, return to statement 11. 

27. Draw graph and plot computed curve data on same. 

28. Increase parameter number by 1. 

29. Check if parameter number exceeds maximum. If yes, continue; if no, 
return to statement 7. 

30. Increase run number by 1. 

31. Check if run number exceeds maximum. If no, return to statement 3; 
if yes, stop. 
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FLOW CHART FOR PARAMETER SIX 




178 




179 



1 1 { 






A A A 




180 




181 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 20 

Cameron Station 

Alexandria, Virginia, 22314 

2. Library 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Prof. G. J, Thaler 5 

Department of Electrical Engineering 

Naval Postgraduate School 
Monterey, California 93940 

4. Naval Ship Systems Command (Code 2052) 1 

Department of the Navy 

Washington, D. C. 20360 

5. Dr. A. G. Thompson 1 



Department of Mechanical Engineering 
University of Adelaide 
Adelaide, Australia 

6. Prof. D. D. Siljak 1 

Department of Electrical Engineering 

University of Santa Clara 
Santa Clara, California 95053 

7. LT George 0. Glavis, USN 1 

Engineering Officer 

USS HENRY B. WILSON (DDG-7) 

% FPO San Francisco, 96601 

8. H. H. Choe 1 

Division of Control Engineering 

University of Saskatchewan 
Saskatoon, Saskatchewan, Canada 

9. Dr. D. R. Towill 1 

Welch College of Advanced Technology 

Cathays Park, Cardiff, Wales 



182 



Unclassified 

S«'«~untv Classification 



DOCUMENT CONTROL DATA - R & D 

^ tassification of title, body of .ihsfr.irf and indexing annotation must be entered when the overall report is classified) 



oN* No ACT viTY Corporate author) 

Naval Postgraduate School 
Monterey, California 93940 



aa. REPORT SECURITY CLASSIFICATION 

Unclassified 



Zb. GROUP 



Frequency Response in the Parameter Plane 



Au TihOR'Si I First name, middle initial, last name) 

Glavis, George Oscar 



6 REPORT DATE 

June 1968 



7a. TOTAL NO. OF PAGES 

183 



7b. NO. OF REFS 
12 



Ba. CONTRACT OR GRANT NO. 



6. PROJEC T NO. 



9a. ORIGINATOR’S REPORT NUMBER(S) 



N/A 



N/A 



10 O' STRI BU TION ST A TEMEN T 
























11 Supplementary notes 


12. SPONSORING MILITARY ACTIVITY 

Naval Postgraduate School 
Monterey, California 93940 



3 abstract 



Algebraic Methods are employed for frequency response studies of physical 
systems. The transfer function of each system's mathematical model is composed 
of at least two independent variables embedded in the coefficients. These, 
together with magnitude, angular frequency, and phase angle are the variable 
parameters to be studied. 

Computer programs are developed for graphically representing these 
variable parameters. Interpretation of these results is made to evaluate the 
utility of such methods in design and analysis studies. 



DD .ir.,1473 

S/N 01 01 -807-681 1 



(PAGE 1 ) 



183 



Unclassified 



Security Classification 




Unclassified 

Soi'uritv Classification 



< E Y wo RDS 



Frequency Response 

Mitrovic 

Sil jak 

Parameter Plane 
Algebraic Methods 
Control Systems 
Compensation Techniques 






ROLE W T 



ROLE W 



-.fix 



DD ,T:.,1473 (BACK) 



184 



Unclassified 



Security Classification 



